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Abstract 


According to the magneto-hydrodynamic theory of sunspots there should be a correla- 
tion between the spots on one hemisphere during a certain cycle and spots on the opposite 
hemisphere during the following cycle. This has been statistically confirmed by Galvenius 
and Wold, and by Whittle. It is pointed out that individual sunspot groups may be 
reproduced at the opposite hemisphere during the following cycle. An analysis of the 
sunspots between — 4° and + 4° seems to support this conclusion, and a number of sun- 
spots are “‘identified’’ as being produced by the same original disturbance in the solar core. 


I. Magneto-hydrodynamic theory of sunspots 


In the magneto-hydrodynamic theory of 
sunspots a rather complicated mechanism 
must be postulated by which the sunspot 
producing whirl-rings should originate in the 
solar core!. The relatively small latitude 
displacement of sunspots shows that the plane 
of the rings must be parallel to the magnetic 
field. The difference in magnetic polarity 
between the northern and the southern hemi- 
spheres indicates that the spots occurring during 
a certain cycle at the northern hemisphere 
cannot have originated at the same time as 
the spots of the southern hemisphere during 
the same cycle. This leads to the conclusion 
that the magnetic lines of force must pass two 


1 Alfvén (1948). Summarized in Cosmical Electrody- 
‘namics, Oxford 1950. Other theories of sunspots have re- 
cently been proposed by Walén and Schwartzschild. Wa- 
léns earlier papers (Ark. f. mat., astr. o. fysik 30 A, No 15; 
31 B No 3 1944) are reviewed in Cosmical Electrody- 
namics. The later development of his ideas, with which 
the present author is not in agreement, are published in 
Ark. f. mat., astr. o. fysik 33 A No 18, 1946, and On the 
vibratory rotation of the sun, Lindstähls, Stockholm 
1949. Schwarzschild’s contribution was presented at the 
Solar physics conference in Florence 1952. 
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activity regions in the solar core, one to the 
north and one to the south of the equatorial 
plane (cf. Fig. 1). The production of a dis- 
turbance e.g. in the northern region gives rise 
to magneto-hydrodynamic waves which about 
40 years later appear on the surface of the 
northern hemisphere and produce sunspots. 
The same disturbance, however, also proceeds 
southwards and passes the southern activity 
region, where it becomes amplified, and 
finally reaches the surface of the southern 
hemisphere where it produces sunspots. These 
spots appear about 11 years after the spots due 
to the same original disturbance have appeared 
at the northern hemisphere. 


2. Statistical confirmation 


Although a mechanism of about this kind 
seems to be an inevitable consequence of the 
magneto-hydrodynamic theory, it is obvious 
that such a complicated mechanism could not 
be accepted without very strong empirical 
confirmation. The first direct support of the 
theory was obtained by an analysis of the 
latitude distribution of sunspots on the north- 
ern hemisphere, compared with the corre- 
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Fig. 1. Magnetic field within the sun. Meridional pro- 
jection. 


sponding distribution of the southern hemi- 
sphere 11 years later(ALFVÉN, 1948). According 
to the theory these patterns should be similar 
because the sunspots of one hemisphere have 
the same origin as sunspots of the other 
hemisphere 11 years later. A statistical in- 
vestigation by GALVENIUS and WoLD (1948) 
showed that the similarity of the patterns was 
statistically significant. 

Nextconfirmation was obtained by WHITTLE 
(1954) who investigated the correlation be- 
tween the time series of spots in the latitude 
belts 16°—20° at the northern and at the 
southern hemisphere. He was able to demon- 
strate in a very convincing way that there 
exists a correlation between the two hemi- 


spheres of the kind required by the theory. 


3. Possible methods for spot correlation 


The statistical correlation between the two 
hemispheres which consequently has been 
demonstrated in two independent ways, gives 
a strong support to the general lines of the 
magneto-hydrodynamic theory. There could 
no longer be any doubt that the spot activity 
on one hemisphere is coupled with the activity 
on the other hemisphere during next cycle. 
So far it has only been demonstrated that there 
is a statistical connection and this may be due 
to some factor which in a general way regulates 
the activity. It is also possible, however, that 
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there is amuch more direct connection between 
the spots of the opposite hemispheres. The 
most extreme case would be that each individual 
spot on one hemisphere is directly connected with 
a certain spot on the other hemisphere one cycle 
later. The purpose of the present paper is to 
investigate whether such a connection exists. 
The work was started as soon as WHITTLE had 
got his positive results and inspired by this. 
Suppose that a disturbance in e.g. the 
northern activity region of the solar core 
produces a whirl ring. This moves out along 
the magnetic lines of force to the surface and 
appears as a spot An in the northern hemi- 
sphere. The same disturbance also produces a 
wave which moves through the core, passes 
the southern activity region, where it may be 
in part amplified and in part reflected, and 
finally reaches the surface and produces a spot 
As on the southern hemisphere one cycle 
later. The problem is: how could we identify 
the spots An and As as being caused by the 
same original disturbance. There are two 
different possibilities: 
1. Identification by position. The two spots should 
occur at the two intersections with the solar 
surface of the same magnetic line of force 
with a time difference corresponding to the 
time of travel through the core (which is of 
the order 11 years). If we could find that point 
of the southern hemisphere which is inter- 
connected by a line of force through the 
interior with a certain point of the northern 
hemisphere we could study whether spots at 
these two points preferably occur with a 
certain time interval (of the order 11 years). 


2. Identification by structure. The two spots which 
originate from the same disturbance in the 
solar core may have some similarity in struct- 
ure. For example, their size, life time, or proper 
motion, may be correlated in some wav. 


4. Regularity of solar magnetic field 


The identification by position would be easy 
if we knew the shape of the magnetic field 
inside the sun. This is not the case: We could, 
however, make a reasonable guess about it 
under the condition that it is rather regular. 
On the other hand if the solar field is very 
irregular, we have little hope to succeed. The 
first question to discuss is whether we should 
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expect the field in the sun’s interior to be 
regular or not. 

At first sight there seems to be an over- 
whelming probability that the solar field is 
very irregular. Zeeman effect measurements 
of the surface field indicate an irregular 
structure and strong variations from one day 
to next. (Whether Zeeman effect measure- 
ments of weak fields are reliable at all is still 
an open question (ALFVEN, 1951). Further the 
structure of the corona, which probably is 
governed to a considerable extent by the 
magnetic field, indicates a variable and irregular 
structure. These arguments are not decisive, 
however, because what is important is the 
regularities of the “general” field in which the 
magneto-hydrodynamic waves proceed. This 
could be thought of as the time average of the 
actual field. It is superimposed by a variable 
field which the waves carry with them. Hence 
the irregularities and rapid changes of the 
actual magnetic field are quite consistent with 
a regular general magnetic field. 

Further, the magnetic lines of force pass at 
least two layers, which according to current 
solar models should be turbulent, viz. in the 
core and in the photosphere. This seems to 
exclude the possibility of a regular field. Again 
it must be observed that the field of interest is 
the time average. A turbulence should be 
described as a series of waves travelling in the 
general field, but it does not contradict the 
assumption of a regular general field. 

Having thus found that there is no decisive 
argument against a regular field we shall now 
examine the arguments in favour of it. Ac- 
cording to Ferraro’s theorem all parts of a 
certain line of force must in average rotate 
with equal angular velocity. As the sun has a 
non-uniform rotation it is likely that this 
brings with it that the rotational axis and the 
magnetic axis must coincide (although this has 
not been rigidly proved). If the solar magnetic 
field is of Cowling’s decaying type it is quite 
reasonable that it is very regular because at 
the decay all irregularities should disappear 
rapidly. LUNDQUIST has discussed the shape of 
“magneto-hydrostatic fields” which can exist 
in a liquid medium without producing a 
hydrodynamic motion (LUNDQUIST, 1951). 
Such fields are in general “twisted”. If the sun’s 
general field is of this type it is quite reasonable 
that it should be regular but “twisted”, i. e. 
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Fig. 2. Twisted magnetic field. 


the field has a longitudinal component, so that 
a magnetic line of force intersects the surface 
at two points which do not have the same 
longitude. In fact, there is no reason to suppose 
that the longitude difference between the 
northern and the southern intersections with 
the solar surface should be zero. It is important 
to observe that Fig. ı shows only the projection 
of a line of force on the meridional plane. The 
real shape is obtained if we rotate this plane 
around the axis when we proceed along the 
line of force. The spiral structure of a line of 
force is shown in Fig. 2. 

Cowling has shown that a general magnetic 
field cannot be sustained if it is regular and 
symmetrical. The models by ELSASSER (1946, 
1947), BULLARD (1949), and others are rather 
complicated and irregular. A regular field 
which does not decay systematically (counted 
over very long periods) is compatible with a 
model proposed some time ago, which assumes 
a decaying field, at intervals reinforced by a 
‘convulsion’ (ALFVÉN, 1950). 

The irregularities and disturbances of the 
general field according to the mentioned 
models need not necessarily affect the lines of 
force which intersect the solar surface. They 
can very well be confined to the lines of force 
which altogether are hidden in the sun’s 
interior. 

Hence theoretical considerations rather sup- 
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port the view that the general field may be 
regular and symmetrical but possibly twisted. 
The conclusions are not decisive, however. 


5. Empirical arguments for a regular field 


There are quite a few empirical arguments 
for a regular field. During a certain cycle 
bipolar sunspots have a certain polarity at the 
northern hemisphere and the opposite polarity 
at the southern hemisphere. The line of 
division is the equator. Only very seldom does 
a spot of the “northern” polarity proceed as 
much as one degree to the south of the equator, 
and vice versa(WALDMEIER). Even independent 
of the magneto-hydrodynamic theory of 
sunspots this indicates that there is a division 
line in magnetic respect very close to the 
rotational equator. With this theory as a 
background we must suppose that waves with 
a certain polarity reach the northern hemisphere 
exactly down to the equator, and that waves 
of opposite polarity reach the southern hemi- 
sphere. Hence the magnetic equator must 
almost coincide with the rotational equator. 

A second argument is obtained from the 
sunspot progression curve. According to the 
theory this curve should be given purely by 
the geometry of the field (Cosm. El. p. 118). 
The fact that the progression curve is very well 
reproduced during all the sunspot cycles, which 
have been studied, speaks in favour of a regular 
magnetic field. 

Our third—and most important—argument 
is derived from the correlation between the 
spottedness pattern of the two hemispheres. An 
excess of sunspots at a certain latitude of one 
hemisphere is connected with a corresponding 
excess at the same latitude of the opposite 
hemisphere during next cycle. Judging from 
Fig. 5.16 in Cosmical Electrodynamics there 
is an agreement between the corresponding 
latitudes within a few degrees. This can only 
be interpreted in such a way that magnetic 
lines of force intersecting one hemisphere at a 
certain latitude, appear at the other hemisphere 
at the same latitude. 

The accuracy can also be judged from the 
fact that if the spottednes is plotted as a function 
of the latitude with latitude intervals of 2°, the 
correlation is poor, but if the curves are 
smoothed with the weight 1:3:5:3:1 a good 
correlation is obtained. This indicates that the 
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accuracy by which a certain latitude is con- 
nected with the corresponding latitude of the 
other hemisphere is a few degrees. 

A further confirmation of the regularity of 
the solar field is given by WauTTLE’s analysis. 
The fact that the time series of 16°—20° north 
is correlated with the series of 16°—20° south 
shows that lines of force must connect these 
belts in a regular way. 

The first ZEEMAN effect measurements of the 
general magnetic field indicated an angle of 
5° between the magnetic and rotational axes. 
It is now generally agreed that the accuracy 
was insufficient to establish this, and there is no 
later confirmation. ABETTI has made an inter- 
esting study of the shape of the magnetic field 
near the poles as found from the structure of 
the corona (ABETTI, 1952). The difference be- 
tween the magnetic pole and the heliographic 
pole which he finds refers to the instant 
structure of the field and not to the time 
average, which is what we discuss. 


6. Identification by position 


Our survey has shown that it is likely that 
magnetic lines of force of the general field 
intersect the two hemispheres at the same 
latitudes. Hence one of the necessary conditions 
for a sunspot identification by position is given. 
On the other hand we do not know the 
difference in longitude between the two “ends” 
of a line of force, so we must introduce this 
quantity as an unknown parameter. It should 
be kept in mind that the period of the solar 
rotation is not very well known. Further the 
rotation is non-uniform. Hence it is difficult 
to keep account of the longitudes. 

We must now ascertain how an attempt to 
identify sunspots should be made in order to 
have some chance of success. 

The regularity of the magnetic field is prob- 
ably best near the equator. This is indicated by 
the fact that the equator strictly separates spots 
of different polarities(WALDMEIER). The pattern 
analysis discussed above (ALFVÉN, 1948) which 
refers to latitudes in the range 10°—30° does 
not give a higher accuracy than a few degrees 
for the correspondence in latitude. 

The non-uniform rotation introduces a 
difficulty into the analysis. Suppose that a 
disturbance travels along a line of force which 
intersects the solar surface at latitude œ and 
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another disturbance travels along a line of 
force through y+/\ 9. Suppose further that 
at a certain time the two disturbances have the 
same longitude. A certain time t later a longi- 
tude difference AA will be produced by the 
non-uniform rotation. Adopting for the 
rotational velocity «© Newron’s formula 
O =W)—w’ sin? 
we find 


AA = o’t [sin? (p + Ag) — sin? y] 
wo't sin 2 p+ Ap 


Putting t = 4,000 days and w’ = 2.8 degrees 
per day we find 


AA = 195 sin 2 g:/\y 


where /\y as well as /\/ is expressed in degrees. 
It is unlikely that we can find or even define 
the latitude of a sunspot group within less than 
a couple of degrees. Putting Ag = 2° we find 
that AA > 90° unless  < 7°. This means that 
for latitudes higher than 7° the uncertainty in 
longitude becomes more than 90° after one 
solar cycle. Consequently the possibility of 
identifying spots by their longitude becomes 
small unless we confine ourselves to very low 
latitudes. 

There is also a third argument for choosing 
low latitudes, viz. that the number of sunspots 
is relatively small near the equator. This 
should make an eventual identification more 
unambiguous. 


7. Construction of spot diagrams for the 
equatorial belt 


Hence we confine the analysis to a small 
belt near the equator. All spot groups which 
in Greenwich “General Catalogue of Sunspots” 
or “Ledgers of groups of sunspots” are listed 
as having a mean latitude o<p< 4° or 
0o>9> — 4° are selected. For each cycle and 
hemisphere à diagram is plotted with the time 
“First seen” as abscissa and “Mean position of 
group, longitude” as ordinate. If the “Mean 
Area, corrected for foreshortening, whole 
spot” is less than so the spot group is repre- 
sented by a dot. Larger spots are marked by 
the symbols shown in Fig. 3. The longitude 1s 
iterated in order to show the pattern even 
near 360°. A recurrent spot appears in the 
diagram as two dots close together. 
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8. Comparison between the patterns 


According to the magneto-hydrodynamic 
theory of sunspots there is a possibility that the 
patterns of opposite hemispheres during the 
two consecutive cycles should be similar. In 
order to investigate whether this is the case we 
should compare e.g. the diagram 15 S (for 
the southern hemisphere during the ısth cycle) 
with 14 N and 16 N. Further we should 
compare 15 N with 14 S and 16 S. When 
comparing the diagrams we do not know a 
priori how the diagrams should be placed in 
longitude because the solar magnetic field may 
be twisted around the axis as discussed in § 4. 
Further it is difficult to keep account of the 
longitudes especially with regard to the non- 
uniform rotation. Hence the longitude displace- 
ment of two diagrams is arbitrary. 

The exact time difference between a spot on 
one hemisphere and the opposite one is also 
unknown a priori, but it should be of the same 
order as the sunspot period, say in general 
between 10 and 12 years. 

For the comparison we take the two 
diagrams (which have been plotted on trans- 
parent paper) and place one on the top of the 
other in such a way that one of the biggest 
spots of one diagram coincides with one of the 
biggest spots of the other diagram. If in this 
position most of the spots coincide the diagrams 
are judged to be similar. 


9. The sequence 13 S—14 N—15 S—16 N—17 S 


A comparison between the patterns is most 
unarbitrary when there are only a few spot- 
groups, and these form a distinct pattern. This 
makes it reasonable to start with the diagrams 
15 Sand 16 N. The biggest spot groups during 
15 S occur in 1919 and at the end of 1922. In 
16 N there is no big spot group later than 
1932. Let us see whether it is possible to iden- 
tify the spot group of 1932 as being caused by 
the same disturbance as the spot group of 1922. 
The difference in time between the two spot 
groups is about 91/, years, the difference in 
longitude about 50°. If we place the two 
diagrams 15 Sand 16N on top of each other we 
find that the medium size spot group in 1921 
falls close to the spot group of 1931. A series of 
spot groups in 1919 and 1920 could very well 
be of the same origin as a similar series in 1929 
and 1930. So far we have only considered large 
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Fig. 4. Comparison between the diagrams of fig. 3. The northern hemisphere is black, the southern hemisphere 
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red. The diagrams are placed in the positions indicated by the values of A and T given in Table 2. 
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or medium size spot groups. There is no 
indication of a general correspondence between 
the small spot groups, and these should not in 
general be used in our analysis. However, in 
1917 there is a small spot group which is 
remarkable because it is the only spot group 
in that year and it introduces the new cycle. 
This spot group corresponds in a notable way 
to the first spot group in the cycle 16 which 
occurs in 1927. In the same way the first 
medium size spot group in 15 S, which occurs 
in the beginning of 1918 corresponds to a 
small spot group in the beginning of 1928. 

The result of the comparison is shown in 
Fig. 4 where cycle 15 S is marked by black and 
16 N by red. The diagrams are placed in 
such a relative position that the best agreement 
between the patterns is obtained. The spot 
groups which are identified are marked by a 
letter. The identifications a, e, f, and g seem 
to be unambiguous. In all these cases one 
group in 15 S corresponds to one group in 
16 N. The groups 15 Se and 16 N g are recur- 
rent and as they appear during two consecutive 
solar rotations they are marked by two dots in 
the diagram. 

The identifications b, c, and d are perhaps 
more doubtful. It is obvious that there is a 
general spot activity in 1919 and 1920 at the 
longitude 100°—250° and that there is a 
similar activity during 1929 and 1930 at 
corresponding longitudes. The big spot group 
in 1919 was accompanied by some small spots. 
It is tempting to identify the two medium size 
spots in the middle of 1929 as being due to 
the same disturbance. This means that we 
identify what we could call a “cluster”? of spot 
groups with another “cluster” or single group. 
Such clusters are b, c, and d in 16 N, identified 
with the cluster b and the spots c and din 15 S. 
The group marked k in the middle of 1929 
might belong to the cluster b, but it is more 
likely that it should be considered as a group 
with no correspondence. Of course we could 
not expect all spots to be reproduced. 

The identifications are listed in Table 1. The 
difference in time AT between a spot in 15 S 
and its reproduction in 16 N is calculated. We 
shall call AT the reproduction time. Also the 
longitude difference and latitude difference 
are given. For clusters the weighted means are 
used, where the weight of the groups in a 
cluster is put equal to the spotted surface. In 
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the figures the displacement in longitude and 
in time is made equal to the average of all the 
identified cluster (inclusive single groups). 

The similarity in pattern between 15 S and 
16 N seems to be too striking to be due to 
chance. We know from Gatventius’ and 
Wozp’s and WHITTLE’s papers that there is a 
statistical correlation between the two hemi- 
spheres. Combined with these facts our results 
give more support to the general lines of the 
magneto-hydrodynamic theory. Of course it 
would be highly desirable to calculate the 
statistical significance of the similarity between 
the patterns. 

The positive result of this comparison makes 
it desirable to compare all other pairs of 
diagrams. 

We next proceed to 16 N and 17 S. There 
are two big clusters in 16 N, viz. fand g. In 17S 
there are three big clusters, one in 1940, one 
in 1941, and one in 1942. The relative time of 
appearance of these clusters agrees rather well 
with the sequence d, f, and g in 16 N, but if 
the diagram is placed in such a way that the 
clusters of 1940 and 1941 coincide with d and f 
there is a disagreement in longitude of 125° 
between g and the cluster of 1942. Further 
there are no other clusters coinciding. Hence 
we must reject this identification and suppose 
that the cluster of 1942 (called “v”) has no 
correspondence, in the same way as the group 
k in 16 N has no correspondence in 15 S. 

In next attempt we assume that the clusters 
of 1940 and 1941 correspond to f and g. If we 
place the diagrams in this way (Fig. 4) we find 
in 17 S spots corresponding to k and din 16 N. 
Further there are two spots which correspond 
to the cluster c. As there seems to be no way 
to obtain a better correspondence between the 
patterns of 16 N and 17 S we accept these 
identifications. In 16 N ‘there are five clusters 
(k, c, d, f, and g) which are reproduced. There 
are three clusters a, e, and b which are not 
reproduced, but two of them consist of only 
one small spot each and were discussed earlier 
because of their peculiar situation in the 
beginning of the cycle. In 17 S there is one 
large spot (v) which is new. (The last spot in 
16 N may be a correspondence.) Further there 
are now some medium size spots (y and five 
spots between f and-g). 

The agreement between 16 N and 17 S is 
not so striking as between 15 Sand 16 N. The 
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identification of 16 N d, f, and g with 17 S 
f, g, and v could not be altogether ruled out 
but is definitely less convincing than the 
identification we have accepted. 

As the spot material for the next cycle is 
not yet available we proceed backwards in 
time and compare 14 N and ı5 S. The big 
cluster 15 S b could correspond either to one 
of the big spots in 1906 and 1907 or to one 
of the big spots in 1909 and 1910. According 
to the first alternative we could identify b, c, 
and d with the three big spots in 1906 and 
1907. No other identifications are obtained in 
this way, and the reproduction time becomes 
as large as 13 years. The second alternative is 
shown in Fig. 4. The clusters 1909 and 1910 
are identified with b, c, and d. The small spot 
a is considered as a reminiscence of one of the 
big cluster of 1907. No identification of e, f, 
and g is obtained. Hence this alternative is only 
slightly better than the first one. It leaves us, 
however, with a more reasonable value (10 
years) of the reproduction time. 

As is evident from the comparison 15 S—16 
N the clusters b, c and d form one big region of 
sunspots: and the division of it into clusters is 
somewhat arbitrary. The comparison 14 N—15 
S shows that the corresponding region of 
sunspots is found in 1909—1910, but the 
subdivision into clusters is somewhat arbitrary. 
An identification of 15 S b, c, and d with 
14 Nm, b, and c is not ruled out although less 
probable than the identification which we have 
accepted. 

Comparing 13 S and 14 N we find in both 
cases two maxima of spot activity (1894—1895 
and 1897 in 13 S and 1906—1907 and 1909— 
1910 in 14 N) separated by about one year 
(1896 and 1908), of low activity. During the 
two maxima of 13 S the spots have a rather 
even distribution in longitude so that the 
pattern is not very distinct. The identification 
shown in Fig. 4 is no doubt the best one but 
it is not absolutely convincing because it leaves 

‚us with quite a few unreproduced clusters. 

The treatment of the sunspot material has 
not yet been extended more backwards in 
time. We have thus treated the whole sequence 
13 S—14 N—15 S—16 N—17S. There is no 
cluster which can be followed the whole way 
through the sequence. Our identifications 
depend to a considerable extent on the spot 
region which in 13 S is represented by k, |, 
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and b with an extension to d. This is reproduced 
in 14 N by k, |, b, ¢, and d. In 15 Sandı6 N 
the same region consists of b, c, and d, whereas 
it fades out in 17 S with c and d as the only 
representatives. 

A new region of spots appears in 15 S with 
f and g, which are reproduced in 16 N. By 
being the only spots in the later parts of each 
cycle they are very valuable as identity marks. 
There are no signs at all of fand g in 14 N. 
Going to 13 S, however, we find in 1901 one 
single small spot which ends the cycle. This 
has the same position as a predecessor of g 
should have, so a genetical connection is 
possible. (The fact that no corresponding spot 
is observed in 14 N does not contradict this 
because there may have been one on the sun’s 
backside.) There is also a possible correspond- 
ence to f in 13 S. The spots f and g develop 
further in 17 S where a number of new spots 
are created between them. Two new big 
groups v and y appear in the late part of 17 S. 

In 13 S there is an extended spot region 
from o and h to n andr. This is also prominent 
in 14 N but it. dies away in 15 S. The only 
remnant is the small spot a which even survives 
to 16 N. 

Consequently in the sequence 13 S—14 N— 
15 S—16 N—17 S there is a tendency for spot 
regions to form in the late part of a cycle, 
proceed to the middle of the cycle where they 
develop their greatest intensity and finally die 
at the beginning of the cycle. 


10. The sequence 13 N—14 S—15 N—16 S—17N 


We shall now make the same analysis for the 
sequence 13 N—14 S—15 N—16 S—17N. 

In 13 N there is a number of spots in 1893 
and 1894. After a lull in 1895 and 1896 there 
comes a long train of clusters «, 6, &, ¢, all at 
about the same longitude. Our earlier results 
indicate that the first spot of the cycle are in 
general not very well reproduced, so we 
should primarily try to find a correspondence 
to a, 6, &, and ¢. A similar long train is found 
in 14 S. This leads to the identification shown 
by fig. 4. A cluster 14 S x seems to correspond 
to the last small spot in 13 N. Further, at 
about the same time as « appears there is a 
number of spots also at other longitudes. An 
identification of the spots 7 and f seems 
possible. The small spots in 1905, 1906 and 
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1907 seem not to correspond to the activity in 
13 N during 1893—94 but the first small 
spots w of each cycle agrees very well. Upon 
the whole the similarity 13 N—14 S is satis- 
factory. 

Proceeding to 15 N we recognize easily the 
long train «, 6, &, €. The cluster y which was 
born in 14 S is very well reproduced. An 
identification of 14 S ß, which is also a new- 
born spot, with a spot in 15 N is indicated. The 
clusters e and € have grown but there is no 
correspondence of x. New-born clusters are ju 
and v. Further there is a strong activity in 
the beginning of the cycle, especially a cluster 
1s N# containing two very big and four 
medium groups. There is no correspondence 
to this in 14 S but the cluster 13 N @& coin- 
cides with it. 

The comparison 15 N—16 S is more difficult. 
The long train has been dissolved. The biggest 
cluster in 16 S occurs at Christmas 1929. This 
could obviously not be identified with the 
biggest cluster ¢ in 15 N because this would 
give a reproduction time of only 8.5 years. 
An identification with the second biggest 
cluster y seems to be acceptable (fig. 4). The 
clusters &, B, y, 6, u and possibly » are repro- 
duced in a reasonable way. Of the big cluster 
e only one medium group remains, and & is 
not reproduced. The big cluster # may be 
reproduced by four small spots which start 
the cycle. A new cluster is &. 

In 17 N the prevailing feature is again a 
long train of spots, even longer than in 13 N 
and 14 S. The best identification seems to be 
to associate this with 9—&—n—e (fig. 4) This 
means that «, ß, 6, and » have died and two 
new clusters g and y have been born. 

In the “Greek sequence” 13 N—14 S—15 N 
16 S—17 N, which we have now examined, 
some clusters have a remarkably long duration. 
In fact, #, 6, and & have persisted during the 
whole time (with an interruption for 9 in 
14 S). If 17 N is placed on top of 13 N these 
three clusters are easily brought to coincide. 
If all the five diagrams 13 N—14 S—15 N— 
16 S—17 N are placed on top of each other 
these three clusters are in about the same place 
in all of them. 


11. Correlation between the two sequences 


We have found a good correlation within 
the “Latin sequence” 13 S—14-N—15 S—16 
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N—17 S and also within the “Greek sequence” 
13 N—14 S—15 N—16 S—17N. It is 6 
interest to see whether these series are inter- 
connected. A comparison of the type (n, N) — 
— (n+1, N) or (n, S)— (n+1, S) or of the 
type (nN) — (nS) shows in may cases some 
degree of similarity. 


12. Check of identifications 


First check: comparison of non-related 
patterns 


It is of decisive importance to investigate 
whether the pattern similarities which we have 
found are significant or not. It is quite possible 
that two patterns are judged to be similar even 
if there is no real correlation between them. 
The first way to check this is to compare 
diagrams which we do not expect theoretically 
to be similar. A difficulty is that if there exists 
a correlation between the two sequences, all 
diagrams without exception ought to show 
some correlation. In order to compare diagrams 
with as small theoretical correlation as possible 
we should choose them from different se- 
quences, and with as long a time interval as 
possible. Hence we select the four pairs 
13 N=16-N, 14 N—17 N, 13 S—16 S, and 
14 S—17 S as material for our check compari- 
son. Proceeding in the same way as earlier we 
place the biggest spot of one diagram on the 
top of one of the biggest spots of the other 
diagram. 

The comparison shows that 13 S and 16 S 
gives a not too bad similarity, 13 S t, |, b, 
and d forming a pattern similar to 16 S y, 6, 
x, &. In the three other cases no similarity at 
all could be traced. 

The negative result in three cases indicates 
that it is probable that our earlier identifications 
are significant but the fact that in one case a 
similarity could be found which mainly must 
be accidental, shows that our method is not 
quite without danger. Hence other checks are 


highly desirable. 


13. Second check: Calculation of the “twist” 
of the magnetic field 


Our second check is of physical character 
and concerned with the calculation of the 
twist of the magnetic field (see § 4). 

Suppose that the magnetic field is twisted in 
such a way that a magnetic line of force which 
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intersects the southern hemisphere at the 
longitude As arrives to the northern hemisphere 
at the longitude 


An = Às ONE (1) 


where © is a constant which we shall call the 
twist. Suppose further that a disturbance travels 
along a certain line of force and reaches the 
southern hemisphere during the cycle n at the 


longitude As. The same disturbance reaches 
the northern hemisphere during the cycle n +1 


after a time interval TS’; (the reproduc- 
tion time). If the latitude at which the spots 
occur rotates with the angular velocity wg in 
comparison to the latitude to which the 
heliographic longitude determination refers, 


the longitude at the northern hemisphere is 
n+1 an N, n+1 
AN =As+@O+ aaTs.'n ae) 
Consider another disturbance which reaches 


the northern hemisphere at longitude AN 
during cycle n and the southern hemisphere 


at longitude ae during next cycle. 
In a similar way we obtain 
as in —O HOsIN +--+ (3) 


Sead. ; ; 
where TN, is the reproduction time of the 
disturbance. We put 


Std n+1 n 
“IN, 1 — Ss —/N Pee 06 asic) (4) 
and 
N, n+1 „n+l n 
As, n = AN. — Às pases. SiG (s) 


The quantities A are the longitude differences 
which are found empirically at the comparison 
between the diagrams. Subtracting (3) from 
(2) we obtain 
I N. n+ 
O os (As. a 


3 ANY) — 5 oT (6) 
with N T= Te = Th ee 

We do not know w, with certainty. We do 
not even know whether it is constant. In fact, 
some of the results of the present analysis 
indicate that this is not the case. The differ- 
ential rotation is usually computed from the 
systematic longitude drifts of sunspots which 
for our latitude belt is 5° per solar rotation, 


1 I ; , 
corresponding to 5 Ma= 0.09 degrees /day. 


From the point of view of the magneto- 
hydrodynamic theory of sunspots it is con- 
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ceivable that systematic drifts may be caused 
also by other effects than a differential rotation. 
However, all these problems will be reserved 
for a later analysis. In view of the uncertainty 
we shall here calculate © in two ways, viz. by 
neglecting the second term in (6) and by 
including this term. The latter value will be 
referred to as Or. 

Table 2 shows the values of © calculated 
with and without the second term. The 
constancy of @ is remarkable. In fact, without 
correction it varies only within the relatively 
narrow limits 46° and 96°. The excessive 
agreement between the corrected values is of 
course only a numerical coincidence. As mean 
we obtain 

O = 66° + m 180° 


where m is an integer. 

The constancy of © gives a strong support 
for our identifications. In fact, if we accept all 
identifications except one, we could calculate 
the relative longitude position A for this one 
from (6). The position would agree with what 
we have found within an error of less than 
about 45°. This shows that in no case a serious 
error can have been made when fixing the 
relative longitude position of the diagrams. 


14. Third check: Correlation between repro- 
duction time and the length of the solar 
cycle 

Having thus shown that the identifications 
must be correct in longitude, it remains to 
investigate whether there is any possibility to 
check the values of the reproduction time. It 

could be suspected that there should be a 

positive correlation between this and the 

length of the solar cycle. In tact, if sunspots 
are reproduced after a. certain time IT we 
should expect that this quantity regulates the 
length of the solar cycle. It is obviously not 
the only factor of importance, because as we 
have seen not all spots are reproduced. There 
is a general tendency that spots in the beginning 
of a cycle are not very well reproduced and 
that new spots are born at the end of the cycle. 
Table 2 shows a comparison between the 


average of Tx and Tea and the length 
of the solar cycle (= time difference between 
two sunspot maxima). The correlation is ob- 
viously positive. This is satisfactory, but on 


the other hand the correlation is not good 


H. ALFVEN 


444 

Table 2. Results of the comparison. 
. Diff. 

S,n+r N,n-+ı S,n-+ı N,n+1ı 

Cycle IA (©) OAI Te AT T'xean | Sunspot 
N,n S,n N,n S,n max. 
13—14 — 47° | + 84 66° 64° 11Y 1484] 11Y 1694 214 11Y 1594 12.9 
14—15 — 11° | + 111° 62° 62° 10Y 724 | 10Y 744 24 10ŸY 734 10.6 
15—16 — 95° | — 3° 46° 033 10% 1334] 9¥ 3054 |— 1934 10Y 364 10.8 
16—17 + 5° + 197° 96° 74° 9Y 1584 | 10Y 424 2494 9Y 2834 9.0 
Mean | 68° | 66° 


enough to give an accurate check of our values 
of the reproduction time. 


15. Identification by structure 


In § 3 it was pointed out that it may be 
possible to identify two spots as produced by 
the same original disturbance by discovering 
some peculiar property which is common to 
both. If for example a sunspot whirl ring 
makes a certain angle with the magnetic field, 
this may cause correlated proper motions of 
one spot of the northern hemisphere and the 
corresponding spot of the southern hemisphere. 

This could preferably be checked by com- 
paring the proper motions of spots which we 
have identified by their position. Simple spots 
like f and g are especially suited for this. A 
preliminary analysis indicates that there may 
be a correlation between the proper motions, 
but the problem must be studied much more 
carefully before any definite conclusions can 


be reached. 


16. Spots at higher latitudes 


The whole analysis has so far been confined 
to the latitude belt 4° S—4° N. An extension 
to higher latitudes is very desirable. It seems 
relatively easy to proceed up to 6° or may be 
even 8° by the same method, but when we 
approach the maximum belt at 10°—15° the 
number of spots gets so large, that it is im- 
possible to distinguish distinct patterns. It will 
be necessary to develop other methods of 
identification before a real attack towards 
higher latitudes could be launched. If the 
identification by structure succeeds, the prob- 
lem will be easier. 

It is of interest to note that the survival of 
spot groups during several cycles which was 
found in $$ 9 and 10 is in good agreement 
with earlier results concerning the latitude 
pattern of spottedness (cf. Cosmical Electro- 


dynamics p. 138). The latitude distribution for 
cycle 14 is very similar to that of cycle 15 and 
16, and there is also a definite resemblance for 
cycle 12. These results refer to higher latitudes 
(10°— 30°) and indicate that the persistency of 
spots may be a general phenomenon. 


17. Theoretical conclusions 


The reproduction of sunspots which has 
been found here is well explicable according 
to the magneto-hydrodynamic theory of 
sunspots. A wave which reaches e.g. the 
southern activity region becomes in part 
amplified, and in part reflected. The reflected 
wave goes to the northern activity region 
where the same phenomenon occurs and a re- 
flected wave is sent out to the southern region. 
This can be repeated several times. 

In the theory it was assumed that when a 
whirl ring passes an activity region it usually 
gives rise to a series of reflected whirl rings. 
None, or one, or many of these may survive 
the destruction by damping and other causes, 
and finally give rise to new spots. The activity 
regions were thought to consist of the upper 
parts of the convective core and to be rather 
thick, perhaps a few tenths of the size ofthe core. 

This idea of the activity regions is well 
reconcilable with the present results except that 
the activity regions cannot be very thick. A 
sunspot group may be reproduced next cycle 
by a cluster of groups, a single group, or not 
at all. If it is reproduced by a cluster of groups 
the spread in time of the spots of the cluster is 
a measure of the thickness of the activity 
region, where they are produced. Suppose that 
a wave pulse with the velocity V passes a 
region with thickness D in which it may be 
partially reflected anywhere. The difference in 
time f between pulses reflected at the beginning 
and at the end of the region is 


t=2D/V 
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In our case f represents the maximum spread 
in time within a cluster. Judging from the 
diagrams it is reasonable to put  — 108 sec. (+ 4 
months). The velocity V is likely to be of the 
order 60 cm/sec (Cosmical Electrodynamics 
p. 100). This gives 3 - 10° cm which is only 
0.004 of the solar radius. 

Consequently the activity region is so thin 
that it more resembles a surface of discon- 
tinuity, which should have the property of 
reflecting and amplifying waves passing it. 
The existence of such a surface is hardl 
reconcilable with current solar models. Gamow 
(1945) and Reız (1948, 1949) have discussed a 
model in which the energy is released in a 
spherical shell. It is conceivable that magneto- 
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hydrodynamic waves passing such an energy 
producing shell should be reflected and am- 
plified in the way our results require. The 
possibility to construct a solar model of this 
type should be seriously discussed. When we 
learn to correlate more closely the appearance 
and structure of sunspots with the phenomena 
in the solar core by which the spots originate, 
sunspots data will give us valuable information 
concerning the interior of the sun. 

I wish to thank Dr. N. Herrorson for numer- 
ous discussions. Mrs. B.-M. Hoımum has 
kindly extracted the sunspots data from the 
catalogues and Mr. Wirt has prepared the 
diagrams. 
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A Critical Examination of Theories of Charge Generation in 


Thunderstorms 


By B. J. MASON, Imperial College, London 


(Manuscript received May 2, 1953) 


Abstract 


A concise review is presented of the more important features of the meteorological and 
electrical behaviour of thunderstorms from which the main conditions to be fulfilled by 
any satisfactory theory of charge generation, are enumerated. The various theories are then 
examined, as far as possible from a quantitative point of view, in the light of these requirements. 
None is found to be entirely satisfactory. It is concluded that the mechanism which most accords 
with the observed facts is the charge generation and separation which, according to the results 
of laboratory experiments, will be associated with the growth of graupel pellets in the super- 


cooled region of the thundercloud. 


I. Introduction 


During the last forty years, at least eight 
different mechanisms for the generation of 
electricity in thunderclouds have achieved 
prominence in the literature. Too often theories 
have been formulated without sufficient regard 
to all the available experimental data on the 
behaviour of thunderstorms, and rarely have 
they been examined for quantitative agreement 
with these data. Any satisfactory theory must 
show quantitatively how charge can be created 
and separated at the rate required by the 
measurements of the electric field-change and 
the frequency of lightning discharges and must 
conform with known facts about the meteoro- 
logical and electrical structure of the cloud, 
its life-history and the nature and time-scale 
of the precipitation processes. 


2. The Electrical and Meteorological Behav- 
iour of a Thunderstorm 


(a) The Life Cycle of a Thunderstorm 


The investigations made during the Thun- 
derstorm Project (see ByErs and BRAHAM, 


1950), indicate that a thunderstorm consists of 
one or more identifiable units or cells which 
are localized regions of pronounced convec- 
tive, microphysical and electrical activity. The 
life of an individual cell can be divided into 
three stages, depending on the direction and 
magnitude of the vertical air currents. The 
cumulus, or growing stage is characterized by 
the existence of an updraught throughout the 
cell; the mature stage by the presence of both 
updraughts and downdraughts, at least in the 
lower part of the cell, the onset of lightning 
activity and the appearance of precipitation 
below cloud base; the dissipating stage by 
weak downdraughts prevailing throughout 
the cell with marked reduction in the intensity 
of precipitation and electrical activity. 

The duration of the cumulus stage reckoned 
from the time of initial detection of the radar 
echo, which appears soon after the visible 
cloud top extends above the 0°C level, is 
between ro and 15 minutes. The mature stage, 
during which most of the precipitation is 
released, lasts for :1$ to 30 minutes. The 
duration of the dissipating stage, measured to 
the time when vertical motions inside the 
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cloud become insignificant, is of the order 
30 minutes. Thus the life of one cell occupies 
about one hour, but the duration of the radar 
echo and its associated electrical activity may 
be confined to a period of less than 30 minutes. 
The larger, more persistent thunderstorms 
which may last for several hours probably 
result from the activation of a group of cells. 
Successive centres develop in turn and more 
than one may be active simultaneously. 


(b) The Meteorological Structure of a Thunder- 


storm 


Most of our information on the structure 
of the thunderstorm again comes from the 
data of the Thunderstorm Project and from 
the observations made at the Observatory on 
the Zugspitze (altitude 10,000 feet) and 
reported by KUETTNER (1950). 

In the early stages, the incipient thunder- 
cloud grows rapidly in the vertical and is 
characterized throughout its volume by up- 
draughts whose magnitude was 7 m/sec. on 
average in the storms investigated by the 
Thunderstorm Project, but values of 15 m/sec. 
are not uncommon and may be appreciably 
exceeded. The horizontal dimensions of the 
cell, as first detected by radar, may be only 
0.5 km with the top extending only about ı 
km above the 0°C level. Thereafter, the 
visible cloud top may grow rapidly towards 
the — 30°C isotherm, or even higher. Byers 
and Braham report observations which in- 
dicated that within the updraught region of 
a single cell there was more than one locality 
of strong vertical ascent; that these regions 
were separated in space or time, or both; and 
that they resulted in the development of in- 
dividual cloud columns appearing as turrets at 
the summit of the visible cloud. Kuettner also 
quotes evidence for this sub-division of a cell, 
the smaller units having horizontal dimensions 
of 100 m to ı km, ie. about one-tenth those 
of the cell. They often occur at different posi- 
‚tions relative to the centre of the cloud and 
successive turrets often penetrate to greater 
heights than their predecessors. These sub- 
cells may perhaps be identified with the 
“bubbles” which are regarded by Scorer 
and Lupram (1953) as the fundamental units 
in the structure of cumulus; a succession of 
rising bubbles would confer on the updraught 
something of the pulsating nature mentioned 
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by Byers and BRAHAM and by Workman and 
Horzer (1942). 

The last-mentioned authors state that, in 
the early stages of cloud development, the 
water particles are small but their size con- 
tinually increases with time. Later, as the 
cloud top reaches colder regions, larger parti- 
cles of snow, graupel and rain appear but 
continue to be carried upwards by the up- 
draught which may now be typically 1.5 km 
in diameter. The greatest concentration of 
hydrometeors occurs at and above the 0° C 
level. The size and concentration of the pre- 
cipitation elements continue to increase until 
they can no longer be supported by the existing 
updraught and they then begin to fall relative 
to the earth. At about the same time, strong 
downdraughts appear in the lower part of 
the cell, particularly in the region where 
precipitation is falling. This marks the onset 
of the mature stage of the cell during which 
it attains its maximum height and may reach 
the tropopause; the horizontal dimensions of 
the radar echo may now lie between 5 and 
14 km at the base, tapering off towards the 
top. There may still be strong local updraughts 
particularly in the upper parts of the cloud 
where the average velocities are about 10 
m/sec., but on occasions may exceed 30 m/sec. 
Precipitation elements consist mainly of rain- 
drops below, and snow or graupel particles 
above, but hail may occur in some cells, par- 
ticularly those containing exceptionally strong 
updraughts. The mature stage is gradually 
transformed into the final dissipating stage in 
which the whole cell contains air moving 
downwards or almost at rest. Precipitation 
continues to fall for some time but gradually 
diminishes in intensity and ceases. 

Kuettner presents some very valuable in- 
formation on the nature of hydrometeors 
found in thunderstorms; his observations being 
made on a mountain top immersed in the 
cloud are more reliable and clear-cut than 
those made from aircraft. He reports that 
solid precipitation elements were predominant 
in the greater part of the thundercloud, being 
found on 93 % of the occasions, while rain- 
drops were present in only 21 % of the clouds. 
Kuettner states that graupel (including snow 
pellets and soft hail) was the most frequent 
form of hydrometeor, being present on 75 % 
of occasions; large hail was relatively rare and 
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was, by no means, a necessary event in thun- 
derclouds. Snow crystals prevailed- during the 
later periods of a storm when the precipita- 
tion became more steady and moderate in 
character and the lightning activity diminished. 

Both Kuettner and WORKMAN and Rey- 
NOLDS (1950 a) state that a necessary condition 
for the development of thunderstorms on the 
Zugspitze and in New Mexico respectively, 
is that the cloud base should be warmer than 
o° C. Furthermore, Workman and Reynolds 
state that the first lightning flash occurs only 
when the top of radar echo approaches the 
- 30°C level; Byers and Braham say that the 
top of the echo must be colder than — 20° C. 
Workman and HOIZER (1942) report no 
lightning from clouds of depth less than 5 km. 


(c) The Electrical Structure of a Thunderstorm 


From his measurements on the change of 
the electric field at the ground during lightning 
activity and the variation in the sign of the 
field-change with distance from the storm, 
Witson (1920) deduced that the thundercloud 
has a bipolar structure with a positive charge 
situated above a negative charge. Wilson’s 
conclusions have been amply confirmed by 
similar measurements by WorKMAN, HOLZER 
and PELSOR (1942) in New Mexico, who 
found some evidence which pointed to the 
upper charge-centre being displaced in the 
direction of motion relative to the lower 
negative centre (this is to be expected in the 
presence of a wind-shear). The positive and 
negative charge centres were, on average, 
2.5 km and 4 km respectively, above cloud 
base, but probably only in the early stages; 
Workman and REYNOLDS (1950 a) report that 
as the storm progresses, the positive centre 
moves upwards while the negative centre 
appears to remain fixed in the zone of heavy 
precipitation. 

MALAN and SCHONLAND (1951 a, b), from 
a study of the electric-field changes in the 
intervals between successive strokes of light- 
ning flashes to earth, conclude that the negative 
charge which is involved in such a flash is 
contained in a nearly vertical column, which, 
on occasion, may be 6 km long, the tempera- 
ture at its base being not much below 0° C 
and its top frequently reaching the — 40° C 
level. Their techniques allow determination of 
the heights of the charges involved in separate 
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strokes; the results obtained from five different 
methods agree in indicating that successive 
strokes tap progressively higher and higher 
regions of the negative column, the average 
increase in height per stroke being about 
0.7 km. The first stroke appears to originate 
from near the — 5° C level, while later strokes 
may extend upwards by more than another 
s km. Malan and Schonland also report that, 
in their South African thunderstorms, 65 % of 
the first strokes of flashes to earth have type B 
leaders indicating that there is a positive 
space-charge below the base of the main 
negative column. The discharge of the negative 
column to earth may therefore, in many cases, 
be triggered off by the strong field between 
the positive pocket and the bottom of the 
negative column. 

WormeLt (1953) has made similar investiga- 
tions but finds that on 75 % of occasions there 
is no appreciable change of field during the j 
portion of the stroke, i.e. although there is a 
transport of charge, the electric moment is not 
changing. This indicates that in 75 % of the 
multi-stroke flashes to earth around Cam- 
bridge, England, the negatively charged 
regions in the cloud which are discharged by 
successive strokes, are all sensibly at the same 
height. In the remaining 25 % of cases the 
form of the field change indicates that the 
heights of these regions increases for successive 
strokes as found by Malan and Schonland. 

The most direct information on the electrical 
structure of a cumulonimbus comes from 
measurements of the electric field made from 
balloons and aircraft. The pioneer work was 
carried out at Kew by Simpson and SCRASE 
(1937) and Sımpson and ROBINSON (1941), 
who used the altielectrograph on a balloon to 
determine the sign of the vertical component 
of the electric field as a function of height and 
also to obtain an estimate of its magnitude. The 
records established beyond reasonable doubt 
that the upper part of a thundercloud is 
positively charged and most of the lower part 
negatively charged. In 7 of 13 different storms 
there was also evidence of one or more localised 
regions in the base of the cloud containing 
positive charge. The upper positive charge 
was usually associated with temperatures 
below - 20°C, and always below - 10°C. 
The temperature at the negative charge centre 
was below 0°C in 13 out of ıs cases while 
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the lower positive charge, in 5 cases out of 7, 
was at temperatures above 0° C. Simpson and 
Robinson stated that their results could be 
represented schematically by a thunderstorm 
model with an upper charge of + 24 coulombs 
distributed in a spherical volume of 2 km 
radius centred at a height of 6 km (- 30° C), a 
negative charge of - 20C in a sphere of 
radius 1 km centered at a height of 3 km 
(—8° C) and a charge of + 4 C in a sphere of 
radius 0.5 km centred at 1.5 km (+ 1.5°C). 
Thus, the upper positive charge and the bulk 
of the negative charge were in regions where 
the temperature was below 0°C. Similar 
results regarding the charge distribution have 
been reported by CHAPMAN (1950) who used 
a radiosonde device, the observations being 
telemetered to the ground. The negative 
charge centres for 7 storms were reported 
to lie with 0.5° C of the o° C isotherm. 

The results of these balloon investigations, 
while of great value in establishing the general 
features of the charge distribution, cannot be 
regarded as very reliable for giving the tem- 
peratures of the charge centres and the actual 
magnitudes of the electric field, because it 
seems probable that on the majority of oc- 
casions, the balloon missed the electrical centre 
of the storm; in fact 15 % of Simpson’s useful 
ascents came to a violent end in a rapidly 
increasing field. The fields reported by Simpson 
et al, barely exceeded 100 V/cm while the 
largest field recorded by Chapman was 
210 V/cm; these are to be compared with 
values of 1,300 V/cm, and on one occasion, 
3,400 V/cm measured by Gunn (1948) from 
an aircraft. The last named author encountered 
the most intense fields near the 0° C level. 

KUETTNER (1950), from his observations on 
the Zugspitze, finds that the lower positive 
charge has a horizontal extent generally less 
than 1 km and is accurately centred at the 
0° C level. The more extensive negative charge 
is centred, on average, near the — 8° C level, 
while the main upper positive charge is still 
higher and even more diffuse. The negative 
charge centre is apparently the starting point 
of the lightning flashes. 

It is of some importance to reconcile the 
deductions of Malan and Schonland and of 
Wormell concerning the distribution of the 
negative charge in a thundercloud. In most of 
the South African storms it appears that the 
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negative charge is distributed in a nearly 
vertical column which is tapped at pro- 
gressively higher and higher levels by out- 
wardly-directed positive streamers during 
successive strokes of a flash to earth. This 
charge configuration might be expected in a 
cell containing only one region of strong 
updraught and one column of precipitation 
and electric charge. However, if there are more 
than one such regions separated in space and 
time throughout the cell, it is possible that a 
particular centre may become discharged by 
positive streamers directed outwards from a 
neighbouring centre which has already passed 
its charge to earth. Thus, the total charge 
dissipated in a lightning flash may be drawn 
from several charge centres of smaller dimen- 
sions than the thunderstorm cell and the dis- 
charge of a new centre by the next stroke 
might involve very little change in electric 
moment or electric field as measured at the 
ground. Such a mechanism may be responsible 
for the 75 % of Wormell’s records. 


(d) Correlation between Electrical Activity and 
Meteorological Processes 


Kuettner states that a high precipitation is 
an indispensable requirement for electrical 
activity in a cumulonimbus, the central light- 
ning area usually being coincident with the 
area of highest precipitation intensity. The 
onset of lightning activity coincides with the 
appearance of heavy solid precipitation within 
the cloud. 

Workman and REYNOLDS (1950) assert 
that separation of electric charge within the 
cloud is closely associated with the precipita- 
tion mechanisms and does not take place until 
precipitation processes are well under way; 
it takes place at temperatures between 0° C 
and — 15°C, most probably between — 5° C 
an — 10°C. The first internal flash occurs 
about 12 minutes after the initial radar echo, 
which by this time has approached the — 30° C 
level and started to descend; at the same time 
precipitation appears at the cloud base. Cloud- 
to-ground discharges follow within a few 
minutes. These authors also report that the 
lightning activity has a periodicity of 25 to 
40 minutes corresponding to the active life of 
a cell. 

Workman and HOLZER (1942) report that 
clouds of greatest vertical depth display the 
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greatest electrical activity. The vertical motion 
often appeared to be pulse-like in character, 
each new surge (or turret) being associated 
with increasing electric fields at the ground and 
a new series of lightning flashes. 
Measurements of field strength from aircraft 
made in conjunction with ground radar 
observations during the Thunderstorm Project 
showed that, during growth of the cell, the 
positive field (positive charge above, negative 
below), became strong in the precipitation 
area. In the dissipating stage the positive field 
was replaced by a negative one. A positive 
charge was often found in the rain core of the 
cloud near the cloud base. The maximum 
frequency of lightning preceded the onset of 
maximum intensity of rainfall at the ground. 


(e) The Electric Moment Destroyed in a Lightning 
Flash—Rate of Charge Dissipation 


From measurements of the change in the 
vertical component of the electric field at the 
ground during a distant lightning flash, one 
can calculate the electric moment 6M = 2HöQ 
destroyed, where ÔQ is the charge neutralised 
and H is the vertical separation of the main 
positive and negative charge centres for an 
internal flash, or the height or the involved 
charge above ground for a discharge to earth. 
WORMELL (1939, 1953) gives 110 C-km and 
130 C-km as the most frequent values for 
the moments destroyed during internal and 
cloud-to-ground flashes respectively, in storms 
investigated at Cambridge. Wırson (1920) 
gave an earlier mean value of 100 C-km, while 
SCHONLAND (1928) quotes 90 C-km as the 
appropriate mean value for South African 
thunderstorms. Taking a mean value of ro 
C-km for all types of flashes, values of H 
varying from 2 to 3.5 km for cloud-to-ground 
discharges and up to 5 km for internal flashes, 
the charges neutralised are calculated to be of 
the order 10 to 30 coulombs. WORKMAN and 
HOIZER (1942) report that the magnitude of 
the charge involved in New Mexico storms 
varies between 10C and 190C, the most 
frequent values lying between 20 C and so C. 
Horzer (1950) states that a multiple-stroke 
flash which may neutralise about 20 C, very 
often originates from several charge centres, 
each of which may contribute from about 
3 to 8 coulombs. He gives the average interval 
between strokes of a cloud-to-ground flash 
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as 0.05 sec.; Wormell gives 0.07 sec. and 
Schonland 0.03 sec. Taking ALLIBONE and 
Meex’s (1938) mean value for the positive 
streamer velocity as 2.4 x 10° cm/sec., and the 
average duration of a stroke as 0.05 sec., we 
arrive at 1.2 km for the average separation of 
the charge centres; this is a very reasonable 
value for the separation of the sub-cells or 
bubbles in the cloud. 

Lightning flashes from a given cell occur on 
average at intervals of about 20 sec., so that 
the current dissipated by lightning is of the 
order I amp. 


(f) Recovery of the Electric Field after a Flash; 
Rate of Charge Generation 


The recovery of the field after the passage 
of a lightning flash which corresponds to the 
regeneration of the electric moment in the 
cloud, is approximately exponential in form 
and according to WILSON (1929) and WORMELL 
(1939) has, on average, a time constant of 
about 7 sec. The initial rate of recovery re- 
presents the rate at which the electric moment 
grows when the internal field is small; the 
slowing down which follows is due to the 
increased electrical forces on the charged 
elements tending to oppose their separations 
as the internal field grows, and by the in- 
creased rate of dissipation by charges on the 
precipitation elements, by point discharge 
currents and ionic leakage to the upper atmos- 
phere. If we take the charge neutralised in an 
average flash to be 20 C and the initial rate of 
regeneration as 1/7 sec.~1, the vertical current 
in the cloud immediately after the discharge 
turns out to be 3 amps; however, the various 
leakage factors mentioned above, subsequently 
restrict the average net charging current to 
about 1 amp. 

The rebuilding of the electric moment, 
however, implies not only a generation of 
charge but a vertical separation of charges of 
opposite sign. If, as is assumed in most theories, 
charge of one sign is attached to the precipita- 
tion particles and that of opposite sign carried 
upwards on ions or cloud droplets, their rate 
of vertical separation cannot exceed the rate 
of fall of the precipitation elements in still air. 
Now, the initial rate of regeneration of mo- 
ment, would if it remained constant, replace 
the whole moment destroyed by the discharge, 
110 C-km, in 7 sec. Taking 8 m/sec. as the 
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fall velocity of the particles, the maximum 
vertical separation of the charges in this time 
would be 56 m; this implies that the magnitude 
of the charges which are separating immedi- 
ately after discharge is about 1,000 C. 

Wilson first suggested that these charges 
would be distributed through a considerable 
volume which is in the main electrically 
neutral, but as separation proceeds through 
the volume, a positively charged region appears 
near the top and a negatively charged region 
near the base. As WORMELL (1953) points out, 
such a bulk charge of about 1,000 C would be 
sufficient without replenishment, to supply 
the charges removed in a considerable number 
of lightning flashes, so that the primary 
charging mechanism need not necessarily 
supply new charge to the region as rapidly as 
it is being dissipated by a brief spell of light- 


ning activity. 


(g) Summary—The Requirements of a Satis- 
factory Theory 
The main features of the thunderstorm 
which have been described in the foregoing 
paragraphs and with which any theory of 
charge generation must be consistent are: 
(i) The average duration of precipitation 
and electrical activity from a single 
cell is about 30 minutes. 


(ii) The average electric moment destroyed 
in a lightning flash is about 110 C-km, 
the corresponding charge being 20 to 
30 C. 

The magnitude of the charge which is 
being separated immediately after a 
flash, by virtue of the fall-speed v of 
the precipitation elements, is of order 
8,000/v coulombs where v is in metres/ 
sec. 

In a large, extensive cumulonimbus 
this charge is generated and separated 
in a volume bounded by the - 5°C 
and the - 40°C levels and having a 
typical radius of say, 2 km. 

(v) The negative charge is centred near 
the - 5°C isotherm while the main 
positive charge is situated some kilo- 
metres higher up; a subsidiary positive 
charge may also exist near the cloud 
base, centred at or below the 0°C 
level. 
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(vi) The charge generation and separation 
processes are closely associated with 
the development of precipitation, par- 
ticularly in the solid form (graupel); 
the precipitation particles must be 
capable of falling through updraughts 
of several metres/sec. 

(vii) Sufficient charge must be generated 

and separated to supply the first 

lightning flash within about 12 to 20 

minutes of the appearance of precipita- 

tion particles of radar detectable size. 


3. The Theories of Charge Generation and 
Separation 


Having reviewed briefly our present know- 
ledge concerning the electrical and meteor- 
ological behaviour of a thunderstorm, the 
two, of course, being intimately related, we 
are now in a position to examine the various 
theories of charge generation and separation 
and determine how far they are consistent with 
the observed facts. 


(a) Influence Theories 


(i) ELSTER and GEITEL (1913) considered the 
effects associated with a raindrop falling 
through a vertical electric field and colliding 
with smaller cloud droplets lying in its fall- 
path. The raindrop will be electrically po- 
larised; in a downwardly directed field such 
as exists in the atmosphere under fine-weather 
conditions, the lower half of the drop will 
carry a positive charge and the upper half a 
negative charge. Elster and Geitel assumed 
that on collision, the cloud particles would 
rebound from the raindrop and carry away 
some of the charge from its lower half. Thus, 
with a positive field, the raindrops would 
acquire a net negative charge, while further 
gravitational separation would enhance the 
original field. Nowadays, this theory appears 
untenable; the coalescence of drops possessing 
differential rates of fall plays an important part 
in modern precipitation theory, and there is 
direct experimental evidence that it takes 
place with high efficiency when drops of 
raindrop-size fall through a cloud of droplets 
of radius about 104 (GUNN and HITSCHFELD, 


1951). 


452 


(ii) Wilson’s Process of Charging by Selective 
Ion Capture 
Wırson (1929) pointed out that under 
certain conditions, an electrically polarised 
raindrop in falling through a cloud of ions or 
charged cloud droplets, could by a process of 
selective ion capture, acquire a net charge. If 
the drop falls more rapidly through a down- 
wardly-directed field than the downwardly- 
moving positive ions, the latter are repelled 
from the lower half of the drop and deviated 
to one side, while the negative ions are attracted 
to it. Hence the drop acquires a net negative 
charge which tends to augment the pre- 
existing field. The mathematical theory of 
the process has been worked out by WHIPPLE 
and CHALMERS (1944) who showed that, if the 
fall velocity of the raindrop is large compared 
with the drift velocity of the ions under the 
existing vertical field, (this will generally be 
the case if the ions are charged cloud droplets), 
(«) the final charge acquired by a drop of 
radius a, independent of its initial charge, in a 
vertical field of intensity X, is 


Qinax Se (3 —2 V2) Xa? = — 0.52 Xa? (1) 


(B) if its initial charge is zero, the initial (i.e. 
maximum) rate of charging of a drop is 


(=) = — 30 Xa°] _ (2) 


where A_ is the polar conductivity for the 
large negative ions, 

(y) the time required for the drop to acquire 
half of its final charge is 


T = 0,04/A_ sec. (3) 


As the above mechanism proceeds, the field 
increases, so that the equilibrium charge 
towards which the drops are tending is in- 
creasing with time. However, we can estimate 
the magnitude of the maximum charge which 
can be generated by this process, and the rate 
of charging. 

Let us suppose that the precipitation process 
has developed to a stage at which the precipita- 
tion rate in the active part of the cloud reaches 
2.5 cm/hr with a corresponding mean volume 
drop radius of 1.5 mm and precipitation-water 
content of 1 g/m?; the concentration of drops 
under these conditions will be 70/m?. Taking 
the average intensity of the vertical field to be 
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1,000 V/cm, the maximum charge that can be 
developed per km? by the Wilson process is 
then 0.9 C. But, we now have to consider 
whether the drops could approach their 
equilibrium charge during the life of a thunder- 
storm cell. There appears to have been no 
reliable determination of polar conductivity 
in cloudy air, but if we substitute a conservative 
value of 4x 107% e.s.u. for A_ in (3), we find 
that it would take 10? sec. for the drop to 
acquire half of its maximum charge. Hence, 
the value of 0.9 C/km? is a considerable 
overestimate of the charge which could be 
built up by the Wilson process during the 
life of an active cell, and even so, is much 
smaller than the value of about 16 C/km? 
required by conditions (iii) and (iv) of Section 
2 (2). 

a WorMELL (1953) points out, the maxi- 
mum possible rate at which the charging 
process can proceed, cannot exceed the rate at 
which ions are being produced in the lower 
atmosphere, i.e. 10 ions cm sec.71,l0r 6 
C-kmhr-?. But the raindrops sweep out 
only a fraction of the space in which ions are 
being produced; in the above example the 
fraction would be only 3.5 x 107%, so that if 
the raindrops catch all the cloud particles en- 
countered in their fall-paths, the charging rate 
would be only 2x 10°? C-km-* hr-!. This is 
to be compared with the average rate of 1,000 
C in 20 minutes generated in a volume of 
about 60 km’, i.e. so C-km~*hr-? required by 
conditions (iii), (iv), (vii). 

It therefore appears impossible that the 
Wilson process can of itself separate the 
observed charges during the life time of a 
typical cell. It may well play a minor role in 
the charge separation process during the later 
stages when copious ionisation and strong 
fields have been produced by the primary 
process. It is almost certainly of importance 
in determining the magnitude and sign of the 
charge on precipitation elements reaching the 
ground after they have fallen through the 
space-charge blanket existing between cloud 
base and ground. 


(iii) Walls Theory 
WALL (1948) does not accept the view that 
the direction in which the electric field builds 
up can be controlled by the normal fine- 
weather field—a fundamental assumption of 
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Wilson’s theory. He also attempts to account 
for the fact that charge separation occurs 
mainly at temperatures below 0°C, by 
suggesting that ice crystals, rather than rain- 
drops, are involved in selective ion capture. 
The initial polarising field is now attributed 
to the polar properties of the hexagonal plate- 
like crystals which, in falling with their princi- 
pal axes vertical, are assumed to become dipoles, 
presumably with their undersides positively 
charged. These crystals, originating at tem- 
peratures below -10°C and having become 
polarised by their own internal field, are then 
assumed to separate the charge in much the 
same way as proposed by Wilson. 

This theory rests on the assumption that 
ice is polar; it should therefore exhibit piezo- 
and pyro-electricity; Wall deduced its polar 
nature, not from direct experimental tests, but 
from observations of asymmetrical develop- 
ment along the c axis of crystals grown from 
the melt. Experimental tests on crystals from the 
melt have all failed to detect piezo- and pyro- 
electricity. Recently, Mason and Owston 
(1952) using very sensitive tests, have failed to 
detect these effects in asymmetric crystals 
grown from the vapour; indeed the molecular 
structure of ice and the randomised motion 
of the dipoles militates against the lattice 
possessing a resultant dipole moment. These 
are serious objections to Wall’s theory which, 
in any case, is subject to the quantitative limita- 
tions of the Wilson theory. It appears to have 
few, if any, redeeming features. 


(iv) Frenkel’s Theory 


FRENKEL (1944, 1946, 1947) considers that 
the electrification of clouds arises primarily 
from the orientation of the molecular dipoles 
at the surface of a water drop or ice crystal, 
giving rise to an electric double-layer with 
the negative ends of the dipoles directed out- 
wards. By virtue of this surface configuration 
the drop (or crystal) is supposed to capture 
negative ions preferentially from an ionised 
atmosphere. In a state of statistical equilibrium, 
the drops are assumed to acquire a net negative 
charge q and a corresponding potential & equal 
to the potential difference across the electric 
double-layer; thus g=aé, where a is the droplet 
radius and £= —0.25 volt. The corresponding 
positive charge is left in the air in the form of 

large ions surrounding the drop. Gravitational 
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separation of the positive ions and negative 
drops then produces a positive field, the growth 
of which is opposed by the migration of ions 
under the field setting up a depolarising cur- 
rent. Frenkel considers that these two effects 
eventually become balanced to produce a 
steady electric field, the total current density 
through the cloud then being zero. The in- 
tensity of the equilibrium field can be written 


4, Àg ote 
X=9 / + = 
N, 8 doa NGeD (4) 


where m, a and v are respectively the mass, 
radius and fall velocity relative to the air, of 
the drop, g the acceleration due to gravity, 
N the number of drops per unit volume 
(assumed all of the same size), and A the total 
conductivity of cloudy air. 

If we take the liquid-water content of the 
cloud to be2 em A= 410° €s'u., 2025 
V "10 -*"e.5.0. werhave that, when g= Jorn, 
X=1,200 V/cm. This value is not seriously 
affected by assuming a higher value of A, e.g. 
3x 10°? e.s.u., the average value for the lower 
dry atmosphere at Kew. 

The time constant of the field is approxi- 
mately BA TGA, ME 0.2102 SCG. OL 82 10 cc 
depending on which of the above values of A 
is adopted. These figures suggest that Frenkel’s 
mechanism, if it occurs, could produce a 
positive field of several volts/cm in the carly, 
non-precipitating stages of the cloud; measure- 
ments by GuNN (1952) suggest that these 
fields are always less than 10 V/cm. 

However, as Frenkel points out, for large 
drops in strong fields, the capture of ions due 
to the polarisation charges on the drop is a 
more important phenomenon. He attributes 
the acquisition of charge by this latter mech- 
anism merely to the difference between the 
polar conductivities due to positive and nega- 
tive ions, the equilibrium charge then being 
positive and of magnitude 


Fe À 
avers ah Xa*=0.13 Xa", (7 13] (5) 
This positive-charging mechanism therefore 
predominates over the negative charging 
process described above if Xa > 6.4 x 107%. 
Thus, a drop of radius 100 u would eventually 
acquire a positive rather than a negative 
charge in a field of 200 V/cm, while the 
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corresponding field for a droplet of 10 y radius 
would be 2,000 V/cm. It then appears that 
while sedimentation of cloud droplets might, 
perhaps, set up an initial positive field by the 
Frenkel process, precipitation particles would 
become positively charged in the fields which 
exist in thunderclouds and that these, in settling 
out, would tend to destroy the original field. 

Frenkel argues that this influence mechanism 
does not affect materially his primary charge 
separation process but may explain the fact 
that rain from thunderstorms is predominantly 
positive in sign. Nevertheless, one would 
expect that, on his picture, the build up of the 
positive field would reach a maximum in the 
Cu-congestus stage, before an appreciable 
fraction of the liquid water becomes trans- 
formed into precipitation elements which will, 
according to his influence mechanism, oppose 
its growth. This is inconsistent with the 
observation that the development of precipita- 
tion is essential to the onset of lightning 
activity. Moreover, as the fall velocities of 
cloud particles are so small, the charge which 
must be separated in order to produce a 
lightning flash must be correspondingly large. 
For droplets of a=1ou, for which v=ı 
cm.sec.-!, generation of an electric moment 
of 110 C-km in 7 sec. would require separa- 
tion of 8x 10°C of charge. The maximum 
negative charge which a 10 u droplet will 
aequite by Frenkel’s process is 10° €.s.u.; 
assuming a concentration of 500 droplets/cm?, 
we find that is a volume of 60 km? the maxi- 
mum possible charge would be 10° C, but the 
actual charge built up during the life of the 
cell will be an order of magnitude less. Thus, 
quantitatively, Frenkel’s theory is incapable 
of accounting for the rate of charge separation 
required in a thundercloud. With regard to 
separation of charge leading to a positive 
field, the maximum charge predicted by 
Wilson’s theory exceeds that given by Frenkel 
for drops of a > sojmin a field of 100 V/cm. 


(b) The Drop-Breaking Theory 


The electrification associated with the 
rupture of water drops on colliding with a 
solid surface was first investigated in detail by 
LENARD (1892). SIMPSON (1909), in a series of 
careful experiments in which he largely elim- 
inated this “splashing” or “Lenard” effect, 
established that breaking of a drop in a strong 
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vertical air jet could also produce a considerable 
electrification. The results were found to be 
sensitive to traces of impurity in the water, but, 
for distilled water the large fragments of the 
broken drops carried positive charges, while 
the surrounding air contained ions of both signs 
with an excess of negative charge. Ruptured 
drops of diameter about 8 mm produced, on 
average, a charge of 5.5 x 10°? e.s.u., i.e. about 
2.3 x 1072 e.s.u./cm?; with less violent shatter- 
ing in conditions approximating more closely 
to those under which a raindrop might break- 
up in a thundercloud, the charge produced was 
slightly less — 1.5 x 10°? e.s.u./cm?. 

ZELENY (1933), using highly purified water, 
found that the rupture of drops in a 20 m.sec.-! 
borizontal air jet was accompanied by charges 
of about 2x 1072 -e:s.u./cm°. The negative 
charge communicated to the air increased 
quite rapidly with increasing velocity of the jet. 

CHAPMAN (1952) allowed drops of distilled 
water of diameter 4 mm to fall into a vertical 
jet of speed 17.3 m.sec.-!. The violent disrup- 
tions were accompanied by separation of 
charge amounting on average to about 0.3 
SU. DER drop, 1.6: 10 ¢.s;u./em* However 
when he allowed two drops to coalesce and 
the resultant unstable mass to break-up into 
large fragments in a steady upcurrent of 8 
m.sec.-1, the charge generated was less than 
one-thousandth of the previous figure, i.e. 
of the same order. of magnitude found by 
Simpson and by Zeleny. 

It then appears, that the charge separated 
by breaking drops depends markedly on the 
violence with which they are shattered, i.e. 
the number of fragments which result. Accord- 
ing to Simpson and Zeleny, the more pure the 
water the higher the charge generated; Chap- 
man, however, found evidence for a slight 
increase in the charge when solutions of con- 
centration 5 x 10-4 N. were used. 

It is generally accepted that falling raindrops 
cannot exceed a diameter of about 5.5 mm. 
At this size they become distorted from the 
spherical, become unstable and break up into 
several large, and probably many more small 
fragments. These conclusions are based on 
the studies of LENARD (1904) who balanced 
drops of various sizes on vertical air jetsso 
appropriate velocities. Some recent work by 
BLANCHARD (1950) has shown that the size at 
which drops rupture depends largely on the 
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scale and intensity of the microturbulence in 
the jet, and that the maximum stable diameter 
may be greater, or less than 5.5 mm, depending 
on these factors. 

It is likely that a raindrop in falling through 
an updraught region of a cumulus will en- 
counter small-scale eddies and be subjected 
to appreciable accelerations, but probably not 
as violent as those obtaining in Chapman’s 
experiments with the high-speed jet. Even the 
accelerations in the experiments of Simpson 
and Zeleny appear severe compared with 
those which might be experienced by a 
raindrop in a thundercloud, where, according 
to the Thunderstorm Project, velocities of 
sharp-edged gusts rarely exceed 8 m.sec.-1. 

We shall therefore assume the values quoted 
by Simpson and by Zeleny of about 2 x 10-2 
e.s.u./cm as an upper limit for the charge 
separated by the rupture of water drops in a 
cumulonimbus, and enquire whether this 
could be a significant factor in the generation 
of thunderstorm electricity. 

It is clear that if a positive charge resides on 
the large fragments of the broken drops and 
a corresponding negative charge is communi- 
cated to the air by ions (which may soon 
become captured by cloud droplets), gravita- 
tional separation will tend to confer a negative 
polarity on the cloud, i.e. of opposite sign to 
the observed polarity. It is of interest, however, 
to determine whether the charges generated 
by drop-breaking are significant in magnitude 
compared with those involved in lightning 
discharges. 

We may take 4 gm as an average value for 
the liquid-water content of a large cumulus, 
and if this is all involved at least once in drop 
breaking, the charge generated would be 
3 x 10-2 C/km?, assuming Simpson’s value of 
2.3 x 10-2 e.s u./gm. It is clear, then, that for 
the process of drop-breaking to be significant 
in the electrical budget of the thunderstorm, 
repeated rupture of the same mass of water 
would be necessary. Now, if one traces the 
growth of a raindrop by the coalescence 
process in a cumulus cloud having a base 
temperature of 10°C, a steady updraught of 
8 m.sec.-! and containing the theoretical 
maximum liquid-water content, one finds 
that it is theoretically possible for it to reach 
a diameter of 5 mm before reaching cloud base. 
If, at this stage, the drop is assumed to break 
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into fragments of radius 1 mm, the minimum 
time taken for each of these to grow again 
to break-up size by coalescence, is about 4 
minutes. This can be seen from the following 
simplified and approximate calculation. The 
growth rate ofa drop by coalescence is given by 


da/dt= Ewv/40 (6) 


where a,0,E are respectively the radius, 
density, and collection efficiency of the drop, 
v its fall-speed relative to the cloud droplets, 
and w the liquid-water content. Thus assuming 
E and w to be constant during the growth 
from a=0.I cm to 0.25 cm, we have 


where ¢ is the time required for growth to 
break-up size. Putting E=0.9, w=4x 10-8 
g.cm-*, g=1, and v the average fall speed 
during growth=700 cm.sec-!, f turns out to 
be 239 sec. Actually, this is an underestimate 
of the time required, because when there are 
a considerable number of growing rain-drops 
all competing for the available cloud droplets, 
the effective value of w will be less than that 
assumed here. 

We have seen that only about 12 minutes 
elapses between the appearance of the initial 
radar echo from a thunderstorm and the 
first lightning flash, so that there is time for 
about only three successive ruptures of the 
same mass of water during this interval. The 
maximum possible charge which will be 
generated in’ 1 km® is therefore-9x 10-2 CG; 
which is two orders of magnitude smaller 
than that required by (iii) and (iv) of Section 
2 (g)". 

It has been suggested by Simpson and others 
that the drop-breaking mechanism might 
account for the appearance of the pocket of 
positive charge in the base of some thunder- 
clouds, generally in the region of heaviest 
precipitation. Simpson describes the positive 
charge as having magnitude about 4 C dis- 


1 Tt must be pointed out, however, that no experi- 
ments have been reported which provide information 
on the charges liberated by drops breaking in the 
presence of a strong electric field. In these circum- 
stances the separated charges might be considerably 
larger than those measured by Simpson and by Zeleny 
in view of the polarisation charges 3 Xa?/4. 
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tributed throughout a volume of about 0.5 
km, ie. 8 C/km?. This again is about 100 
times larger than we should expect to be 
generated by the drop-breaking process, so 
that it appears necessary to find another mech- 
anism to explain the origin of the lower 
positive charge centre. It may well arise by 
capture of positive ions by the drops from 
the upwardly-directed point-discharge current 
which prevails beneath thunderclouds. 


(c) Simpson’s Theory of Charging by Collision 
of Ice Crystals 

SIMPSON and SCRASE (1937), SIMPSON (1942) 
suggested that the main charge structure of a 
thundercloud could be accounted for by the 
production of frictional electricity during the 
collision of ice crystals which were assumed 
to receive a negative charge, the compensating 
positive charge being carried upwards on the 
cloud droplets. 

Unfortunately, there is practically no direct 
information as to the sign and magnitude of 
the charge released by colliding ice crystals. 
Pearce and Currie (1949) have eroded a 
block of snow with an air blast and found the 
large eroded fragments to carry a negative 
charge, the air receiving a positive charge. 
CHALMERS (1952) has rubbed two handfuls of 
snow together and allowed the fragments so 
produced to fall into a collector. On every 
occasion a negative charge was recorded, the 
generation of which he attributed to friction. 
The results of these experiments, cannot how- 
ever, be applied quantitatively to the atmos- 
pheric problem, and since we can make no 
reliable estimate of the frequency of collisions 
between ice particles in a turbulent cloud 
mass, it is not possible to demonstrate that 
such a mechanism could generate and separate 
charge at the required rate. 

There is, however, some laboratory evidence 
to show that when air currents of a few 
cm/sec. are allowed to flow past a deposit of 
frost grown by sublimation, small splinters 
are broken off the delicate dendritic crystals 
and carry away charges predominantly of 
one sign, leaving those of opposite sign on the 
parent crystal. Such a splintering process might 
well occur in clouds when snowflakes com- 
posed of delicately branded crystals collide or, 
they may be torn off by frictional drag with 
the air. FINDEISEN (1940, 1943) found the 
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splinters from a growing deposit to be pre- 
dominantly negatively charged and the deposit 
positively charged; with an air stream of 35 
cm.sec.-! the rate of charging of the deposit 
was 4 x 10715 C/cm? sec. KRAMER (1948) found 
a charging rate of the same magnitude and 
sign in the early stages of a growing frost layer 
but found that later, the polarity became 
reversed. Kum (1951) measured the charges 
on individual splinters and found positive 
splinters to be about 7 times as numerous as 
negative ones, so that the frost deposit acquired 
a negative charge. The average charge carried 
per splinter was much the same as found by 
Findeisen, i.e. about 10-15 C. 

Let us take Findeisen’s value of 4 x 10-16 C 
cm-? sec.-! for the rate of charging and make 
a crude estimate of the rate of charge produc- 
tion by splintering in a cloud. We shall assume 
a frozen-water content af 1 gm”? composed 
of snowflakes of average mass ı mg and effec- 
tive area 1 cm?, (a high value). The charge 
generated per second in a volume of ı km? 
would then be 4x 10-4 C, ie. about 0.5 C 
would be produced over a period of 20 
minutes. If the fall-velocity of the flakes is 
assumed to be I m.sec.-!, generation of an 
electric moment of 110 C.km in 7 sec. would 
require separation of 8,000 C of charge, or 
about 130 C per km?. Thus, unless the splinter- 
ing of ice crystals in the atmosphere is more 
efficient than in the above laboratory experi- 
ments, by at least two orders of magnitude, 
it cannot play a significant role in the genera- 
tion of thunderstorm electricity. Furthermore, 
it is doubtful whether dendritic snow crystals 
form a major part of the frozen-water content 
during the active growth stage of the thunder- 
cloud, when falling crystals collect large 
numbers of supercooled droplets which freeze 
on contact to form pellets of graupel or hail. 


(d) The Dinger-Gunn Effect 


Dinger and GUNN (1946) have shown that 
when ice containing air melts the air bubbles 
are released, and upon breaking the surface 
they transfer to the adjacent air a negative 
charge, while the melted water retains an 
equal and opposite positive charge of magni- 
tude 1.25 e.s.u./gm. This process may have 
relevance during the melting of graupel pellets 
and snowflakes. It is of the wrong sign, how- 
ever, to account for the observed polarity of 
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thunderclouds, but may be partly responsible 
for the predominance of positive charge 
associated with steady rain and may, perhaps, 
contribute to the local positive charge found 
in the bases of warm thunderclouds. In a 
cloud where solid precipitation elements in a 
concentration of 2 gm? are approaching the 
0° C level, a charge of 0.8 C.km-* would be 
developed according to the results of Dinger 
and Gunn. 


(e) The Workman-Reynolds Theory 


Workman and REYNOLDS (1948, 1950 b) 
have made extensive investigations of the 
electrical effects associated with the freezing 
of water and dilute aqueous solutions. They 
discovered that, during freezing, a potential 
difference developed across the ice-liquid 
interface, the sign and magnitude of which 
depended on the nature and concentration of 
the solute. For the majority of solutions tested 
the ice became negative with respect to the 
liquid, with important exceptions in the 
case of ammonium salts. The largest potential 
difference, —232 V, (potential of liquid with 
respect to ice), was obtained with a 5 x 1075 N 
solution of ammonium hydroxide. Sodium 
chloride and calcium carbonate solutions of 
concentration 10-4N produced potential dif- 
ferences of +30 V and + 20 V respectively, 
while doubly-distilled water carefully freed of 
traces of ammonia showed practically no electri- 
cal effects. It appears that, during the freezing 
process, solute ions of one sign are preferentially 
incorporated into the ice, those of opposite sign 
being left in excess in the liquid. 

Workman and Reynolds have attempted 
to erect a theory of charge generation and 
separation in thunderstorms on the basis of 
these laboratory experiments. The kernel of 
their arguments is the wet hailstone—an ice 
pellet which, on reaching a critical size and 
fall velocity determined by the temperature 
and liquid-water content of the cloud, com- 
mences to collect supercooled cloud droplets 
at a rate faster than these can freeze, and so 
acquires a liquid coat. The limiting factor is 
the rate at which the latent heat of fusion can 
be dispersed to the environment, the theory 
of which has been worked out by LupLam 
(1950). Workman and Reynolds assume that 
hailstone will first become wet at temperatures 
between - 10°C and —15° C and that, there- 
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after, only a fraction of the impinging water 
will become frozen, the rest being flung off 
in the form of small drops. It is further assumed 
that at the ice-water interface, negative ions 
will be preferentially incorporated into the ice, 
(as in the laboratory experiments with solu- 
tions of NaCl and CaCOs), so that the water 
drops flung-off will carry away a positive 
charge leaving the hailstone with a net nega- 
tive charge. Gravitational separation of the 
negative hailstones and small positive drops 
will then lead to a charge distribution of the 
observed polarity. In the non-supercooled 
region of the cloud, the water shed by the 
hailstone will carry away its excess negative 
charge and the authors suggest that these 
negative drops may be carried up to levels 
where they may be captured by negative 
hailstones, whose charge is thereby enhanced. 
The build up of the electric field thus becomes 
a cumulative process. 

The two main advantages which the authors 
claim for their theory is that it affords a natural 
explanation for the centre of the negative 
charge being located around the — 10° C level 
and that the charges transferred across the 
ice liquid interface during the laboratory 
freezing experiments (a maximum value of 
9 x 10* e.s.u. per cm? of frozen NaCl solution) 
are of sufficient magnitude to account for the 
charges appearing in thunderstorms. 

The assumption that hailstones will not 
become wet until they approach the - 10° C 
level may be challenged. For a cloud of given 
temperature and liquid water content, the 
temperature at which an ice particle becomes 
wet depends largely on the density with which 
it is assumed to grow. Calculation shows that 
in a cloud with base temperature 10° C, the 
maximum theoretical water content and a 
steady updraught of 10 m.sec.-1, an ice pellet 
originating at —5°C and growing with a 
density of 0.3 will first become wet at - 25° C 
—i.e. 10 !/, minutes later. 

Whether the Workman—Reynolds effect 
has relevance to thunderstorm electrification 
will depend upon the contaminants present 
in the precipitation elements being able to 
give a positive water-negative ice charge 
distribution during freezing. WORKMAN and 
REyNOLDS (1951) collected small hail pellets 
from near the 0° C level in an active thunder- 
storm and found them to contain CaCO; 
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and NaCl, the concentrations of which were 
not, however, determined. The electrical 
effects associated with the re-freezing of the 
melted particles were erratic. When stored at 
room temperature, the water became negative 
as freezing began and assumed zero potential 
at the end of the freezing period. When 
stored at 0° C before freezing, the water again 
assumed a negative potential when freezing 
began, but became positive after about 20 
sec.; thereafter the potential fell to zero at the 
end of the freezing period. These erratic 
effects were attributed to the absorption of 
CO, and its reaction with the calcium car- 
bonate to produce soluble calcium bicarbonate. 
Certainly the electrical effects obtained with 
solutions of CaCO, and Ca(OH), were very 
sensitive to the amounts of dissolved CO. For 
NaCl solutions the charging phenomena were 
very sensitive to the concentration of salt; 
while maximum charge separation was ob- 
tained with a 10-4 N solution, no effect at all 
was produced when the concentration was 
increased to 5x 10-4N. With doubly-distilled 
water the water became negative with respect 
to the ice during freezing, an effect which has 
been confirmed by Grit and ALFREY (1952) 
Workmann and Reynolds attribute this to the 
presence of small traces of ammonia. 

The fact that the Workman—Reynolds 
effect appears to be so sensitive to the con- 
centrations of salts and carbon dioxide and to 
traces of ammonia in the cloud water makes 
it unattractive as the main mechanism of 
charge generation in thunderstorms. Also, an 
ice particle becomes wet only when it has 
attained a radius of several millimetres; it 
then acquires a clear coat of ice to form a true 
hailstone and takes longer to melt than the 
low-density graupel. This leads us to suppose 
that if the Workman—Reynolds effect were 
the most important charge generating mech- 
anism, thunderstorms would generally be ac- 
companied by hail. The observational evi- 
dence, although admittedly not well docu- 
mented, does not show this to be the case. 

The cumulative build-up of the electric 
field in the manner visualised by the authors 
also seems open to doubt. It seems likely that 
most of the water accumulated by a wet 
hailstone will be flung off from the top-side, 
in the form of large drops rather than as small 
ones; thus their velocity of separation from the 
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hailstones and the probability of their capture by 
hailstones at higher levels will be correspond- 
ingly reduced. Furthermore, in the presence of 
a positive field, the separation of liquid water 
from the negatively polarised upper half of 
the hailstone will tend to diminish the field. 


(f) Charge Separation Associated with the For- 
mation of Rime 


A number of experimenters have, in recent 
years, reported that when supercooled droplets 
are allowed to impinge and freeze on an ice 
surface, the latter acquires a substantial charge. 
This effect is to be distinguished from that 
which obtains when non-supercooled water 
splashes off an ice surface. In the latter case it 
was shown early on by Faraday and later by 
Sohncke that the ice becomes positively 
charged, the rebounding water drops acquiring 
a negative charge; the same result has been 
obtained recently by Girt and ALFREY (1952). 

FINDEISEN (1940, 1943) formed a rime layer 
by spraying water droplets onto a cold metal 
surface and found that it acquired a positive 
charge of order 3)x 10.1? Gem sec ie 
charging ceased if the surface became smooth 
and glassy, or if it became wet. 

A very careful re-investigation of these 
phenomena has been made by Kramer (1948) 
who found the ice deposit acquired a negative 
charge which varied as the impact velocity of 
the drops; for a speed of 0.5 m.sec.-1 the 
charging Late was 52.2102. ©. em esc rmee 
Findeisen assumed that the compensating 
charge was carried away on splinters in much 
the same way as for his frost deposits, (see 
Section 3 (c)), and some evidence for the 
production of charged splinters during riming 
was obtained by Kramer. 

LUEDER (1951) made experiments in natural 
supercooled clouds on a mountain top in 
order that the contaminants present in the 
water should be those occurring in nature. 
The rime was deposited on a slowly rotating 
metal rod which was enclosed for a part of 
each revolution by a Faraday cage. Its charge 
caused a periodic variation in the potential of 
an electrode. Lueder states that the growing 
deposit of rime acquired a very strong negative 
charge, an equal positive charge being com- 
municated to the air, probably on the parts 
of the drops which were flung off without 
freezing. The charge generation was more 
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pronounced at lower temperatures, i.e. the 
more efficient and rapid the freezing process. 
It is not, however, easy to deduce the actual 
magnitude of the rate of charging from these 
experiments. 

MEINHOLD (1951) measured the electric 
field strength at the surface of the fuselage of 
an aircraft flying at 80 m.sec.-! through a 
supercooled Cu-congestus. The deposition of 
rime was accompanied by a rapid rise in the 
field strength from about 200 V/m to 5,500 
V/m in a sense which showed that the aircraft 
skin was acquiring a strong negative charge. 
From the rate of increase of the field intensity 
Meinhold calculated the rate of charging to 
be 5 x 10-12 C.cm-? sec.-! and attributed this 
high rate to the high speed of impact and the 
rapid freezing of the droplets. 

The balance of the evidence from these 
experiments points to the acquisition of a 
negative charge by a growing layer of rime, 
the rate of charging being higher the velocity 
of impact of the drops and the more efficient 
the freezing process. The outstanding contra- 
diction is provided by Findeisen’s experiment 
in which he claimed the rime deposit acquired 
a positive charge. Some difficulty may arise 
in the interpretation of the experimental 
results because of the complications connected 
with the rebound of some of the droplets 
from the ice—a process which appears to 
communicate a positive charge to the ice. It 
seems important, therefore, to separate the 
electrical effects associated with the freezing 
and with the rebound of droplets in these 
investigations. 

In this respect the most satisfactory and 
conclusive experiments have been carried out 
by WEICKMANN and aufm KAMPE (1950 and 
private communication). Their cold body was 
a metal rod of diameter s mm placed in a 
coldoom kept ‘eitherrat»— 5° C ors—12° C. 
Water droplets of diameters 5 4! to 100 u were 
sprayed on to the rod at velocities varying 
from 5 to 15 m.sec.-! The droplets in the 
immediate vicinity to the rod were slightly 
supercooled and the water content of the 
cloud in this region was about 4 g.m~*. Thus, 
conditions relevant to graupel pellets or small 
hailstones falling in a cumulonimbus were 
fairly well simulated. Singly-distilled water, 
tap water, weak NaCl solutions and strong 
solutions of ammonium hydroxide and sodium 
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fluoride were used in different experiments, 
but in each case the rate of charging of the 
rod was the same, for a given air velocity. 
The charge generation increased with the 
velocity of the air-stream, i.e. the rate of 
accretion of the rime, and for a velocity of 
Is m.sec.-! attained a value of 5 x 10-12 C. 
cm~*sec.-' The ice deposit on the rod showed 
the clear and opaque layers characteristic of 
hailstones, but charging stopped when the ice 
became wet. When water at temperatures 
slightly above freezing was sprayed on the 
rod it acquired a slight positive charge. 


4. Conclusions 


We have seen that all the theories of charge 
generation and separation which have been 
advanced during the last forty years, and 
which have been discussed above, are open to 
objection on quantitative grounds and/or 
because they do not fit the known facts about 
the meteorological and electrical behaviour of 
thunderstorms as listed in Section 2 (g). The 
only other mechanism suggested by laboratory 
experiments and which appears to fit, in a 
qualitative manner, the observed facts, is the 
charge generation associated with riming. The 
onset of lightning activity appears to coincide 
with the appearance of heavy solid precipita- 
tion, especially in the form of graupel. We 
recall that Malan and Schonland gave evidence 
for the negatively charged column extending 
up to, but not beyond the — 40° C level; this 
is consistent with the theory of charge genera- 
tion accompanying the growth of graupel 
particles because supercooled droplets may be 
found at temperatures down to, but not 
below -40°C. The fact that laboratory ex- 
periments show that during the formation of a 
rime layer, a negative charge is acquired by the 
ice and presumably an equal positive charge 
by the air, indicates that the growth and fall- 
motion of graupel pellets would bring about 
a charge distribution of the observed polarity. 

It remains to be shown, quantitatively, that 
this process is capable of generating and sepa- 
rating the required quantity of charge in the 
time suggested by the electric field-change 
and radar data. This will form the subject of 
a future paper.? 

2 See Mason, B. T. “On The Generation of Charge 
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Abstract 


Investigations into the O18-abundance in the precipitation from a warm front reveal a 
distinct process of fractionation explainable by the unequal pressures of the saturated vapours 


of H,O# and H,0}8, 


An explanation is given of the low O!8-abundance in glacial waters. 
A method of experimental investigation into the mechanism of the process of precipita- 


tion is discussed. 


Research upon the abundance in nature 
of the heavy oxygen isotope, O18, has been 
carried on persistently for several years. 
Before World War II the most important 
method of measurement was by density tests 
made on water. Of late, however, excellent 
mass-spectrometers have become available, 
and these have made it possible through a 
relatively simple measuring technique to 
achieve an equally good reproducibility as by 
the density tests; and it is possible, moreover, 
by a suitable procedure to obviate the cor- 
rection for the deuterium content in the water 
—a not inconsiderable correction, which is 
lacking in several of the densicy tests published 
in the relevant literature. 

The O!$ tests on water, recorded in this 

aper, were made with a mass-spectrometer 
of the Nier type (1). For technical reasons, 
e.g., owing to the vacuum, it is not possible 
to introduce the water directly into the appa- 
ratus. It is necessary, therefore, either to electro- 
lyze the water and introduce the oxygen—a 
process which is not quite satisfactory, because 
of the undesirable action of the pure oxygen 
on the filament of the ionic source—or to 
shake the excess of water with carbon dioxide, 
as in the tests here under discussion. According 
to ConN and Urey (2), oxygen atoms are 
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exchanged between CO, and H,O, to the 
effect that equilibrium 


CO; +H,O# = CO18O18 + H,O18 


is achieved after the lapse of 3 hours. Next, 
the O18abundance in the carbon dioxide is 
measured and made equal to the Ot$-abun- 
dance in the water multiplied by the coefficient 
of equilibrium of the process stated above, 
which at 25°C is equal to 1.080. From this 
the abundance of O18 in the water is computed. 
In similar manner the abundance of Ol!” may 
be found, though this is rarely of any signifi- 
cance. It is proved that the abundance of the 
three oxygen isotopes 


O16 : O17 : Ol8 = 99.76 : 0.04 : 0.20 %, 


depending to some extent on the origin of the 
water. Technically it is of importance that 
one never is faced with the need to measure 
these absolute values directly but only to 
determine differences, e.g., between O1- 
abundances in two different samples of water. 
The absolute values can be found by deter- 
mining the difference between the sample 
and a standard of known O18-abundance. 
Difference tests on samples with O18-abundance 
approximating the normal value of 0.20 % 


0/ 


show a reproducibility of + 0.0001 %. 
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Dore, VROOMAN, a.o. (3, 4, 5) have deter- 
mined the relative OY-abundances in atmos- 
pheric air and water of various origins. The 
tests are set out in Table 1, where the O18- 
abundance in the fresh water from Lake 
Ontario is made equal to 1.000. 


Table ı 


O!8-abundance 


Ref. 
relative] % 

Water: 
Fresh water, Lake Ontario] 1.000 |0.1972| (5) 
Atmospheric water vapour 

from above Lake Ontario} 0.991 | 0.1954 (5) 
Glacier water, Lake Louise] 0.977 | 0.1927| (5) 
Dead Sea water ........ 1.020 | 0.2011] (5) 
Oceans Wateranec er 1.009 | 0.1990} (4) 

Air: 
Atmospheric oxygen...... 1.033 |0.2033| (3) 
Atmospheric carbon dioxide] 1.040 | 0.2053] (5) 


The variations in O!%-abundance, given ex- 
pression in Table 1, are explainable principally 
through the fact that the ratio of the pressures 
of the saturated vapours of H,O"% and HO}, 
pi$ and p18, is more than 1. By experiment, 
RIESENFELD and CHANG (6) have determined 


px /p% as: 
log1o(po [po ) =2.74/ T - 0.0056 (1) 


where T is the absolute temperature. In Fig. 1 
(x) is plotted in a graph. At the temperature of 
0° C pi6/pt8 = 1.01. Conversely pi/p} is within 
the same temperatures approximately 0.99. 
This means that the O!%-abundance in vapour 
is obtainable ay multiplying with 0.99 the 
O*-abundance in the water in equilibrium 
with the vapour. This directly explains, 


253 258 263 268 273 278 283 268 293 298 is 


Fig. 1. The ratio of the pressures of the saturated vapours 
of H,O1® and HO, po6/pot8, plotted against the 
absolute temperature T. 
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first, why the O18-abundance in the water 
vapour above the Lake Ontario is about 1 % 
less than in the water from the lake proper; 
second, that in fresh water, mainly a con- 
densation of ocean water vapour, the O18 
abundance is about 1 less % than in ocean 
water; and third, that it is a fact, obviously, 
that the Dead Sea, where the affluent water 
is exclusively compensated by the high rate 
of evaporation, will attain a very high O!%- 
abundance. On the other hand, the low 
value for the Glacier water is not immedia- 
tely explainable in a similar way. 

Information on precipitation is sparse. 
Harapa and Tirant (7) have recorded that by 
observations on rains of long duration the 
density was higher in the beginning than at 
the end, and that snow-water is somewhat 
lighter than rainwater. RIESENFELD and CHANG 
(8) explain this phenomenon by equation (1), 
assuming that the O18-abundance in the snow 
is equal to the O18-abundance in the vapour, 
and that the latter value is dependent on the 
temperature of evaporation only. However, 
the variations mentioned in this paper with 
reference to the O!8-abundance in the pre- 
cipitation exceed by several times the magni- 
tudes that can be explained exclusively from 
the assumptions made by RIESENFELD and 
CHANG. 

Teıs (9) found the following variations in 
the density of rain falling within the same 


period: 


Time cc A d(y/cc) 
12bo0— 1315 75 +0.4 
1321 5—14530 150 — 2.3 
14530—1545 75 a 
ısb4s—ı8hıs 60 — 2.3 
185] 5—19®15 40 — 0.3 


Although it is not possible to distinguish 
between the H2- and O18-variations the results 
indicate a decrease in the O18-abundance in 
the beginning and in contradistinction to the 
results obtained by Harada and Titani, and 
increase towards the end of the period. 

A certain kind of precipitation frequently 
occurs at the so-called warm front. So far, 
no explanation in detail of the mechanism of 
this process can be given with any certainty. 
However, in principle, it so occurs that warm 
and humid air (Fig. 2) flows upwards over 
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Fig. 2. Vertical section perpendicular to a warm front. 


relatively cold air and is chilled so that some 
of its vapour is condensed. In the case of a 
cold front the cold air glides forward under 
the warm and humid air, which is carried 
aloft and chilled so that some of its vapour is 
condensed. 

On June 22nd, 1952, a very strong warm 
front passed over the city of Copenhagen. In 
advance, some showers had fallen, and a 
cold front followed immediately upon the 
passage of the warm front. In the upper part 
of Fig. 3, plotted along the horizontal axis, 
are shown the atmospherie conditions pre- 
vailing at those times and in those places where 
the rainwater samples used in the present study 
were gathered. In the lower part of Fig. 3 
the sections of the horizontal line indicate 
the time intervals within which precipitation 
occurred. The horizontal line-sections are set 
out above the time-axis at levels corresponding 
to the percentage excess of O in the precipita- 


6 12 18 0 6 12 18 0 6 12 18 
Jun. 21 22 23 


Fig. 3. The upper part shows weather conditions and 
structure of the atmosphere at the place and time of 
collections of rainwater for analysis. The horizontal 
line-sections in lower part of the figure indicate the 
time intervals during which precipitation occured. 
These sections are placed above the time-axis at distances 
which correspond to the percentage excess of OT in 
the corresponding rainwater sample as compared with 
a certain standard. 
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tion in question compared with a certain 
standard. It will be seen that the rain showers 
which preceded the warm front show an 
excess of O18 of about 0.0023 %,, whereas the 
rain from the upper part of the warm front 
is considerably lighter. Through the front 
(proceeding to the right in the figure) the 
O}8-abundance of the precipitation (hereinafter 
designated ay) increases, and the rain falling 
at the point where the warm air commenced 
its upward motion is heavier even than the 
preceding showers of rain. On the arrival of 
the cold front a similar phenomenon occurs. 
However, because of the violent turbulence 
developed here, the conditions are so compli- 
cated that quantitative discussion is out of the 
question. 

Following the movement of the air upwards 
through the front, the drop in a, (to the left 
in the figure) can be qualitatively explained 
through equation (1): The first rain, precipi- 
tated after the upward motion of the air had 
begun, according to (1), has a higher O18- 
abundance than the vapour from which it was 
condensed. Hence, the remaining vapour has 
become lighter and the rain condensed there- 
from is consequently lighter again than the 
first rain, and so on. There is, however, a 
slight increase in a, in the rain originating 
from the topmost lavers of the warm front. 
This might possibly be attributed to the fact 
that this rain had been subject to vigorous 
evaporation during the relatively long fall 
through the cold air. 

This fact also goes to explain qualitatively 
the low figure for O18-abundance in the 
glacier water, Table 1, as it may be assumed 
that this water exclusively came from high 
altitudes above the ocean surface, thus cor- 
responding to the topmost part of the warm 
front. 


This interpretation seems to be more likely 
than an earlier statement made in the literature 
where the low O18-abundance is ascribed to 
an isotopic fractionation assumed to take place 
in the melting of the ice during which “the 
Oï6-rich molecules, being the most energetic 
of the possible isotopic combinations, would 
be prevalent in the early fractions of the 
liquid phase” (Sırverman (10)). Regardless of 
the fractionation that may occur during the 
first part of the melting the water, for a 
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considerable period of time, must average an 
O18-abundance equal to or higher than that 
of the precipitation, since accumulation of 
H,O1$ can not take place. The only reason 
possible for any difference between these two 
values must be ascribable to evaporation 
which would cause the O18-abundance in the 
melting water to rise above that of the ice 
and the precipitation. 

The explanation put forward in this paper 
find support in measurements of the O18- 
abundance in six chunks of ice from Green- 
land glaciers. These chunks of ice, picked 
up from the sea, originated from the inland 
ice-cap, where the precipitation presumably 
falls at temperatures down to — 50°C. The 
O18-abundances relative to that of fresh water 
from lake Ontario were found be in the 
vicinity of 0.977 to 0.984 in good agreement 
with the value of glacial water given in 
Table 1. The ice investigated was taken from 
the central part of the chunks, which should 
preclude any possibility of exchange with sea 
water. Such exchange, otherwise, would con- 
tribute to an increase in O18-abundance, 
which would also take place in the event of 
evaporation. Hence, the cause of the extrem- 
ely low values must be looked for, exclusively, 
in the very great loss of moisture occurring 
in the air previous to the formation of the 
precipitation involved. 

Before attempting a quantitative explana- 
tion of these conditions, it is necessary to 
consider more closely the mechanism of the 
warm front process: In the process of con- 
densation primarily cloud droplets of diameter 
10 u are formed. From these droplets and/or 
from the vapour, rain drops of diameter 1 
mm are formed. Observations indicate that 
at least in our latitudes precipitation is not 
likely to occur until the clouds have reached 
an altitude where the topmost part is at a 
temperature below - 10° C. Consequently, 
BERGERON (11) has assumed that at this tem- 
perature some ice-crystals are formed which, 
on account of the supersaturation of the sur- 
rounding vapour, would rapidly increase in 
size at the expense of the water droplets and in 
consequence drop down. There are, however, 
rain-clouds which do not extend above the 
o°-isotherm, showing that rain drops can be 
formed in other ways than from ice-crystal- 
nuclei. FINDEISEN (12) and Lancmurr (13) have 
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investigated the possibility of producing rain 
drops through collisions of cloud droplets. 
If the droplets are not of equal size they are 
likely to collide in falling because of their 
differing velocities. It is probable, however, 
that a certain minimum size is required for a 
collision to result in the merging of two drops. 

With regard to the first mentioned theory 
Houcuton (14) has calculated that the time 
required for an ice-crystal of diameter 20 y te 
grow by sublimation to a diameter of only 
200 u would be 5 to 10 minutes in the most 
favourable circumstances, and that it would 
require several hours for its growth to the 
size of a rain drop (diameter 1 mm). However, 
the same investigator (15) has also calculated 
the time required for a droplet of diameter 
30 u to grow exclusively by collision in air 
containing 20 «-droplets of the order of 1 g/m?. 
The result shows that 105 min. are required 
for the droplet to grow frow 30 u to 20 4,0 
and only 18 min. to grow from 200 u t theo 
diameter of 1 mm. From these facts Houghton 
concludes that it is most likely that the process 
of sublimation predominates in the beginning, 
and that the process of collision gradually 
assumes a greater role as the droplet or the 
erystal increases in growth. Since the volume 
increases directly with the third power of the 
diameter it may be concluded, further, that 
the greater part of the volume of water which 
reaches the ground as rain is formed by the 
process of collision. 

On this basis, a picture of the formation of 
precipitation of the warm front, in best 
agreement with our present knowledge, may 
be drawn. Thus, AB in Fig. 2 represents the 
boundary between the warm air and the 
cold air, the warm air crossing from the right 
side upwards to the left side. When the cloud 
is chilled to slightly below the dew point 
(the humidity being about 102 %), the 
formation of cloud droplets begins. Apart 
from a very slow settling velocity due to 
gravity the cloud droplets will follow the 
course of the air. During the following 
upwards glide an increase in the number of 
cloud droplets occurs as a result of the con- 
tinued cooling of the vapour, while a decrease 
in number occurs because some of the droplets 
either evaporate on to, or collide with, the 
heavier particles falling from above (ice- 
crystals, or else water drops). Most of the 
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precipitating particles develop near the — 10° 
isotherm. At that level ice-crystals are formed 
which during the downward drop increase 
in size; to begin with mainly owing to diffu- 
sion and sublimation of the vapour which is 
highly super-saturated relative to ice. Later, 
when the particles have grown considerably, 
their further increase is mainly due to colli- 
sions with the surrounding cloud droplets. 
When the cloud is dense a rain drop may grow 
to such an extent that it becomes unstable and 
splits up into two or more drops. If the air 
below the boundary is relatively dry, a con- 
siderable evaporation of the rain drops will 
take place after their departure from the cloud. 

On factor of importance to the O18-abun- 
dance (a,) in the rain is the relationship be- 
tween the cloud droplets and the surrounding 
water vapour. On the surface of the droplet 
evaporation and condensation are constantly 
going on. These two processes mutually 
maintain equilibrium, if the size of the droplet 
remains constant. At all events, however, 
they will cause an exchange of water molecules 
between the vapour and the droplet. If this 
exchange is extremely rapid, the O1%-abun- 
dance in the droplet must correspond to that 
in the surrounding vapour. Whereas if the 
exchange is extremely slow the O18-abundance 
in the droplet will correspond with that of the 
vapour at the place where the droplet was 
formed, i.e., because of the oblique movement 
of the droplet, at a considerably lower level 
of the front. In reality the exchange goes on 
at a final, probably not very high, rate and 
consequently a precipitating particle while 
falling will accumulate cloud drops with an 
O!8-abundance corresponding to that of the 
vapour which is still at a lower level than the 
precipitating particle. Fig. 4 shows hypotheti- 
cal tracks of particles falling through a cloud. 
Owing to the fact that the tracks are diverted 
by the wind, especially while the precipitant 
particle is small, the lines drawn are not vertical. 
The curve d represent the course to be taken 
by certain precipitating particle in order to 
acquire the same O8-abundance as that of an 
actual particle, supposing the exchange between 
the cloud droplets and the vapour to be 
complete and immediate. Thus the rain col 
lected at C has followed the track of c, while 
its O18-abundance is representative of the 
H,O18-value of the cloud along the track of d. 
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Fig. 4. Imaginary tracks of precipitant particles falling 
through a cloud. 


In this connection it may be worth while 
to mention a method which, at any rate 
in principle, might permit us to determine 
what proportion of the precipitation is due to 
condensation and what proportion is due to 
collisions with cloud drops. At a certain 
locality in the cloud the O18-abundance of 
the cloud droplets, the vapour and the rain 
drops are designated a,, a,, and ay, respectively. 
Then a,, according to the above, will be bigger 
than 1.010: a, at 0°C. If the precipitation is 
formed exclusively through collisions, a, will 
be equal to a.. If a, is less than a,, the fraction, 
x, due to the condensation of vapour along 
the track, can be determined by the equation: 


dp =x 1,010+4,+ (1 —x)a,, 
or 
= a) 


A — 
1.010 * dy — 4, 


Hitherto, technical difficulties by. collecting 
separately water from cloud droplets, vapour 
and precipitant perticles have impeded an 
experimental determination of x. 


\ 


Since, inter alia, x in equation (2) is un- 
known, a quantitative determination ex- 
pressive of a, is not possible. It is possible to 
solve the problem only by assuming that the 
distance between the tracks of c and d (Fig. 4) 
is o. If so, the expression for a, will be a 
function, first, of the O18-abundance in the 
vapour before the upward glide starts (a,°); 
second, of the average temperature 


T= fr di 
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along the track of the first released precipi- 
tating particle (t being the temperature of 
the small element of water and water vapour, 
dq, taken up by the particle in an infinites- 
imal section of its track, and Q being the 
total volume absorbed); and, third, of the 
corresponding temperature, T, along the 
track of the observed precipitation. Leaving 
out of consideration the assumption of the 
distance between c and d being 0, an expression 
for ay may be deduced to give correct values 
for variations in the a, along the front, if 
1. in each of the tracks of precipitation in the 
cloud the air in rising from d to c is chilled 
in equal degree; 
2. no evaporation takes place from the drop 
during its fall from the cloud to the ground. 
1. presumably is satisfied in sufficiently 
close approximation, while 2. rarely or never 
is realized. Since the evaporation increases 
with the distance from cloud to ground, the 
observed difference will be less than the 
computed value. Hence the expression renders 
the maximum value for variations in dy. 
With the partial pressures of H,O and 
H,O1* termed p!® and p!*, Raoult’s Law of 
equilibrium between vapour and cloud drop- 


lets makes 
a 
1) ER. 18 
i 100 f° | 


and 
16 Clee »| 


100 /? 


(3) 


where pi and pi? represent the vapour pressures 
of pure H,O' and pure H,O1% respectively, 
and where the index, c, may be changed into 
p on the supposition that the part played by 
the process of collision is considerably greater 
than the part played by the condensation 


process. 
Hence the O18-abundance in vapour is 
18 
p 
a AT OO 
16 18 
Peake 
and after inserting (3) 
18 
Ay ; Po 
Aust: 16 ig 100% 
(100 - Ap) pe + app} 
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100 Pi 
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whence 


and 


Within the interval 258 < T < 293 equa- 
tion (1) may be approximated by the straight 
line 
Po _ 1.03139 - 7.80:10-5-T (6) 


i 
Po 


Next re-calculating (5) we get 


100 
Pi 16 18 
100 „Po\ ‚Po 
I 18 16 
dy Po Po 
16 
oh 100 „PR 
16 18 
100 (8-1 )ay Po 
Po 
16 16 
IOoo Po Po 
d = “a O. 8a, 
ee Dr 9999 ie 
whence 


dp = (1.03137 - 7-80-10 *T)a, (7) 
With A representing the total volume of 
H,O% in the M g of water vapour contained 
in one m? of saturated air, we have 
_ 20 Ma, 
18 100 5 


and 


11e 20 
1800 


(M: da, + a, - dM) (8) 


However, dA being equal also to the con- 
densed quantity of H,O18, we get 


20 


TEN ui 
Ann 1800 (9) 


1800 ? 
- (1.03137 - 7.80 10°5T)a,dM 
From (8) and (9) we get 
Mda, + a,dM = (1.03137 — 7.80 - 10-5T)a,dM, 
whence 
Mda, = (0.03137 — 7.80: 10 T )adM= 
=(c-c2T) a, dM, 
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and 
(to) 


Since M is equal to 100 times the density of 
saturated vapour or 


where E represents the pressure of saturated 
vapour in mm Hg at T° K, we get 


im—2s8(F AT) = > (REIT) (12) 


Teer N oie! 
By Clapeyron’s formula 
LdT 
dE= An (1 3) 


where v, and v, respectively represent the 
volume of 1 g water in the vapour phase and 
in the liquid phase, while L represents the heat 
of vaporisation. Since v,<v, and 
ER 

Eee 
(13) may be written 


18L 
Te 


Vi 


EdT 


which inserted in (12) gives 


yet & = 1) EAT 


T2\RT (14) 


and from (10) together with (11) and (14) 
we get 


da, 
dy 


(6-0 (zT Us) 


L however is dependent on T. In the inter- 
val 258 < T < 293 we may with good ap- 
proximation put 


L=745.7—0.550 T=(c-c,T) cal/g. 
Thereby (15) is transformed into 


da, 1 8C 18¢, I 
ER = (4-61) RT? e oh 1) „er 


HSE 18 I 
a (TR Ts \R (ciel G, at 


CO 
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By integration of this expression from a, to 
a,, respectively from T, to T, we get 


Al, 
ee EX) rs 
n- ar, 7 + Cy Ina t+ GT T,) 
or 
Gr Ci 
We. Tire lee TT Ge 


whence 
Ci 


a,=K-e T.T-G.eGT (16) 


Finally from (7) and (16) we get 
log a, =log (I+c, - QT) +log K-C, log T+ 


= (ar-à) 0.43429 (17) 
From this formula (17) a, can be determined as 
a function of T when once the constants have 
been calculated from the meteorological facts 
of the problem under consideration. 

Based upon the Nier value, 0.2039 % 
for the O!8-abundance in atmospheric oxygen 
(16), we obtain the values shown in the third 
column of Table 1; thus from column 2 


1.009 _ 2 
Rae CARO + for 


and for vapour from ocean water, 0.1990 - 
-0.991=0.1972 %. However, because of the 
variation with temperature of pif/pl, the latter 
quantity to some extent depends on this 
variable factor. In equation (3) substituting o 
for c and regarding a, as the Of-abundance 
in ocean water, equation (3) becomes valid 
for the partial pressures respectively of H,O™ 
and H,O?*; and substituting o for p in (4) this 
equation becomes valid for vapour from ocean 
water. Upon the insertion of a,=0.1990% and 
the elimination of less significant terms, we get 


Sey (18) 


0 
dy=0.1990— 5 % 
Po 
Table 2 shows a, at various vaporisation 
temperatures: 


0.2039 ocean water, 


Table 2 
u: DRS a 
O 273 0.1970, 
6 279 0.1970, 
13 286 0.19719 
19 292 0.1972, 
25 298 0.19735 


mean 0.1972 
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It will be seen that a, varies only slightly 
with the temperature of the sea, and if we 
assume the evaporation from the oceans, to be 
a,= 0.1972 %, no great error is committed. 

For the warm front earlier discussed we in- 
sert in equation (17) Ty =285°K, a, =0.1972 %, 
corresponding to the conditions existing in 
the warm air at the beginning of the upward 
motion taken together with the values in- 
troduced in the computation, viz: 


Gj =0,03137 
Enr DO 110, 2,grade 4 
(,= 745.7 cal/gram. 
¢,=0.550 cal/gram.grad. 
R=1.9864 cal/grad.mol. 


and further the values derived therefrom: 
ke Cw Gy and Ce Hand we ger 


log ay =log (1.03137 —7.80-107>T) + 1.44555 — 
— 0.76066 log T+ (4.92 : 10-4 T-225.89/T) - 
"0.43429 (19) 


The variation of a, with T is shown graphically 
ms Bios, 

The altitude of the warm front, whereby 
is understood the distance from the ground 
to the topmost level of the boundary from 
which precipitation reaches the ground, most 
likely was 3.5 to 4 km in the instance 
mentioned. At that altitude the temperature 
must have been in the vicinity of -7 to 
-11°C, assuming a saturation adiabatic rate 
extending upward (along the front) from an 
elevation of 200 m and a temperature of 
+12°C at that level. According to Fig. 5, 
a fall in temperature from +12°C to-7 or 
—11°C should cause a maximum change in 
a, from 0.1990 % to 0.1966 or 0.1960 %; 
hence, between 0.0024 or 0.0030 %; which is 
in satisfactory agreement with the measured 
value in Fig. 3: 0.0016, %, especially con- 
sidering the fact that in the calculation the 
evaporation from the rain drops in falling from 
the cloud to the ground was not included. 
Drawing a straight line in Fig. 3, through the 
levels of collections of rain, Nos 12 to 7, we 
get, by a line intersecting No. 6 vertically, the 
value for the O48 excess at No. 6, +0.0004 %, 
to some extent corrected for evaporation. Thus, 
the measured a, change is corrected to 0.0022 °/,, 
per cent which is in good agreement with the 
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Fig. 5. The calculated variation of the O!8-abundance of 
precipitation plotted against the temperature of con- 
densation. 


calculated value for the maximum change in 
Ay, VIZ: 0.0024 to 0.0030 %. 

However, the measured absolute values of 
a, differ considerably from the values cal- 
culated here. In the warm front here discussed 
the a, was found to vary between 0.1976 and 
0.1959, %; whereas the a, according to the 
calculations should have changed from 0.1990 
to about 0.1963 %. 

A reason for this is that the circulatory 
process of the water in nature is far from 
being so ideal as assumed in the above. Thus 
a, is not solely dependent on T, but also on 
the pre-history of the air. For example, the 
air may have absorbed some of its water 
vapour over solid ground, i.e., from fresh 
water of an O}8-abundance about 0.0018 % 
lower than in ocean water. In air of this 
description a, will be less than supposed in the 
case discussed in the above. However, if 
maritime air only is considered deviations of 
this nature should be of slight importance. 

More important is the fact that air rarely 
is chilled so much that the absolute humidity 
practically becomes to. This means that the 
general circulatory process, which consists of 
the following three processes: 
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1. Saturated air at Ty (with a? as stated in 

Table 2) will liberate rain by cooling to T°; 
2. The air is again heated to T,°; 

à: The air is saturated with water vapour at 

0 >» 

will alter the a, considerably, bacause the a, 
remaining in the air after process 1. will be 
less than aj. When the air in process 3. has 
absorbed a certain quantity of water vapour of 
O!8_abundance aÿ, the a, for the total quantity 
of vapour must also be less than a}. 

Such a fall in a, naturally cannot proceed 
indefinitely. Since the air continually passes 
through circulatory processes, an equilibrium 
has been established where the final a, on an 
average is just so much below the values 
stated in Table 2 that the decrease in a, during 
process 1. is balanced by the increase in a, 
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during process 3. The final mean value of a, is 
dependent on the climatic and geographic 
conditions. 

A third reason for the disagreement in the 
absolute values is, no doubt, the fact that the 
warm air in question had given off moisture 
before starting the upward movement at 
Copenhagen. 

However, the significant part of this in- 
vestigation is that these considerations have 
no effect on the agreement between the com- 
puted and the measured, corrected changes in 
dy. Actually, the insertion in (17) of a value 
for a, differing 0.0014 %, from the value 
earlier used makes no essential difference in the 
calculated a, changes. 

Further investigations into these problems 
will be discussed in a later paper. 
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A Note on the Chemical Composition of Rainwater 


By W. A. MORDY, Honolulu, T. H. 


(Manuscript received June 2I, 1953) 


Abstract 


A limited number of observations of nitrogen, sodium and chloride in rainwater made 
in Hawaii are presented. Because of frequent and high rainfall, calculations of the annual 
amounts of nitrogen received in various Island localities are rather large. In one storm 
nitrogen was brought down by rainfall at a nearly constant rate although rainfall intensity 
varied considerably. 

In general, the concentrations of these elements in Hawaiian rains agree with observa- 
tions in other oceanic areas. 


In view of the current interest in atmospheric (ammonia and nitrate) was determined by 
nitrogen, Dr C. G. Rossby recently asked distillation in the presence of alkali and 
me to submit to the readers of Tellus some Devarda’s alloy, reducing the nitrate to 
observations of nitrogen and other chemical ammonia. Ammonia in the distillate was 


constituents in rainwater in Hawaii. 
During March and April 1953, several 

samples of rainwater were collected and 

analyzed for nitrogen content. Total nitrogen 


nesslerized and the color was measured on an 
Evelyn photo-electric colorimeter. 

The results of these determinations are 
shown in Table r. It appears that these figures 


Table 1 


Elapsed] Rate of | Rainfall 


Volume} Total 


Total N ? 


Ean : , 
Date re Time |Collection|Intensity|. |, Rain | Ki -no,.n IN/min) Na® CIS 
Sample |for Day 

min | cm?/min | mm/h cm? cm mg/l g/min | mg/l | mg/l 

3/11,53 | 9: 34—0: 37 | 3 220 1622 660 .16 35.2 

9:38 3 4 

9:43 1/2 4—3/4 97 5:9 460 1.19 Ses 10.6 

3/27/53 14: 24—14:27| 3 188 11.3 566 47 -08 15.1 
4/22/53 |14: 45—15: 10] 25 38 1.8 740 -5I 15-1 2.8 4.8 
15: 12—15: 43 31 23 1.4 715 .39 9.0 1.4 3.0 
15-44 15555 or 65 3:9 720 .10 6,6 .6 .6 

16: 08—16: 25| 17 35 2.1 591 -24 8.4 
16: 28—17: 30] 62 THE 7 707 35 .66 7.5 2.0 2.0 

| Av =.28 


Sodium Determined by Flame Photometer. 
4 Chloride Determined by Titration. 


Hawaiin Standard Time=150th Meridian Time. 
Nitrogen Determined by Distillation and Colorimetry. 


Published with the approval of the Director as Technical Paper No. 218 of the Pineapple Institute of Hawaii. 
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Fig. 1. Relationship of rainfall intensity to rainwater 
nitrogen content for 8 samples of rainwater collected in 
Honolulu for March 11 to April 22, 1953. 


are in good agreement with the figures found 
elsewhere in ocean climates, Erıksson (1952). 
They are somewhat lower than the figures 
that Äncström (1952) indicated should be 
found for rainfall from tropical air masses, 
and the relationship which he suggests con- 
cerning the decrease in the rate of collection of 
nitrogen as the rainfall endures, is not borne 
out in the one test of it here on April 22 
where five consecutive samples were obtained. 

An inverse relationship between nitrogen 
content and rainfall intensity appears to exist 


21°30’ 


LOCATION OF 
OBSERVATORY 


158°00° 


Fig. 2. Location of meteorological observatory Experi- 

ment Station of Hawaiian Sugar Planters’ Association, 

where rainwater samples were collected. Rainwater 

was collected orographic rainfall 4 miles (6 km) leeward 
from the crest of the Koolau Mountains. 
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in these data. Only one observation, the first 
on March 11, 1953, lies away from the line 
sketched in on the scatter diagram, Fig. 1. 
If the nitrogen is absorbed by the rainwater 
in a washing action of the rain an inverse 
relationship might be expected. Such a rela- 
tionship also might be expected to show 
better in these data from Hawaii due to the 
nature of the rainfall. These samples were of 
orographic rainfall collected at the meteor- 
ological observatory on top of the Hawaiian 
Sugar Planters’ Association Experiment Station 
4 miles (6 km) leeward from the crest of the 
Koolau Mountains, Fig. 2. Because the rain- 
fall was orographic it can be considered 
essentially stationary with new air pouring 
into the rainy zone with the speed of the wind. 
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Fig. 3. Rate of accumulation of nitrogen (curve 1) with 

the accumulation of rainwater (curve 2), on April 22, 

1953. The rate of accumulation of nitrogen per unit 

time appears to be quite constant, while the rate of 

accumulation of rainwater as a result of changes in rain- 
fall intensity is quite different. 


If the air contains only a fixed amount of 
nitrogen which can be brought down by 
the rain, then the rainwater from the higher 
intensity rains will contain less nitrogen per 
unit volume of water than the lower inten- 
sity rain. Such reasoning would also require, 
of course, that the washing efficiency of high 
and low intensity rains be nearly the same. 
The amount of nitrogen accumulated in the 
rainwater per unit time has been examined 
with this reasoning in mind. In Fig. 3 are 
graphs of the accumulated nitrogen and the 
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accumulated rainfall after a given elapsed time 
from the beginning of the rainfall. It is in- 
teresting to see that the rate of nittogen accu- 
mulation (slope of curve 1) seems fairly steady 
for a two hour period, while the rate of rain- 
fall (curve 2) accumulation varies significantly. 
This is interpreted by the author, somewhat 
intuitively, to mean that the nitrogen exists 
as a limited supply in a given volume of the 
atmosphere and is almost completely brought 
down from the air through which the rain falls. 

A further point of interest lies in the rela- 
tionship between the nitrogen and sodium or 
chloride determinations. Sodium and chloride 
determinations have been made on a number 
of rainfall samples in cloud physics studies 
during the past ycar. Variations in the latter 
elements have been accounted for by postu- 
lating the role the sea salt particles in the 
atmosphere play in the formation of raindrops, 
(Wooncock 1951, BLANCHARD 1952). It is 
interesting therefore to compare the amounts 
of nitrogen with the sodium and chloride 
determinations to determine: 


1. If the nitrogen in the rainwater could be 
directly from the sea as it has been reasoned 
the salt is; 


VErEMZAEMENERZAZSZIEOENED 
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HONOLULU W.B.O. 
1909 - 1947 


Fig. 4. The amount of rain at the Honolulu Weather 

Bureau office for the years 1909 to 1949 inclusive, which 

resulted from rains of 0.2 in. (0.51 cm) or less, 0.4 in. 

(1.02 cm) or less, I in. (2.54 cm) or less per day, and 

the total rain received for the 39 year period by months. 

Rains of less than an inch per day account for more than 
half the total rainfall. 
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2. If the nitrogen and salt exist as particles in 
the air which might be washed down 
rather than the salt being the nucleus for 
raindrops. 


If the nitrogen and sea salt get into the rain- 
water by the same route from the sea, then 
the ratio of salt to nitrogen should be relatively 
the same in rain and sea water or at least 
the amounts of nitrogen and chlorides should 
be strongly correlated. This is not the case 
since the nitrogen to chloride ratios in the 
rainwater exceeds roughly 1,000 times the 
nitrogen to chloride ratio in sea water 
(SVERDRUP 1946). 


(0) 


JFMAMJIJSAS OND 


NUUANU DAM 4 
1905 - 1947 


Fig. 5. Rainfall at Nuuanu Dam No. 4 (average rainfall 

150 in.) which has resulted from rains of various amounts 

per day. Rains smaller than ı in. (2.54 cm) per day 
count for approximately half of the total rainfall. 
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Secondly, if the nitrogen and salt were 
part of the same particles, or separate particles 
subject to the washing action of the rain 
the ratio of nitrogen to salt in the different 
rainfall samples taken within relatively short 
periods of time should be fairly constant. In 
the four April 22 cases in Table 2 where this 
comparison can be made, this is not the case. 


The amount of nitrogen which arrives in 
the soil as a result of rainfall is a main point 
of interest. In view of the apparent inverse 
relationship between rainfall intensity and 
nitrogen content, some rough calculations 
can be performed which are likely to be con- 
servative, since the rainfall was in each case 
heavier than the median rainfall and simi- 
larly of greater intensity. 


To perform these calculations it is helpful 
to look first at monthly frequency distributions 
of two rainfall stations, the Honolulu U.S. 
Weather Bureau office with 30 inches (76 cm) 
average annual rainfall, and Nuuanu Dam No. 4 
with 150 inches (381 cm) average annual 
rainfall, on the Island of Oahu. These distribu- 
tions are different from the usual represen- 
tation of such data in that they show the 
amount of the total rainwater which is re- 
ceived by each category of rainfall (Figs 5, 6). 
It can be seen that daily rains which are one 
inch or less comprise over half the total 
rainfall for the thirty inch station, and rains 
of one inch or less similarly comprise over 
half the total rainfall at the 150 inch station. 


In another study the author has examined 
the rainfall from a number of such stations 
and found that in general the increase in 
rainfall in the wetter areas of the island is due 
to increased frequency of smaller rains parti- 
cularly, although of course the frequency of 
all intensities of rain is greater. 


These statements are intended to help show 
that if we calculate nitrogen arriving on the 
surface of the land by using the figures from 
the graph in Fig. 1 or by using the average 
nitrogen figure of .28 from Table 1, we are 
likely to arrive at a result somewhat less than 
the actual, the average nitrogen/unit vol. figure 
producing the lower value of the two for 
nitrogen per unit area. 


Using the average value of nitrogen/unit 
vol. from Table 1, the results obtained were: 
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Table 2 
Fe : 
verage Rainfall Calculated 
Kur a Kg N/Ha/year 

30 76.2 2 
100 254 7 
200 508 aA 
400 I O16 28 
600 I 524 43 


Later, more meaningful calculations can 
be made when more data have been assembled 
showing the relationship of nitrogen to rain- 
drop size, rainfall intensity, and daily rainfall. 

However, in lieu of these calculations it is 
interesting to examine some unpublished ob- 
servations of monthly amounts of nitrogen 
in rain water made here at the HSPA Experi- 
ment Station (38 inches annual rainfall) and 
the HSPA Manoa Arboretum (160 inches 
annual rainfall), Fig. 2, by L.R. Smith in 
1935—36. These data are shown in Table 3. 

The methods which Smith used in his 
determinations are subject to some criticism. 
1) The use of zinc in his alkaline distillation 
was likely to reduce some of the nitrate 
present to ammonia. The ammonia values 
therefore may be too high. 2) No provision 
was made to remove chlorides. In the pres- 
ence of chlorides and upon the addition of 
concentrated sulphuric acid in the phenol- 
disulfonic acid reagent, nitrate may have been 
lost by reduction to gaseous nitric oxide. 
Nitrate values may then be too low. The 
magnitude of the gain in ammonia and the 
loss in nitrate is of course speculative. Assuming 
that these errors are compensatory, then the 
total nitrogen values may not be too greatly 
in error. It can be seen that the estimates of 
nitrogen from the average determinations is 
low for the Makiki case and high for the 
Manoa station. 

As a final point of discussion, one factor 
stands out in Eriksson’s table of values of 
nitrogen determinations for various parts of 
the world. This is the near uniformity of the 
low values which have been obtained in 
isolated oceanic regions. There appears to be 
little difference between areas of frequent 
storms where electrical discharges would not 
be uncommon and areas such as Hawaii where a 
thunderstorm or lightning discharge is very 
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Table 3 
u m [EEE —— 
Manoa Total N Makiki Total N 
1936—1937 Arboretum Rainfall for Experiment station | Rainfall en 
th 
NH,N | NON end ANELNEENO,S Oe 
Lk ce eee | re ee | ee ne ee le u de 
Kg/Ha | Kg/Ha Cm Kg/Ha | Kg/Ha Kg/Ha Cm Kg/Ha 
September er: .I4 27 53.59 AI .28 .05 7.01 33 
October ee 53 25 52.83 -79 23 .15 29.26 37 
November ee 97 — 19.33 .97 .15 .06 6.99 .20 
December .19 .I4 52.02 .33 89 -09 17.73 .98 
Januar. ee .55 .36 70.28 .92 21 21 21.36 .43 
HÉDITAL EEE EEE == — 41.63 — .80 .23 18.85 1.03 
Total 5 months 3.42 Total 6 months 3.34 


infrequent. Although the proportion of ni- 
trogen in the rainwater is fairly constant from 
one ocean area to another, the amount and 
frequency of rainfall, of course, is very differ- 
ent. Hence the distribution of nitrogen 
brought down by the rainfall over different 
parts of the world should be very uneven 


indeed. 


A continuing program of nitrogen obser- 


# 


vations is planned to attempt to relate nitrogen 
with drop size, rainfall intensity, character, 
and the general weather synopsis. 
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On the Production of Mean Zonal Currents in the 


Atmosphere by Large Disturbances 


By HSIAO-LAN KUO, Massachusetts Institute of Technology! 


Manuscript received July 22, 1952. Revised manuscript received June 2, 1953 


Abstract 


The problem of the production of mean zonal currents by large-scale horizontal disturbances 
in the atmosphere is attacked by computation of time-tendencies from the vorticity equation. 
It is shown that the earth’s rotation, acting through the variable Coriolis parameter, generally 
produces a mean northward transport of westerly momentum, with the maximum transport 
occurring south of the middle of the zonal belt occupied by the disturbances, and therefore 
creating mean westerly and easterly currents in the central and outer parts respectively. Its 
effect on the disturbance may either be a damping or an amplification, depending on the asym- 
metry of the mean flow profile. 

Six different mean flow profiles have been investigated to determine the nature of the in- 
teraction between a disturbance and the mean flow. For three profiles with inflection points, 
the interaction depends upon the wave length of the disturbance. If the wave length is longer 
than a critical value, an amplifying effect on the disturbance is produced accompanied by a 
divergence of westerly momentum in the region of maximum mean flow, and if shorter, a 
damping effect with a convergence of westerly momentum. On the other hand, the inter- 
action produces only a damping effect for three mean flow profiles without inflection point. 

Since the zonal currents in the atmosphere are seldom very strong and not far from symmet- 
rical, the total effect of the earth’s rotation and of the mean flow leads to horizontal damping, 
in which the contrast in the mean flow at different latitudes is increased and a mean northward 
transport of westerly momentum is produced. This is in qualitative agreement with observa- 
tions, indicating that these factors are of real importance in the development and maintenance 


of the mean zonal flow. 


I. Introduction 


In previous studies, the author has discussed 
the mechanism of the mean meridional 
transport of westerly momentum by large- 
scale atmospheric disturbances (KUO, 1951 a, 
1951 b), based on the properties of solutions 
of the linearized vorticity equation for horizon- 
tal, inviscid, barotropic, non-divergent flow. 
This will again be our subject. Although 
vertical motion and baroclinicity in the 
atmosphere are of paramount importance in 
some other aspects of the meteorological 


problem, it is believed that the omission of 
such factors does not significantly distort the 
average interactions between the purely hori- 
zontal components of the large-scale dis- 
turbances on the one hand and the mean flow 
and the earth’s rotation on the other, as can 
be demonstrated by integrating the equations 
along the vertical with respect to pressure, or 
by comparing the simple vorticity equation 
with the general differential equation for 
quasi-geostrophic flow. In previous studies, 


1 The research reported in this document has been made possible through support and sponsorship extended by the 
Geophysics Research Division of the AF Cambridge Research Labs, AMC, U.S. Air Force, under Contract No. 
AF 19—122—153. It is published for technical information only and does not represent recommendations or con- 


clusions of the sponsoring agency. 
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the classical eigenvalue (perturbation) method 
was used, where the disturbances are assumed 
to be harmonic both in the coordinate x and 
in time f, so that each fundamental harmon- 
ic is represented by a function of the type 
f(y) exp (ik(x - ct)). The problem is then to find 
the permissible values of c (the eigenvalues) 
and the function f(y) (the eigenfunction) for 
different wave numbers k and a given mean 
flow profile U(y). Using this method it was 
possible to find some general properties of the 
disturbances related to the nature of the mean 
flow without actually solving the eigenvalue 
problem for each k and U. It was shown that 
if the absolute vorticity of the mean zonal 
flow has extreme values in the zone considered, 
and increases with latitude where the mean 
zonal velocity U is large and decreases with 
latitude where the mean zonal velocity U is 
small, the linear vorticity equation permits 
solutions with non-vanishing real exponential 
time factor (c is complex). In this case, dis- 
turbances of wave length shorter than a 
certain value will decrease exponentially with 
time and are therefore “damped”!, and those 
longer than the critical wave length will 
increase exponentially, i.e., “amplify”. The 
effect of the damped shorter disturbances is, 
in the mean, to transport zonal westerly mo- 
mentum from regions of weaker to regions of 
stronger mean zonal flow, thereby accentuating 
the contrast in the existing zonal flow, leading 
to the development of a sharper jet type mean 
flow profile. The effect of the amplifying 
disturbances is in the reverse direction, flat- 
tening the existing mean flow profile. 

If, on the other hand, the absolute vorticity 
of the mean flow does not have an extreme 
value, then no such eigenvalues exist for the 
inviscid vorticity equation. The failure of the 
eigenvalue method for this case is due to the 
assumption of a particular harmonic time 
dependence, which is a very restrictive assump- 
tion. If the problem is one relating merely to 
the stability or instability of steady mean 
flows which can be established in a real fluid, 


1 The term ‘damped’? and ‘“amplified” are used, 
as in the previous investigations on horizontal non- 
divergent motions, merely to describe the interactions 
between the horizontal disturbances and the mean flow 
or the earth’s rotation; they need not be identified 
with actual changes of any disturbance in the atmos- 
phere, since these disturbances may be produced by 
the factors that have been neglected. 
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it is natural to take viscosity into consideration, 
whereupon the difficulty is removed and 
definite results can again be obtained by the 
eigenvalue method. Thus in the case mentioned 
above, most of the small disturbances will be 
damped in a flow of very high Reynolds 
number. 

However, the inclusion of viscosity seems 
somewhat artificial when treating the large- 
scale motions of the atmosphere, since the 
Reynolds number is very large and the motion 
of a nonviscous fluid should also remain 
determinate. Furthermore, although a mean 
zonal current almost always exists, large 
disturbances are also invariably present; a 
purely zonal steady state flow has rarely 
been observed. Thus, in considering the 
production of zonal motions by large dis- 
turbances, we are confronted by an inherently 
nonlinear process, therefore, the question of 
the validity of the perturbation method arises 
in addition to the difficulties in its application 
mentioned above. Furthermore, in discussing 
this process, we are less concerned with the 
stability or instability of the existing mean 
flow, which are difficult to define when the 


_ disturbances are finite, than with the effects 


of the interactions between the disturbances 
and the mean flow or the earth’s rotation. In 
view of these difficulties confronting the 
perturbation method, it becomes most de- 
sirable to approach this problem by another 
method. 

The method to be used in the present in- 
vestigation is to assume the existance of a 
mean flow and of certain disturbances and to 
compute the subsequent changes from the 
nonlinear vorticity equation. This method is 
similar to the one used by TAYLOR and GREEN 
(1937) in their study of the production of 
small eddies. The same method has also been 
used by PLATZMANN (1952) in a recent paper. 
In principle, the general solution corresponding 
to a given initial condition can be found by 
repeating the process and expanding the 
quantities in a Taylor’s series in t. However, 
there is the question of convergence of the 
series, and it is found that the labor involved 
in obtaining the higher derivatives becomes 
prohibitive even for very simple initial dis- 
turbances. We must therefore be satisfied with 
only one or two terms in the expansion. But 
since we are not particularly interested in the 
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instantaneous rate of change, but in the end 
effect of the interaction, the question arises as 
to whether the single term represents only 
a transient effect. 

Still another difficulty arises in the general 
interpretation of the results obtained from a 
specific model. Since the behavior of a dis- 
turbance is directly and continuously reflected 
in the form of the disturbance, the result 
obtained for any one disturbance is valid only 
for that disturbance. This difficulty is absent 
in some measure in the eigenvalue solutions, 
as both the time dependence of the solutions 
and the forms of the fundamental modes of 
the disturbances are determined by the eigen- 
values and the eigenfunctions, so that all 
effects can be attributed to the character of 
the mean flow and to the scale of the disturb- 
ance. This is precisely what we would like to 
do in our present study. We must therefore 
attempt to eliminate those effects that are the 
result of the particular form of the initial 
disturbance. 

One way of eliminating the effects that 
depend mainly on the initial form of the 
disturbance is to consider arbitrary or random 
disturbances, which may be expressed as a 
Fourier series in x and y. So far as the effect 
of earth’s rotation is concerned, we are able 
to extend the results to the general initial 
disturbance and demonstrate that whenever 
the disturbance is not represented by the 
exactly neutral solution of the vorticity 
equation, it always produces a mean north- 
ward transport of zonal momentum in almost 
all parts of the belt (northern hemisphere) 
except near the northern boundary, where a 
reverse transport may occur. On the other 
hand, no such generalization is possible 
regarding the interaction between the dis- 
turbance and the mean flow, which actually 
depends both on the form of the disturbance 
and on the type of mean flow, and therefore 
may change direction. We therefore choose 
only one simple initial wave-disturbance for 
the study of this effect, which initially does not 
produce a mean momentum transport. How- 
ever, we shall take it to be different from the 
exactly neutral disturbances if such disturbances 
exist, either by taking a different amplitude 
function or by taking an arbitrary wave length 
so that some change will be produced later. 
Since all these initial changes are produced by 
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the linear terms of the vorticity equation, the 
results obtained should be comparable with 
those obtained from the linear vorticity 
equation, at least qualitatively. It should be 
emphatically remarked that as no energy 
source is included in the simple vorticity 
equation, we merely assume the existence of 
an arbitrary disturbance rather than to identify 
it as given by the neutral solutions of this 
equation. 


2. The production of mean zonal currents 
on a rotating earth. 


Let us assume that in the atmosphere 
there exist disturbances at the initial moment, 
presumably produced through some external 
processes such as differential heating or the 
release of potential energy but which are now 
left to be controlled by the simple vorticity 
equation for nondivergent, barotropic, hori- 
zontal flow. 

For motions over a spherical earth this con- 
trolling equation may be written 


a DER OL EN ee (a9) 


at Roos p aa Rad 
where A and are the longitude and latitude, 
f=2m sin y is the Coriolis parameter, R the 
radius of the earth, u and v are the linear 
west—east and south—north velocity com- 
ponents and ¢ is the vertical component of 
the vorticity of the motion, respectively. For 
this motion, a stream function y can be 


introduced, defined by 
u = RA cos Park, | 

(2.2) 

v cos p = Rp cosp = = 


where A and ® are the velocity components 
in angular measure. In terms of this stream 
function and the new independent variable 
n=sin y, the vorticity equation (2.1) takes 
the symmetrical form 


Oy _ dy 0V?p ap aV*y +f) 
dt on 2 À on 


where the symbol V* denotes the two- 
dimensional Laplacian operator in spherical 
coordinates and Y2p is proportional to the 
relative vorticity, 


(2.3) 
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ale. 
X“ PT cos? y OR er \ 
(2.4) 


To avoid complications, we at first study 
the simple initial disturbance with only one 
wave-component and north—south trough- 
and ridgelines, such that no mean momentum 
transfer is produced at the initial moment, 
and then extend the treatment to the more 
general disturbance with several wave-com- 
ponents and with tilted troughs and ridges, 
1.e., to the most general disturbance represented 


by 


Me 


y= LAn(n, dei (2.5) 
0 


m 


where H is generally complex. 

Since the amplitude of the waves varies, we 
shall now consider the proper boundary con- 
ditions that should be applied to the motion. 

From symmetry considerations, it seems 
natural to assume that there is no momentum 
flux across the equator, which can be satisfied 
by requiring either u or v to be zero. We shall 
assume v is to vanish, which requires the 
vanishing of the amplitude functions identically 
along the equator. We note that the parts of 
the north—south velocity from the wave- 
components with m larger than 1 must be 
zero at the pole, since the corresponding 
components of the stream function are sym- 
metric about this point. For the component 
m=1, the corresponding stream function is 
not symmetric about the pole, and there is 
one stream line which cuts through this point, 
therefore v is finite. This represents a single 
vortex centered off the pole. Therefore, we 
must require that the Fourier-component of 
v with m=1 to be finite and all the other 
components with m = 2 to vanish at the pole, 
if this point is to be included in the region of 
consideration. This requires the functions 


co | 
(1-7?) *H,(n, t) to be finite for m=1 and to 
be zero for all m > 2. In this section, we shall 
concentrate on the effect of the earth’s rotation, 
and therefore assume that there is no basic 
current. 


a. A simple wave disturbance. We first discuss 
the motion given by 


Vo = aF,(n) sin mA (2.6) 
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where F,,(y) satisfies the boundary conditions. 
Substituting in (2.3), we get 


2 
yt = - 2wmaF(n) cos mA+ = w(n) sin 2mA 


(27) 


where 
_dF 9 dF) _@ P =] 
vd =, gL | Fe |; 
& 2m? 
= 7 


à À 
Therefore, the solution — must be of the 


ot 
following form 


2 
7 2wmay(n) cos mA + = Z(n) sin 2mA (2.8) 


where y and Z satisfy the differential equations 


d dy m? 
LE Zr ey y F (2.9 a) 


d dZ 4m? 
Le ne | 3 Te az 


and the boundary conditions 


a, -3 
(Len?) 010) 


2-0} (59 
for n„=oandn=1ı | 


Substituting the solution (2.8) in the expres- 


sion for a we find 
dt 
(uv) _ (2 Py op Py 
at On HOA OK Aton) se 
DRM DT pl? in 
2) dn dn 


where the bars denote averages around the 
2 : : O(uv 

latitude circles. It is seen that 1 depends 

only on y and F, that is, depends upon the 

interaction between the disturbance and earth’s 


df 


rotation, as represented by the term v Rie in 

} pP 
eq. (2.1). Since Z does not contribute to 
O(uv) 


ap ake shall find y only. In order to deal 
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with an arbitrary F, it is convenient to use 
Green’s method, writing y in the form 


y = Jn, £) F(é)dé (2.11) 


where G is the Green’s function, satisfying the 
homogeneous part of (2.9 a) and the boundary 
conditions and has a jump of the first derivative 
for n=£. It is given by 


m m m 


_ (1-8 (x40)? (um? 
en) ea a 


for n <& 


my m m 


I-n\2| (1 +8&\2 BSG 
= el ae (=) (3) ; 
for n > & 

(2.12) 


Because of the particular form of the La- 
placian operator (2.4), it is sometimes con- 
venient to express F,,(7) in terms of spherical 
harmonics. It is evident that any sum or 


product of the functions 
fal 


F(n) =apr (m) =a, y(t — 42) ? (2-73) 


satisfies our boundary conditions. For sim- 
plicity, we shall simply take (2.13) as the 
amplitude function of the initial disturbance 
(2.6), and leave the more general case for 
treatment in section (c). Substituting in (2.11) 
and then in (2.10), we fin 


See n—3 | 
R2a?m 26 n°) 2 aye Sf on nn?) ö 


n 


Alu) 
dt > 


en 2 T = é\2 ah 
Efron 
n (2.14) 


It may be mentioned that disturbances with 
their relative vorticity proportional to the 
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stream function, i.e., A’y= -Ky, are exactly 
neutral. For these disturbances, w and therefore 
also Z are identically zero and y is a constant 
multiple of the initial amplitude function 


which is given by p„(n). These disturbances 
merely rotate around the pole without change 
of shape, and no mean transport of momentum 
is produced by them at any time.! Thus the 
case with m=n-ı belongs to this group of 
exactly neutral disturbances. We note that if 
we take a sum or product of the functions of 
eq. (2.13) as the amplitude function, no 
disturbance could be exactly neutral. 

For all values of mÆn-—1, (2.14) gives a 


Cay 
OSITIVE —~— al €Ver omt. Increases from 
t t y point. It 13 


dt 

zero at the equator to a positive maximum 
somewhere south of @=45° and decreases 
to zero toward the pole. The value of this 
function for m=4 and n=3, when divided 
by the constant factor 4v?U/D to make it non- 
dimensional, is represented by the dotted 
curve in Fig. 1. 

Similar results may be obtained for dis- 
turbances with amplitude functions different 
from that given by (2.13), and it is found that 
if the maximum is stronger, then a slight 
southward transport in the region near the 
pole will occur, which also appears to be 
verified observationally. 


b. A single vortex. As has been mentioned 
before, a requirement of v=o at the pole is not 
justified for the motion of a single vortex 
centered off the pole, or a wave motion of 
m=t. Since this motion is also of some interest, 
we shall analyze its effect on the mean mo- 
mentum transfer. As the motion represented 
by p2(7) sin À is exactly neutral, we shall take 
a different amplitude function. The simplest 
initial disturbance which satisfies the conditions 
that v=o at the equator and gives a finite v at 
the pole and which is not exactly neutral is 

1: 


Wy =an?(1 -n?)” sin À (233) 


Substituting in (2.10) we find 


RB ones na) 
ot 60 toe a 21 i 


1 For discussions of these exactly neutral 
disturbances see Neamtan (1946) and other papers. 


0.4 0.5 0.6 07 


Fig. I. 

corresponds to two sets of disturbances for profile (b), 

shifted a disturbance D/3. Full curve is that given by 
(3.3a) and the dotted curve is given by (2.14). 


Theoretical distribution of uv, dashed curve 


= (LONS + isn + nf 1) + (2-29-38) 


mr? on) (2.10 


which is also positive for all values of n, 
and has its maximum around 45°. 


c. Disturbances limited between two parallels of 
latitude. Although the preceding results are 
obtained under the assumption that v vanishes 
identically at the equator and at the pole but 
not at any other latitudes, similar results will 
be obtained whenever the amplitude of the 
disturbance varies with latitude, which may 
be considered as determined by the supply of 
energy to the disturbances in the atmosphere. 
Thus we shall investigate the disturbances that 
are limited to some other belt, and demand 
the vanishing of v along two latitude circles 
@1 and 9,. To satisfy these boundary condi- 
tions for unspecified positions of these parallels, 
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it is more convenient to use the new coordinate 
system 


(2.17) 


Then, after multiplying by cos? y, the vor- 
ticity equation (2.3) becomes 


CU REN ee 
DIN OR oy" 
i Ip 7p _ op IV *y 
dy 0% dx, “Oy 


x=A, y=log (tan y +sec y) 


0 
— 2 sech?y a + 


ee. 


It is seen that the Laplacian operator on the left 
side has the same form as in Cartesian coordi- 
nates. The symmetry of the last two terms on 
the right has been preserved by using the 
operator Y? given by (2.4). We note that the 

earth rotation term takes the form 

dw dy 

2 Auer NEE 

2@ cos 9-20 sec h u. 
in this equation. The origin of y can be 
changed to latitude y,, by changing sech?y 
into sech?(y+y,). Then the boundary condi- 
tions will be v=o for y=oand y=y,, corres- 
ponding to y, and 9. We also assume the 
initial disturbance to be a single wave of the 


type (2.6). Then a takes the same form (2.8), 


dt 
while now y satisfies the differential equation 
d2 
aye My = FO) sec h(y+y)= —R(y) (2.15) 


Both y and F should vanish at y=o and 
y=7Y2. Here we assume y, to be finite, i.e., 
Pa < 90°. Since Z does not contribute to the 
change of the mean momentum transfer, we 
shall not discuss it. 

We shall now specify the amplitude varia- 
tions by assuming that F has its maximum 
near the latitude y=y,/2. Thus we put 

N 


Fy) = A(y) 25,0" 


where O=sin = and A(y) is a monotone 


2 

function of y. Each term of F has its maximum 
near y=y,/2 and satisfies the boundary condi- 
tions. As the term in the vorticity equation is 
linear in y, the effect of the disturbance with 
this amplitude function is similar to that given 
by a single term of it. Therefore we shall 
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simply put F to be proportional to @". To 
simplify the integration, we take A(y) =cs?(y + 
+71), so that R(y) =O". Then y is given by 


n 


y= Ea" Go f à - SM Tan) | 


TT 


A 


2.20) 
cs 


where pee and cs stands for the hyperbolic 
7 


cosine, and the coefficients are given by 


3 Ri Ye a 2 . 2 ORO ac 
‘ fat UE (no) ru: 


I n!} 
(n— 21) +? (n— 21)! 


Gel 2-0 


When n is even, c, is given by the upper 
expression by putting 2/=n; when n is odd, 
€9=0. 


Substituting y in the expression for a we 
find 
om 
cost ee = R’m?a®w(M,+M;) (221) 
where 


M, = Pen cos! > mn er 
lv: Ya — UT 


CS Fr 


2 


NIT ny( um I 
. cos cs— —csm(y —— + 
[ Ya Ya ( 2 Y =) 


= mOSnm(y = a) |} 


M, = 2 sin gO"y 


(2.22) 


Similar results are obtained if we choose A=1; 
then M, will be larger. It can be shown that 
M, is positive for y < y,/2 and negative for 
y > 2/2, therefore, has the effect of producing 
a westerly current in the central part of the 
belt and easterly currents in the outer parts. 
On the other hand M, is positive for all y and 
has its maximum near y=y,/2, therefore 
produces a westerly current north of this 
latitude and an easterly current to the south. 
It can also be seen that for smaller values of 
-n, M, is much larger than M,. For example, 
M, disappears when n=1. As n increases, or 
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as the maximum of the amplitude becomes 
more pronounced, M, becomes larger than 
M,. The sum of these two terms gives a 
northward transport of momentum for almost 
all y except near the northern boundary, 
where a reverse transport may result if M, is 
larger than M,;. The value of the transfer 
corresponding to n=ı (M,=o), with the 
factors sin @ and cos @ omitted, is plotted in 
Fig. 1, as the fully drawn curve. The case for 
n > 1 is something like those given by the 
curves b, c, or d in Fig. 9. 

These examples show that whenever the 
disturbance is not exactly neutral, the interac- 
tion between it and the earth’s rotation pro- 
duces a northward transport of momentum, 
except at very high latitudes where a reverse 
transport may occur. The maximum transport 
always occurs somewhere south of the mid- 
point of the belt in northern hemisphere. 
Figuratively, it may be said that this interaction 
tends to tilt the trough- and ridge-lines in a 
SE—NE direction in almost the entire region 
except near the northern boundary where 
they may take a SE—NW direction. 


d. Disturbances with more wave-com- 
ponents. We now discuss the disturbance re- 


presented by 
vo=ajFıln) sin A+ ayFn(y) sin mA (2.23) 


where / and m are different integers. Substi- 
tuting in (2.3) it can be seen that + must 
take the form 

op 

5, zolaryı (n) cos A+2@maniym (n) cos mA + 


+ Z,(y) sin 21 + Z,,(y) sin 2mA + Zy(n) sin 
(m+ A+ Z,(n) sin(m—l)A (2.24) 


The Z-terms are new wave-components 
while y’s satisfy eq. (2.9 a) or (2.19) with the 
proper F-function on the right. Substituting in 


, Ouv 
the expression for ae we get 


Ae = RP ato vig" = i) + 


2.25 
Al R2m2a2,,0 ( y a I =) 
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It is seen that each wave-component produces 
its own momentum transfer because of the 
earth’s rotation, if it is not exactly neutral. 
Thus the results obtained before can be 
extended to disturbances with any number of 
wave components. 


e. Initial disturbances with tilted troughs. If 
this effect of the earth’s rotation is always in 
the same direction regardless of the initial tilt 
of the disturbances, then we may say that it is 
not a temporary effect. Thus we shall consider 
the initial disturbance 


Wo=4,F,(y) sin mA+a,F,(n) cos mA (2.26) 


d 
with F,ZkF,. For this disturbance, = takes 


the form 
oy 


ES 2mAla,yı(n) cos mA - agyy(m) sin mA] + 


+ Z, sin 2mA+ Z, cos 2mA+ Z, 


in which y, and y, satisfy eq. (2.19), with the 
proper F. Substituting in (2.10) we find 

I uv _ 2 dF, as 

Rmo a (71 dn ‘on 


Since y; is related to F, in the same way as 
y, is related to F,, we see that both the sine and 
the cosine components produce an increase of 
the transport irrespective of the initial transport. 
Thus our result also applies to the most general 
disturbance given by (2.5), and the effect will 
persist. 


3. Modification of the mean zonal current 


Let us now assume that there exists a mean 
zonal current U=u(y) and study its inter- 
action with the disturbance. Writing the 
streamfunction as p(y) + p(x, y, t) and u for 
pair wel) , it can be seen that the terms in 

R cos? 9 dy 
(2.18) depending on the mean flow are given 


OV ?p dp 0 (ou : 
az ox à ax dy (5 ar qu), where ke 
2% 
now stands for the operator — + a Ylnkthe 
a AN 
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following, we shall neglect the term 2 sin @- u 
and the variation of cosœin the nonlinear 
terms, and consider the coordinates x, y as 
has been multiplied by R. Equation (2.18) then 
becomes 


2 2 UG ie 4 Ind 
oe el: | 


Me 
cos? p\dy dx dx dy 


(3-1) 


where the primes denote differentiation with 
respect to y. This equation has the same form 
as the equation for plane motion, except for 
f'. Since it is only an approximation to (2.18), 


we may still take f’ as given by = cos? gy. In 


what follows we shall investigate the interac- 
tions between different mean currents and the 
initial disturbance 


(3.2) 


where /=z/D and D is the width of the belt. 

As the general effect of earth’s rotation is 

independent of the presence of the mean zonal 

current and is already known, we shall find 

only its expression for this disturbance, by 
AR 


putting Yan + cos 29) =? (1 + x cos ly), 


Yo = a sin ly sin kx 


with «<1. The increase of the mean mo- 

mentum transfer produced by this effect is then 
equal to 

duv CON Ce 

dt 2R I u?+ 


in? ] ; 
u y (3598) 


which is positive everywhere and is very 
similar to that given by (2.14). It may be 
noted that if the amplitude of the disturbance 
is the nth power of sin/y, then an antisymmetric 
uv ; 

Er which 
originates from the constant part of f’ and is 
given by M, of (2.22). From (3.3 a) we find 


part will also be present in this 


du Mere k2a? 
Ot R 4(u? + 4) 


(3.3 b) 


sin ly sin 2ly 


These quantities may be added to the corre- 

sponding changes due to the mean current. 
For the particular disturbance (3.2), the 

nonlinear terms on the right side of (3.1) 
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disappear and the solution a takes on the 
form 


— =aky(y) cos kx 


where y is to satisfy 
y" —key=[(ke + P)u tu" —f"]sin ly (3.4) 


and the conditions y(o) =y(D) =o. 
The rate of increase of the mean momentum 
transfer is given by 


Au) | a?k? 
die 


(ly cos ly—y" sin ly) (3.5) 


From this we get the second derivative of the 
mean zonal current 
Au Puv) a2k 


2 
— 3 [77 2 
at ems 1, 20.0) 


Since the present problem is related to the 
stability of the mean flow, we may also write 
down the energy equation for the mean flow. 
Since u,=0 initially, E,=o and 

D 

Eu = ed) (3.7) 
Thus the disturbance might be said to be 
“amplified” when E, is negative, “damped” 
when E,, is positive and neutral when E, is 
zero. All these quantities are related to each 
other and their determination depends upon 
the solution of the vorticity equation. How- 
ever, it can be shown that u,, not only occur 
in E,, but also in the higher derivatives. 
This can be seen from the increment of the 
kinetic energy of the mean flow relative to 
its initial value, 


D 
AE=, | (ut - Undy- | 
0 
; 3 | (3.8) 
= [unauay +, fau“ 
D D 


If AU is identified with Un irene be 


shown that the first term corresponds to 
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RE in the Taylor expension of E and the 

5 ! {4 me Thus 

a! ott d 

even when all the higher derivatives of U are 
4 


DE ES : : 
zero, a is still different from zero. Since this 


term is positive if Au is not identically 
equal to zero, it represents a damping effect 
upon the disturbance. Thus, E, as a measure 
of the initial rate of amplification or damping 
of the disturbance generally underestimates 
the damping effect and overestimates the 
amplifying effect as compared to the informa- 
tion given by U,. 


second term corresponds to 


4. The mean zonal flow U= U, + (U,,— U,) 
sin? ly 


Let us discuss this symmetrical mean zonal 
flow. In the expression U,, is the maximum 
value of U at the axis of the flow. The profile 
has two inflection points (where U’ =o) 
midway between the axis of the flow and the 
walls, and also critical points (where u” =f’) if 
Um is higher than a certain value. According 
to the previous investigation (Kuo, 1949), the 
inviscid linear vorticity equation will yield 
solutions with positive and negative real 
exponential time factors, so that there exist 
amplified and damped disturbances in the 
strict sense. We shall discuss this case in some 
detail to find whether the present method will 
also give comparable results. 

When this expression for the mean flow is 
substituted into (3.4), the function Y can be 
found by simple integration. Thus we find 
with the part depending on f’ omitted 


4y en sin 3ly — | 
k2 + ol En) 
(3k? — 22) À + 4(k + Plus Par | 
x +P ? 


The second time derivative of the mean flow, 

expressed in non-dimensional form, is then 

given by 

k2a?n2 A(u? - 3) 
(u? + 9) D? 


Uy = sin ly sin 3ly (4.2) 
where u=k/land A is used as an abbreviation 


for Un 7: (Ola. é 
From (4.2) it may be seen that the interaction 
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between the disturbance and the mean flow 
increases the mean flow in the central part of 
the zonal band from ly = 60° to ly=120° and 
decreases the mean flow in the rest of the band 
if the wavelength L=2x/k of the disturbance is 
shorter than 2D/V3. The change in the mean 
flow is in the reverse direction if the wave- 
length is longer than this critical value. 

When U, and the initial mean flow U are 
substituted into (3.6) we find 


4D 702 


FE (u2=3)u8 
a(k2 + 2) A® 


2 (u?+1)(u?+9 


(4.3) 


Since a?(k?+])/4 is the average disturbance 
kinetic energy per unit area, this equation with 
a minus sign expresses the normalized second 
time rate of change of kinetic energy of the 
perturbation. From (4.3) it will be seen that 


the wave length L=2D/V3 is the neutral 
wave length which separates the longer 
amplifying disturbances from the damping 
shorter disturbances. It may be noted that this 
rate of change has a maximum negative value 
(=- 0.496 A®/D?) for the wave length L=2.1 D 
so that this disturbance may be termed the most 
unstable disturbance. In a flow of U,,=20 m 
sec-1, of width D=5x10%m, the most 
unstable disturbance will increase in energy 
by but 2.96 % during the first day. However, 
if the mean flow is stronger and contained in 
a narrower belt, the change may be very 
rapid. Thus in a flow of U„=som sec~? of 
width D=2x10° m, this most unstable dis- 
turbance will double its energy during the 
first day. It will also be noted that the earth’s 
rotation does not contribute to E, with this 
symmetric mean flow profile. 


5. The higher time derivatives of the mean 
flow 


As the kinetic energy of both the mean 
flow and the disturbance is finite for a finite 
disturbance, there must be an upper limit to 
the energy attainable by the disturbance when 
it is amplified or by the mean flow when 
the disturbance is damped. Whether this limit 
is actually attained and the rate of approach 


to the limit depends not only upon U, or Ey, 
but upon higher derivatives as well. Thus, a 
discussion of the higher derivatives in one 
case is of interest. 


lel Silla (OY SILAN IN| 1 (O@ 


2 DE 
Generally, in order to find Se it is neces- 


sary to compute all the time derivatives of y 
up to the order (n-ı) which while possible 


in principle, would require an undue amount 


: Fu . ES: 
of labor if n is large. Since 38 is zero initially | 


4 
let us find only ce for the mean flow dis- 


cussed in the preceding section. For simplicity, 
we shall take f’ as a constant (=f) and put 
U, equal to zero. Then the fourth derivative 
is given by 
4(u?+9)?D? ou 
Braque = 3)A®n® 26 — | 


_ {e+ 9~ 2B) (3m? — 1-43 


| 2(u?+1) 
(„+9 28)? (#249) (3¢" = 7-4) 
2(u? +9) 8(u?+ 1)? 


(+ 5)(m? + 21) | aP(aut- 3 su? — 207)|, 
8(u? + 25) (u? + 4) A? | 
3(u?+ 5)(u?+ 9) 

x (cos 4ly - cos 2ly) + . 
a re 
‚ju?+9=2f’ u®+25=2P’_ 3u-1-4ß 

iP ees m?+25 w+1 


4a2l? (24 — 7u?- 69) 
(u? + 5)(u? + 4) A? 
le 3)(u? + 5) 
4(u? + 25) 


in which u=k/land 6’ =$/I?A. The expression | 
(5.1) is too complex to discuss generally. 
However, there are a few points worthy of 
mention. It may be noted that for the dis- 
turbance with u?=3, these derivatives are 
zero so that the disturbance is neutral at 
4 
d _ The factor À comes 
into this equation raised to a higher power 
than in the expression for U,, so the influence 
of the higher derivative is much higher with 
a stronger mean flow than with a weaker one. 
Computation shows that this fourth derivative 
generally has a sign opposite to that of Uy, 
indicating that the time required by the system 
to reach the limiting condition is longer than 
that estimated from U,, or the change may 
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x (cos 6ly — cos 4ly) + 


x (cos 8ly-cos 2ly) (5.1) 
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Fig. 2. The variations of U,, for the profile U= A sin? ly, 
with A= 20 m sec-1. (a) is for the damping disturbance 
u? = 5, (b) for the amplifying disturbance u? = 1.5. 


“y 
ot: 
its effect must be small for disturbances of 
long wavelength. 

To illustrate the effect of this fourth deri- 
vative, a few computed rates of change of the 
mean flow U, are shown in Figures 2 and 3, 
expressed in units of a?k®AD-*. Since the term 
containing & in (4.2) is neglected, the change 
is a symmetrical function of y so that only 
half of the flow is plotted. Figure 2 shows the 
effects of a relatively weak mean zonal flow 
(A=20 m sec-!) and Figure 3 shows the effect 
of a stronger mean flow (A=5so m sec-1). 
These examples show that the parabolic 

2 


be reversed. Since is proportional to kt, 


: Sn t 
extrapolation u —#)=—v, cannot be used for 
2 


time intervals longer than 24 hours. 

To estimate the order of magnitude of the 
interaction between the disturbance and the 
mean flow, we may take the example plotted 
in Figure 2a, which is for the damped dis- 
turbance with u?= 5, al= 10 m sec-1, A=20 m 
sec-1, and D=sxiofm. Taking the mean 
value of U, for a time interval of one day, we 
find an increase of the mean zonal flow of 
about 2m sec-! in the middle of the stream 
and a decrease in the outer parts. These changes 
are of the same order of magnitude as those 
produced by the earth’s rotation. 

As mentioned above, the disturbance with 
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(a) (b) 


Fig. 3. The same as Fig, 2, with Aso m seca 


RTE Saat cooled 
u?=3 is neutral at least up to — or =>. Since 


fg Gh 
it is very difficult to find the higher derivatives, 
it appears desirable to examine the exactly 
neutral disturbance from another point of 
view, if one can be found. For the present 
case, let us consider the total stream function 


A AU 
p= (ot S) y+ sin ay + 163 


+ a sin k(x — ct) sin ly 


which represents a system of alternating 
cyclonic and anticyclonic vortices embedded 
in a sine-curve mean zonal flow, propagated 
with a phase velocity c without changing 
shape, and therefore is exactly neutral. At 
t=o, this reduces to the initial motion we 
have been discussing. In order that this stream 
function satisfy the vorticity equation (3.1) 
with f’=, the wave number k and the phase 
velocity c must satisfy the following relations: 


2— 3/2 


Auden ie 
2 +P 


(5.3 a) 


C= Uy (5.3b) 


It is seen that the wavelength and phase 
velocity of this neutral disturbance are deter- 
mined by the mean zonal current and the 
width of the belt, so that it cannot be con- 
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sidered as representing a continuous spectrum 
of neutral disturbances unless A is zero, or 
unless the mean flow is uniform. For the latter 
case, the condition (5.3 a) is unnecessary. 

Since the stream function (5.2) reduces to 
the initial flow we have discussed, it may be 
concluded that the disturbance of the type 
(3.2) with w?=3 is exactly neutral for this 
mean zonal flow if the variable part in f’ is 
neglected. 


6. The mean currents 
U = U, sin! ly 
The simple sine profile U=U, sin ly has 
no inflection point within the walls, therefore 
is stable for small disturbances according to 
the linear theory. With this mean flow, the 
solution Y of (3.4) is 


ad [> + Snk(D - y) 
ur +4 SnkD 


U=U,, sinly and 


u? 2 y | 
1 - — sin?ly 
(6.1) 
Substituting into (3.6), we obtain 
Di Ir = Tor (he ar 1) 


aU,  u?+4 
r Snky + Snk(D - y) 3u? sin ?ly 
co SnkD © 2(u2+ 1) | 
(6.2) 


It may be seen that for any wavelength, this 
expression is positive in the central part of the 
zonal belt and negative in the outer part, 
indicating that the interaction between the 
disturbance and the mean flow is to accentuate 
the mean flow while the disturbance is damped. 
In this case, the energy integral E; is given by 


4D 


2 E = 
a2(k2+P)U,, 
2,,2 2 
= (2 | NT =| 
u? + 4Lu(u? + 4)z0 2 16(w?+1) 2 


(6.3) 


which is positive for all values of k, indicating 
that all the disturbances of type (3.2) are 
damped. 

To find whether this is really the case, we 
might try to find an exactly neutral disturbance 
in this mean zonal flow. The only disturbance 
that can propagate with a constant phase 
velocity without change of shape is the one 
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Fig. 4. The distribution of U,, for the mean flow U = U,, 
sin ly. 


of k=o, that is, only an infinitely long dis- 
turbance can remain exactly neutral and per- 
fectly periodic. Since it has been shown that 
disturbances of finite wavelengths are damped, 
it must be concluded that all disturbances of 
the type (3.2) must be damped in this mean 
flow. This is in agrement with the result of 
the linear theory. 

The change of the mean zonal flow produced 
by the interaction between the disturbance and 
the mean flow, with L=D, is given by the 
curve U, in Figure 4, while the other curve re- 
presents the initial mean flow expressed as the 


ratio (U- U)/ Us, where U,, is the maximum 


and U the average value. From these two 
curves, it may be seen that the kinetic energy 
of the mean flow is increasing and the dis- 
turbance is being damped. For a disturbance 
with L=D, the relative second derivative of 
the kinetic energy given by (6.3) is equal to 
0.5776 U;,/D2, so that for U,=20 m sec-! 
and L=D=sx10° m (about 45 degrees of 
latitude) the kinetic energy changes by 3.45 % 
during the first day. 
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Fig. 5. The distribution of U,, for the mean flow U = U,, 
sin@]y. 


As a further example of a sine-curve profile, 
we consider the mean flow U=U, sint ly, 
which has two inflection points at y= D/3 and 
y=2D/3, and therefore may yield both 
amplifying and damping disturbances ac- 
cording to the linear theory. It can be found 
that for this profile U, and E, are given 
respectively by 


eee sin | 
kaU, 4 7 
ae 2 
[peas HET ME aly | 
[PRP ZS, ur9 
4D Te 1182-945 
3 En Bris 2 2 
aX(k2+2)U2 32° (uw? +1)(u?+9)(u? +25) 
(6.5) 
According to Ey, the disturbance with 


wavelengths longer than 0.8322 D are ampli- 
fied and those shorter than this value are 
damped. The most unstable disturbance is the 
one with L=1.7D, for which the energy 
integral has the value —0.5. The changes of 
the mean flow produced by the interaction 
Tellus V (1953), 4 
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with the disturbances with L=D (amplified) 
and L=2 D/3 (damped) are shown in Figure 5. 


7. The polynomial profiles 


If the mean zonal wind exhibits a profile 


of the third degree 


U=a,+a,y+ Agy? + 4379 62) 


and a disturbance has the form (3.2), the 
solution of Yin (3.4) is 
AT) De 
ee sin ly — 
2 21 I (7-2) 


l 
~ eg BOW) 008 I~ ED 
-[G(D) Snky - G(o) Snk(D - y)] 


, 3k? ri lo m 
where G(y)=U’(y) Fe aU (y). We shall 
apply this solution to three particular cases: 
a. The constant shearing flow U= U'y (Plane 
Couette flow), where U’ is a constant. Since 
the profile has no inflection point, according 
to the linear theory, it is stable for small 
disturbances. In our development, 


-[Snzzu(1-n) — Snuren] sin sen 


where U„=U’D, the maximum value of U, 
n=y/D and u=k/l. Then the energy integral 
Ey, is given by 

4DE, pe ‘ 
a2(k? fe 208, Qu? = 1)3 4 Lu os 1)(u? 2 3) ur 


4a D°u? 
+i)(+4aU, 


+ he coth od a (7.4) 
7 21 3(u 

which is positive for all values of x, indicating 
that all disturbances are damped. In (7.3) it 
will be seen that if U’ is positive, the effect is 
to increase the zonal wind in the northern 
half of the belt and to decrease it in the 
southern half, the same direction of change as 
that produced by the earth’s rotation. It will be 
noted that with this profile « also enters into 
the energy integral E,, increasing the kinetic 
energy of the mean flow. Since this profile is 
of small interest, no numerical example will 
be given. 
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0.4 The distribution of U, produced by the 
interaction between the disturbance and the 
mean flow, given by (7.5), is represented in 
Figure 6 for L=1.155 D. The relative second 
derivative of the energy of the disturbance 
given by (7.6) is 0.4875 Umn?/D® for a dis- 
turbance of this wave length. It may be noted 
that this case is quite similar to that of the 
simple sine curve profile discussed in section 6. 

c. The profile v=27 u„n?(1-n)/4. This is 
an unsymmetric mean flow profile which has 
only one point of inflection, at y=D/3, as is 
shown in Figure 7a by the curve marked 
(U-U)/U,. The perturbation theory shows 
that if the effect of the earth’s rotation is 
neglected, or if the mean flow is strong, both 
damped and amplified disturbances exist, 
separated by a neutral disturbance. Because of 
the presence of a singularity in the linear 
vorticity equation, even the neutral disturbance 
is difficult to obtain, making it most desirable 
to study this case from the present point of 
view. 


With this profile, we find 


Dit ei 277 6 Wee ai 
PAU, 8 [re Ta 


-2.0 -1.0 fe) to. 0. “2:0 


Fig. 6. The distribution of Uy, for the mean flow 
U= 4U yy (I—7)- 


+ sin 270 re, — 3n) sin? 

b. The parabolic profile U=4U,n(1-n) Der LI Has per 
(Poiseuille flow). According to the linear an k — 
theory, this form of profile is also stable IR ns um + — a ae 
for small disturbance, since U” is equal to 4 Snr ror (i? He) 

-8 U,,/D? and there is no inflection point in [Sine See Ge 


the profile. From (7.2) we find 


DU, a4 The energy integral is given b 
Br = 2n(I — 27) sin 277 + au — 3) 3) > Sy 5 8 Y 
Ur wert DET (: =) + (u? — ı)n2 
i DEN 5 
* sin? gon + es [Snur(1 — 9) + Snuxn] RU, 3 16 2m*}  360(u? + 1) 
4 (7-5 Faß) | 
The energy integral E, is given by 27? 27. 3(u? + 1)? 
SDE as (1 BR: Ga. . E He + 4 re TS 
a:(k?+B)U), 3 eu (He ENT)E 167? mn? (u?+1)5 6(u2+ 1)? 


2 De 6 2, — 22e. 
(x ee D” +38 =] SH AGE Le 
TT 


= m? (u2+1)8 (u? + 1)? (ea 
256 u° HS cschum + © De | + a 
. a + — = . 
a (7.6) eS we LT nue + ne 
which is positive for all values of u, so that all en awD* (40 
disturbances of the type (3.2) are damped. Du 35R(u? + 4)\ x v8 (7.8) 
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Fig. 7. The distribution of U, for the mean flow U = 27U,,n? (1—1)/4. 


(a) is for the longer disturbance L= 3.464D, (b) for L= D. The curves I represents the results when the 
earth rotation factor @ is neglected, curves 2, when the effect of w is included. 


It will be noted that the parameter w of the 
earth’s rotation enters into the last term of 
this energy integral because of the lack of 
symmetry of the mean flow. This term is 
always positive and therefore represents a 
damping effect. 

Two examples worked out for this mean 
flow are plotted in Figures 7 a and 7 b. Figure 
7 a is for a longer disturbance of L=3.464 D 
and Figure 7 b is for a shorter disturbance of 
L=D. To demonstrate the relative importance 
of the mean flow and of the earth’s rotation, 
the quantity D?U;,/k?a?U,, has been computed 
both without and with the part given by 
(3.3 b). The results are represented by the 
curves I and 2 respectively. Curve 1 in Figure 
7 a shows that when this part is neglected, the 
interaction between a long disturbance and 
the mean flow diminishes the mean flow in 
the northern part of the flow where U is 
relatively high, and increases the mean flow 
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in the southern part where U is relatively low, 
and so abstracts kinetic energy from the mean 
flow and flattens it. The relative second deri- 
vative of the kinetic energy of the disturbance 
is —0.1062 U„/D?. When the effect of the 
earth’s rotation is included, this amplifying 
effect may be reversed. Thus, with D =5 x 10°m, 


“xO 
= MO io m 


R 
sec-!, the last term on the right side of (7.8) is 


equal to 0.2308, so that the total relative second 
derivative of the kinetic energy is now equal 
to 0.1246 U,,/D?. The corresponding modifi- 
cation of the mean flow is shown by curve 2 
of Figure 7 a, which shows a general accentua- 
tion of the contrast of the mean flow for 
different zones. 

For the shorter disturbance with L= D, both 
the mean flow and that of the earth’s rotation 
have a damping effect so that the inclusion of 
the term with m only magnifies the damping 


sec-!, and U,,=20 m 
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Fig. 8. Same as Fig. 7 for the mean flow U=27U,n (1—n)?/4. 


effect, as indicated by the curves 1 and 2 of 
Figure 7b. The relative second derivative of 
the kinetic energy of the disturbance is 0.6955 
Un/D? without the parameter m and 1.0955 
U;,/D* with this parameter of the earth’s rota- 
tion. 

It should be pointed out that the earth’s 
rotation may have an amplifying effect if the 
mean flow is unsymmetrical and is stronger 
in the southern part of the flow and weaker in 
the northern. To demonstrate this, let us 


27 Un 


take the mean flow U= (1-1)? which 


is the image of the profile (c) about the line 
n=— The change of mean flow is obtained 


by replacing 7 in (7.7) by 1-7. Figures 8 a 
and 8b show the result of this substitution 
for L=3.464 D and L=D respectively. The 
energy integral is obtained from (7.8) by 
changing the sign of the last term containing 
the factor w. For these disturbances, the energy 
integral including the last term is — 0.337 and 
0.295 respectively. 


8. The mean meridional transfer of zonal 
momentum 


The mean transfer of zonal momentum 
produced by the large scale atmospheric 
disturbances has been investigated by a number 
of workers in recent years, and its meridional 
distribution is now roughly known. It will 
be most interesting, then to compare the 
theoretical results with that obtained from 
observational material. 


In this section we shall discuss the mean 
momentum transfer corresponding to the 
disturbance (3.2) and the mean flow profiles 
(a) U=U, sin*ly, (b) U=U, sinly, (c) 
U=27 Unn?(1-n)/4 and (d) U=U, sin ly. 
For these four mean flow profiles, the rate of 
increase of the mean meridional transport of 
zonal momentum is given by the following 
equations: 


D dw (u?-3)x, . 
FEU, ot” Bue eo) (2 sin 277 — sin 477) 
(8.1) 
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To each of these we should add the transport 
produced by carth’s rotation, given by (3.3 a). 
For the mean flow profiles (a), (c), and (d), 
the nature of this interaction depends upon 
the wave length of the disturbance. However, 
since the mean zonal flow in the atmosphere 


‘ds generally not very strong, this interaction 


represents in the mean a horizontally damping 
effect on the large atmospheric disturbances, 
therefore we shall only discuss the damped 
disturbances with L=D for (a), (b), and (c) 
and L=2 D/3 for (d). It may also be pointed 
out that for these large-scale disturbances, this 
average meridional Be Buion of the tänsfer 
is not much effected by the occasional presence 
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Fig. 9. Theoretical distributions of uv for the mean 
flow profiles (a), (b), (c), and (d). 


of amplifying disturbances, as it is determined 
largely by the earth rotation effect. Taking 


0 
IDF OG MOP? a Ra 1.1x 107m! SCC 


and U„=20 m sec-!, the nondimensional 
rates of increase of the mean momentum 
transfer given by these equations are as 
represented by the curves (a), (b), (c), and (d) 
in Figure 9. Similar transports can be produced 
by earth rotation alone, if the disturbance is 
more concentrated as have been discussed 
before. 

Since the mean transfer of zonal momentum 
is initially zero for the type of motion we are 
considering, we may take the rates of change 
as representative of the actual transfer and 
compare them with the results obtained from 
observations. 

In Figure 10, the curve 1 is the six months 
average (January through June, 1950) total 
meridional transport of zonal momentum 
across unit length of the latitude circle, from 
earth’s surface to the 100 mb level, based on 


492 


-8 8 16 24 32 


10° 1 
-4 fo) 4 8 12 


1 
16 X10" cogs 


Fig. 10. Observed distribution of uv: I, 6 months average 

(January through June, 1950). Hemispherical. II, February 

1949, American Continent. II, January 1949, Hemi- 
spherical. 


the results obtained by Starr and White from 
actual wind studies along the latitude circles 
13° N, 31° N, 42° N, 55° N and 70° N (Starr 
and WEITE, 1051; 1952; 1952 b; and un- 
published computations). Curve 2 is the 
average total transport over the American 
Continent in February, 1949 (STARR, 1951), 
and Curve 3 is obtained from Mintz’s (1951) 
results of geostrophic transport for January, 
1949. The units for the two curves 2 and 3 
are twice as large as for the curve 1. If we 
put the maximum point of the curves in Figure 
9 at the maximum points of the Curves 1 and 
2 of Figure ro, it can be seen that they are in 
rough agreement. The magnitude of the 
transfer represented by the curves in Figure 9 
is also correct if we take the mean north—south 
velocity to be about 5 to 7m sec-1, and let 
the increase take place for a period of one day. 
For the month of January, 1949, Mintz’s mean 
geostrophic zonal wind profile indicates that 
two jet streams of westerly current were 
present, a stronger one around 30°N and a 
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weaker one around 50° N. This suggests that 
the curve 3 in Figure 10 may be interpreted 
by superimposing two sets of disturbances 
with the one having a smaller intensity shifted 
about 20 latitudinal degrees to the north. 
Thus, if we take the curve b in Figure 9 and 
add another one with an intensity equal to 
two-thirds of that of b and shifted a distance 
equal to D/3 toward north, we get the dashed 
curve in Figure 1. It is seen that this curve fits 
the curve 3 in Figure 10 nicely. The momen- 
tum transfer produced by the interaction be- 
tween the disturbance and earth’s rotation 
alone without any mean flow, is also shown 
in Figure 1 by the full curve. 


9. Summary and Conclusions 


The problem of the production of mean 
zonal flow by atmopheric disturbances is 
attacked by the computation of tendencies. 
The investigation shows that the effect of the 
earth’s rotation is always to produce a mean 
northward transport of westerly momentum 
and a tilt of the troughs and ridges in the 
SW—NE direction over almost all the regions 
except near the northern boundary, with the 
maximum transport and maximum tilting 
occurring south of the mid-point between the 
two latitudinal walls assumed. This interaction 
may either be a damping effect or an ampli- 
fying effect for the disturbance, depending on 
the asymmetry of the existing mean flow 
profile. For the symmetric mean flow profiles, 
this process always results in a slow damping 
of the disturbance, and appears in the fourth 
time derivative of the mean flow energy. 
From the nature and distribution of this 
effect, it would appear that it must have 
some bearing on the meridional distribution 
of the mean northward transfer of zonal 
momentum observed in the atmosphere. 

In regard to the interaction between the 
disturbance and the mean flow, the six examples 
examined indicate that the presence or absence 
of an inflection point in the mean flow profile 
is an important and decisive factor, so far as 
the particular disturbance is concerned. Thus, 
in the three cases in which U” is of one sign 
throughout, all disturbances of the type 
discussed are damped, while in the cases in 
which U” changes its sign, both damped and 
amplified disturbances exist. The disturbances 
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of wavelength longer than a critical wave- 
length abstract kinetic energy from the mean 
flow and flatten the profile, and are therefore 
amplified, while disturbances of shorter wave- 
length are damped, accentuating the mean 
flow. The three sinecurve profiles discussed 
above also indicate that the range of the ampli- 
fying wavelength increases as the inflection 
point of the profile moves away from the 
boundaries. 

Since the mean zonal current in the atmos- 
phere is generally not very strong, the in- 
teraction between the disturbance and the 
mean flow generally represents a horizontal 
damping process similar to the effect of 
earth’s rotation, while the possible instability 
of the horizontal motion associated with a 
stronger mean zonal current happens only 
occasionally. 

To compare our theoretical results with 
observational facts, the meridional distribution 
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of the mean northward transfer of zonal 
momentum are computed by this method for 
four different mean velocity profiles. These 
examples show that the effects of earth’s 
rotation and that of the mean flow are of the 
same order of magnitude for the large-scale 
atmospheric motions. When taken together, 
the combined effect is to produce a mean 
northward transport of zonal momentum 
corresponding to a general SW—NE orienta- 
tion of the troughs and ridges in almost all 
parts of the belt except a possible reversal 
near the northern boundary. It is shown that 
if the mean north—south velocity is about 
5 m sec-! and the increase of these transports 
takes place for a period of one day, then the 
computed transports agree roughly with the 
results obtained from observational material, 
indicating that the processes discussed are of 
real importance for the development of the 
mean zonal current. 
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Note Concerning the Nature of the Large-Scale Eddies in the 


Atmosphere 


By VICTOR P. STARR, Massachusetts Institute of Technology, Cambridge, Mass. 


(Manuscript received August Io, 1953) 


Abstract 


The rate of generation of kinetic energy of mean zonal motions by the large-scale horizontal 

eddies in the atmosphere is computed from three independent sets of hemispheric wind data. 

The results are larger than certain analogous dissipation estimates would require. Meteorological 
implications of the material are discussed. 


General discussion 


The aim of this brief article is to discuss 
further, with the aid of new additional data, 
some essential properties of the large-scale 
horizontal eddies in the atmosphere, which are 
involved in their relation to the mean zonal 
motions present. The particular physical 
principle to be reviewed has been treated 
previously mainly by my colleague Dr. H. L. 
Kuo (1950, 1951 a, 1953) both theoretically 
and observationally, although it appears that 
certain additional remarks and graphic illustra- 
tions might be of value, since some residual 
misunderstandings seem to be prevalent. 

The mean zonal motions (i.e., averaged 
with respect to longitude) of the atmosphere 
may be pictured as comprising a large vortex 
about the polar axis in which the angular 
velocity is a function of elevation, latitude and 
time. As an example the solid curve in (a) of 
the figure gives the distribution with latitude 
of this angular velocity relative to the earth, 
averaged also with respect to pressure in the 
vertical from 150 to 900 mb and with respect 
to time over the six-month period January to 
June 1950 inclusive (N. H.). Most of what is 
to be said concerning this profile can be made 
to apply equally to an instantaneous profile at 
some given level. One may now ask what the 


effect of genuine lateral friction alone on such 
a flow might be. Quite obviously the effect of 
such viscosity would be to retard the zones of 
most rapid rotation and to increase the angular 
velocity of the less rapidly rotating ones, so 
as to cause the whole to assume a more nearly 
uniform angular velocity. That is to say, 
lateral friction would cause a flow of angular 
momentum northward and also southward 
away from the zone of most rapid rotation. 

One may inquire next as to whether the 
large-scale eddies in the atmosphere act upon 
the mean zonal motions in a manner analogous 
to friction and hence may be regarded as 
merely giving rise to an eddy viscosity. This 
is outstandingly not the case, for according to 
any compilation of data at all adequate for the 
purpose (and there are a number of such 
independent ones now available) there exists a 
strong northward eddy transport of angular 
momentum to the south of the jet of westerlies 
(N. H.) and therefore a fortiori to the south of 
the zone of maximum angular velocity, on 
the average. This is exactly contrary to a 
frictional effect and would be compatible only 
with a negative virtual viscosity coefficient, if 
it were desirable even:to use such terminology. 
The mere existence of this state of affairs is in 
itself most remarkable and by no means to be 
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lightly dismissed, for regardless of further 
considerations it is a manifestation of a funda- 
mental process which requires explanation 
ultimately. The associated angular momentum 
transport by eddies is given by the dashed line 
in (a) of Fig. ı. 


(b) 


Fig. 1. Relative angular velocity in 1077 sec~ as function 
of latitude given by full curve in (a). Total angular 
momentum transport by eddies in 10% c.g.s. units 
given by dashed curve in (a). Production of zonal kinetic 
energy using mean wind and total transport (full curve), 
using 300 mb wind and transport (dashed curve) and 
using 200 mb wind and transport (dotted curve) shown 
in (b) in terms of arbitrary units. 


Kinetic energy considerations 


It is convenient and instructive to study 
these effects of the eddies upon the kinetic 
energy of the mean zonal motions. Considering 
the northern hemisphere as acting independ- 
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ently, the net influence on this total kinetic 
energy is to increase it, if the (positive) action 
to the south predominates and to decrease it 
if conditions are the reverse. The rate of 
production may be written as 


, t Ÿ [u] 

2 2 

az | fo[u'v'] cos Bl; Ar AL (x) 
if small approximations such as the horizontal 
uniformity of density are allowed in the 
definition of mean zonal kinetic energy. The 
notation is that of Starr and WHITE (1951, 
1952), u, v, being eastward and northward 
wind components, o density, r earth’s radius, 
® latitude, z elevation, the square brackets 
denoting averaging along complete latitude 
circles, and the primes indicating deviations 
from such averages. This expression may 
either be set up directly from a consideration 
of the work done by net forces due to the eddy 
stresses on elementary rings of air, or it may 
be obtained by writing the proper form of 
the balance equation for zonal kinetic energy 
as done, for example, by Kuo (1951 a). 


Data and computations 


Among the available data for investigating 
expression (1) the following may be listed. 

(1) The first six months of the year 1950 
already referred to. These data are based upon 
(daily) actual wind observations similar to 
those of Starr and White in the references 
already given, for Latitudes 13°, 31°, 42°, 55°, 
and 70° N. The integrand of (1) evaluated 
from the six-month averages of [u] and [w’v’] 
is given for the 300 and 200-mb levels by the 
dashed and dotted lines respectively in (b) of 
the figure. The full line results from using 
the vertical pressure averages of wind and 
eddy transport (i.e., corresponding to the 
curves given in (a) already discussed). Since 
in (b) the area under each curve measures the 
production of zonal kinetic energy, it is to 
be noted that in each case this area is positive 
by a wide margin. 

(2) The second six months of the year 1950. 
These data are exactly analogous to those for 
the first six months and were treated in the 
same way. The integrands (not shown) were 
again positive without question, i.e., do not 
represent small differences between large 
quantities. 
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(3) One month of data (Jan. 1949) published 
by Mintz (1951). These represent geostrophic 
evaluations obtained by punched card methods. 
As in the case of the curves (b), here also there 
is no ambiguity as to the sign of the areas 
beneath the curves (not shown) in all three 
instances. 

Using the vertically averaged winds and 
transports the data (1) give 4.2 x 107° ergs per 
second for the production over the hemisphere. 
The data (2) give 4.6 x102° units. Data (3) 
give 10.5 x 102° units. These magnitudes are to 
be compared with appropriately defined and 
properly estimated measures of the dissipation. 
When the vertical averages are used one deals 
essentially with the kinetic energy of the mean 
rotation of the atmosphere, and the dissipation 
becomes merely the work done against the 
zonal component of ground stresses in virtue 
of this mean relative rotation, save for minor 
other effects. Further discussion of this subject 
is given by Srarr and WHITE (1953). The 
zonal ground stresses may be obtained in the 
form of empirical climatological estimates as 
given by PRIESTLEY (1951), for example. The 
stresses deduced from the angular momentum 
balance considerations for the particular period 
in question as calculated for example by 
Mintz (1951) would not give an independent 
check for the present purpose. Use of Priestley’s 
stresses with the data (1) gives 3.6 x 102° units 
for the hemisphere. The data (2) give 1.9 x 1020 
units, while data (3) give 4.0 x 1020 units. 

If the angular velocity curve in (a) of the 
figure is assumed to be representative of the 
entire depth of the atmosphere, we may 
compute from it the corresponding kinetic 
energy of mean rotation for the northern 
hemisphere, which turns out to be about 
6.8 x 1076 ergs. Upon appeal to the figure 
given for the mean rate of production from 
data (1), it is seen that this amount is normally 
generated in about 19 days. The corresponding 
period for data (2) is 12 days and for data (3) 
it is 17 days. Presumably if the production 
were suddenly to cease while the dissipation 
were to maintain its normal rate, the motions 
of the air here involved would be abolished in 
about two weeks. 


Discussion of results 


Relegating certain further technical matters 


and qualifying remarks to the next section, 
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we may enumerate certain important inferences 
which may be drawn from the general char- 
acter of the results, if these latter are correct in 
main outline, as seems most probable. These 
are: 

1. The large disturbances in the atmosphere 
produce kinetic energy of mean zonal motions, 
and not the reverse as has been with few 
exceptions generally supposed in the past. 

2. Since the magnitudes of this generation 
are of the order of the estimated dissipation, 
and could regenerate the motions involved in 
about two weeks, it may be stated that any 
theory for the general circulation of the 
atmosphere which omits this effect is crucially 
deficient, lacking an essential correspondence 
to reality, however interesting it may other- 
wise be. 

3. Through the years there has been current 
in meteorology an assortment of general 
circulation theories which do not presuppose 
this action of the eddies. The singular sterility 
and lack of promise which have characterized 
these can perhaps be ascribed largely to the 
shortcoming under discussion. 

4. It is increasingly clear that a valid formal 
theory of the general circulation showing how 
the incoming radiational energy supply results 
in the observed motions (rather than others) 
will probably be evolved only slowly because 
of its inherent complexity. 

The problem of the mean zonal motions and 
the large-scale disturbances is one and in- 
separable, because of the now demonstrated 
interdependence of these circulations; hence 
there exists a natural and direct connection to 
such other modern developments as those in 
cyclone theory and similar topics, which latter 
will first have to be not only better understood 
but also fitted into the more general picture. 
For the time being it may be necessary to 
content oneself with the theoretical examina- 
tion of only portions or phases of the total 
process. Indeed much of what now may 
appear as composed. of separate and distinct 
subjects in theoretical meteorology must 
ultimately form a coordinated and harmonious 
whole, displaying in true light the proper and 
natural relationships among these parts. 

5. The theoretical efforts of my esteemed 
friend Dr. H. L. Kuo may be mentioned 
specifically as an example of what has just 
been said (see Kuo, 1950, 1951 b, 1953). These 
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studies are directed toward the elucidation of 
the transformation of eddy kinetic energy into 
mean zonal kinetic energy through the 
application of the vorticity equation. This 
procedure may be interpreted also as an effort 
to apply the general principles, if not the 
specific techniques, used in what is now called 
numerical forecasting to the problem of the 
general circulation. It is to be noted that in 
this connection, and also from the general 
results stated in this paper, the seat of the 
kinetic energy processes is placed in the 
cyclonic and anticyclonic systems in the 
atmosphere, which directly or indirectly 
account for the continued replenishment of the 
eddy kinetic energy. 

6. It may with profit be again contemplated 
that the positive nature of the production 
process as measured from data is not contingent 
upon the attainment of some petty accuracy, 
but is at once fixed by very gross aspects of 
the zonal-wind and angular momentum trans- 
port profiles. The end result could, in fact, 
hardly be otherwise if the maximum north- 
ward angular momentum transport is found in 
the vicinity of the mean latitude of the jet 
stream. 

7. Most of the so-to-speak simple theories 
of the general circulation start from the assump- 
tion that the mean zonal motions are main- 
tained by mean meridional circulations. From 
time to time thinking men have expressed 
dissatisfaction with such a framework to 
various degrees, e. g., JEFFREYS (1926), RossBy 
(1947, 1949), although the direct measurement 
of the mean meridional circulations cannot be 
carried out with great enough confidence for 
a direct refutation which would be reliable. 
Once, however, very extensive independent 
compilations of data consistently point to the 
existence of a sufficiently potent alternative, it 
is quite unscientific to concentrate merely on 
the older approach. The writer’s opinion in 
this regard does not exclude the possibility 
that sufficiently extensive observational studies 
may show that mean meridional circulations 
perform certain other important functions, or 
even modify to some extent the eddy-produced 
zonal flows. In any case it would be highly 
improper not to give recognition to those 
indications and suggestions deriving most 
clearly and unambiguously from purely objec- 
tive data as already expounded. 
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Critique and further computations 


During the course of the work performed a 
number of added questions and matters of 
thoroughness were examined and treated in 
much greater detail than can be fully entered 
into for the purpose of this note. Some critical 
statements and results of corresponding checks 
in some instances may nevertheless be made 
for purposes of better orientation. 

1. The use of vertical averages for the zonal 
wind profile is a rather specialized device 
which introduces valuable simplifications, but 
general interest still demands that the results 
be compared with what would be obtained 
from calculations at individual levels which are 
then integrated vertically as a final step. Such 
alternative figures were calculated giving for 
datar(T)13.9%x 102%, Moridat (2) sboro and 
fondat Mo so here sec-i; 

2. À second check carried out was to cal- 
culate the production using a vertically 
integrated zonal-wind profile but applying the 
technique to the instantaneous daily conditions, 
the long term mean being then formed by 
averaging the daily figures for the production. 

By this procedure an opportunity is also 
afforded to compute confidence limits defined 
as twice the standard error of the mean of the 
daily values, indicating roughly the 95 per 
cent confidence range. The data (1) were 
subjected to such a treatment, although at the 
time this work was done the observations for 
70° N latitude had as yet not been tabulated 
and hence were not included. The figure 
obtained was (4.5--1.0) x 102? erg sec-!. 

3. A still further elaboration consists of the 
use of the analysis by levels, not with the long 
period averages, but with instantaneous daily 
data. The long term mean may then again be 
formed by averaging the daily production 
figures. The data (1) with the 70° N latitude 
circle omitted were used also for this pur- 
pose. The figure which resulted was 9.8 x 102° 
ere SEC" 

4. In the computations it was arbitrarily 
assumed that [u] and [w’v’] are zero at the 
equator. What the exact values of the momen- 
tum transport are at the equator is not known, 
but the use of more proper negative values of 
[u] would tend to increase the anticyclonic 
wind shear immediately to the north where 
the transport is positive and hence would 
probably increase the production. 
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5. Selectivity in the availability of wind 
observations favoring lower speeds probably 
leads to an underestimate of the production. 
It is to be noted that the integrand of expression 
(t) involves a triple product of quantities in 
the nature of wind velocity components, 
hence underestimates of the latter could 
produce a larger percentual underestimate here 
than in such quantities as, for example, the 
kinetic energy itself. 

6. It should be most carefully observed that 
if the effect of the eddy processes on the mean 
zonal kinetic energy is studied for a narrower 
belt of latitude, terms representing the merid- 
ional transport of this energy across the 
bounding latitudes should be considered, and 
the mere adjustment of the limits of integration 
in (1) does not provide for this requirement. 
In the case of the whole hemisphere it is 
doubtful, however, that such a transport 
across the equator is anything more than a 
slight correction, which is here neglected (and 
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would necessarily vanish for the entire atmos- 
phere). 

7. The surface zonal stresses given by 
Priestley are according to most proper inter- 
pretation values for the oceans. Doubtless the 
extensive land masses present in the northern 
hemisphere give rise to larger such stresses as 
suggested by the angular momentum con- 
vergence study of Mintz (1951), thus helping 
to bridge the great disparity found here be- 
tween the eddy production and the dissipation 
of the energy of mean rotation. 

8. The expression (1) is actually nothing 
more than the spherical polar form of the 
corresponding term in the familiar theory of 
Ossorn Reynorps (1895) concerning the 
nature of eddy processes. A review of this 
material is given also by Lams (1932) in his 
textbook on hydrodynamics. Indications of the 
importance of the process here measured by 
expression (1) have been obtained also by 
BLACKADAR (1950) and by VAN MIEGHEM (1953). 
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Zur numerischen Wettervorhersage mittels Relaxationsmethode 


unter Einbeziehung barokliner Effekte. II’ 


K. HINKELMANN, Deutscher Wetterdienst, Bad Kissingen 


(Manuscript received January 20, 1953) 


U. Lösungsmethoden 


Es sollen 2 Lösungsverfahren kurz behandelt 
werden, die zur Berechnung der im Abschn. 
III. gezeigten Lösungen benutzt wurden. Jedoch 
sindauch andere Lösungsverfahren anwendbar.? 


1.) Relaxationsmethode nach R. Southwell: 

Zur Lösung des Differenzengleichungssy- 
stems (16), in welchem die obere und untere 
Randbedingung (7, 8) bereits enthalten sind, 
ist besonders -ür Handberechnungen die Re- 
laxationsmethode von Southwell geeignet. 

Man geht bei diesem Lösungsverfahren von 
einem an sich beliebigen bequemen Lösungs- 
ansatz, im allgemeinen von t =o (Lösung 
nullter Naherung) aus. Dano geben die rechten 
Seiten von (16) den Fehler R an, um den die 
Gleichungen durch den Ansatz t = onicht er- 
füllt sind. Diesen Fehler bezeichnet man all- 
gemein als Residuum. 

Durch schrittweise Verbesserung des t-Fel- 
des versucht man, das Residuum gleichmäßig 
im gesamten Gitterfeld gegen o konvergieren 
zu lassen. 

Als erste Nahrung t benutzt man gewöhnlich 
I 


Wea (0G ha OG 


1 0 . 
Tu ÔTk = Aghg mit Ak = 


In diesem Falle ist das zum Ansatz tT = 0 ge- 
hörige Residuum R mit dem inhomogenen 
Term der Gleichungen (16) identich. Da die 


1 Continued from Vol. s, No. 3. 


2 z. Bsp. die Verwendung des modernisierten Gauß- 
schen Algorithmus, vgl. dazu die Berichte Nr. s, 6 des 
Instituts für Praktische Mathematik, Darmstadt. 
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Gleichungen (16) einen Funktionswert 7 mit 
den Funktionswerten zi#"J, 14/41, 7%, , funk- 
tionell verkoppeln, bewirkt eine Änderung 
Ôt eine Anderung des Residuums nicht nur 
im Aufpunkt rj selbst, sondern auch in den je 
um eine Maschenweite vom Aufpunkt ent- 
fernten Nachbarpunkten. 

Im vorliegenden Falle ergibt eine Änderung 
dt, folgende Anderungen der Residuen. 


ORE =—Aj otf, OR,= = ÖRWEI= def (18 a) 
ORY = AR OÔTIN ORE =. drt (18 b) 


wie man sich leicht klarmacht, wenn man den 
Gitterpunkten rd, ri +, ru '=1, r,, die zu- 
gehörigen Gleichungen (16) zuordnet. Anstelle 
von (18 b) liefert jedoch entsprechend Gleich- 
ung (16 b) eine Anderung 

ét# im Punkte k='!/, die Residualänder- 
ungen ÔRŸ, ,= af, ,- dt; OR 4= 0 und analog 
Gleichung (16c) eine Anderung 

été im Punkte k=n—1/, die Residualänder- 
ungen 

DR ort; OR = 0 


Die durch eine Änderung von ör entstehende 
neue R-Verteilung kann wiederum zur Ver- 
besserung der t-Werte benutzt werden, ent- 
sprechend erhält man als zweite Näherung 
9 5 1 1 . 
t=t+ dr; mit ör= RA und gem. obigem 
Schema (18) neue Residuen R usw. 
iL 12 : 

Die obige Annahme Ôt} = AZ Ri bewirkt 
also gerade, daß das Residuum in Punkte (rj) 
exakt verschwindet. Allerdings wird durch 
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Verfügung über die t-Werte in der Nachbar- 
schaft von (rf) wiederum ein Rj erzeugt. 

Beim Relaxieren muß deshalb angestrebt 
werden, die Änderung von t im Aufpunkt 
(rÿ) und in den Nachbarpunkten so abzustim- 
men, daß R# erst durch die vom Aufpunkt 
und von den Nachbarpunkten herrührenden 
Änderungen möglichst klein wird. 

Eine schnellere Konvergenz wird z. B. er- 
reicht, wenn man anstelle der Abschätzung 


I 
ore ApR, die dem Lösungswert näherkom- 
mende Abschätzung 


N 1 Der I 
OT, = ApRyt+ AD Al Ry” + oe + Ages Revit 


(19) 


benutzt, vgl. dazu die Operatorenmethode in 
Abschn. II, 2 

Jedoch braucht man sich bei Handberechnung 
nicht an ein starres Schema fiir die Verbesserung 
der t-Werte zu halten. Durch geschicktes 
Schätzen der Funktionswerte t bzw. deren 
Zuschläge zur vorangehenden Schätzung unter 
Berücksichtigung möglichst vieler Residual- 
werte in der Nähe des Schätzpunktes kann der 
Lösungsprozeß u. U. erheblich beschleunigt 
werden. 


1 
+ &%- Ag Re) 


N 
Als Resultat erhält man schließlich 7 = 0. 

1=0 
wobei der ProzeB abgebrochen werden kann, 
wenn N einem Residualfeld entspricht, das den 
Genauigkeitsansprüchen genügt, d. h. dessen 
Residualwerte R sich innerhalb der Fehler- 

grenzen der inhomogenen Terme h halten. 
Die Behandlung der seitlichen Randnachbar- 
punkte richtet sich nach den speziellen seitlichen 
Randbedingungen. Im Falle, daß die t-Werte 
auf den seitlichen Randflächen vorgegeben 
werden, erübrigt sich selbstverständlich eine 
Residualwertkontrolle für die seitlichen Rand- 

punkte. 


2. Operatorenlösungsmethode. 


Mit etwa dem gleichen Rechenaufwand kann 
die Lösung von (16) auch durch Anwendung 
der Operatorenlösungsmethode herbeigeführt 
werden, wie sie z. Bsp. für die zweidimensional 
barotrope Vorticity-Gleichung von R. Fjor- 
Torr (1952) — vgl. Gl. (16)’ — benutzt wird. 
Gleichung (16 a) 


Te — Ag (E th? + or Terı + de Tei = Aghe 
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schreibt sich, wenn folgender Mittelbildungs- 
operator 


D=A (Dal Ho: pea pn) 
eingeführt wird: 
(t— D) t= Ah mit der formalen Lösung 


I 
iD) 


T= Ah (20a) 


Unter der Voraussetzung, daB die Koeffizienten 
der Iterierten von D hinreichend rasch nach 
o konvergieren — ob diese Voraussetzung zu- 


trifft, muß die Lösungsfolge erweisen — ist 
t= (I+ D+ D?+...+D!...4+) Ah= > D'Ah 
1-0 (20b) 


Wir setzen abkürzend ae = D!Ah; dann ist 
= Dör und öT =D (6) 
= 


Man erhält damit: 


One Ade 
1 i Ue 
6%} = Dy (6t) = Ap (E Oth" + 
1 l 
Xp OTp + fae ay ÔTk-1) 


NaturgemäB konvergiert dieses Operatoren- 
lösungsverfahren ebenso rasch wie der Relaxa- 
tionsprozeß, wenn man sich bei letzterem bei 
der »Schätzung« der t-Werte starr an das 
Schema (18) hält. 

Eine raschere Konvergenz wird durch fol- 
gende Entwicklung bei etwas veränderter Be- 


deutung der bt erreicht: 
ör= (1+ D) Ah 
ör = D? (1+ D) Ah= D? (ôr) 
[+1 


Ô T = D2(62) 


Als Lösung erhält man wiederum 


T 


I 


© | 
2.ö7 (21) 
= 

Das Lösungsverfahren erfordert die Bestim- 
mung des Operators D?. Durch Anwendung 
von D auf D; y — Ak (Cv se Ogi Wr+i Sf 
+ Op-Wr—1) erhält man 
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[Ar (Di)? pl = EAN (Zvi + af, plat 
+ of Va) + + Ansa (Dvir t+ ees Vere + 
+ +1 Pe) + &p- Apes (Spel + Ors Yet 
+ &-1- Pe-s) 


und nach Zusammenfassen von w-Werten 


gleicher Gitterpunkte. 
a Diy f= (EAM + AY a alla + 
+ Alta of.) yl ALES pa à 
ee ee anne 
ve (Ait TI + AE) pt hart (Ais hit 
us Aces) wo aie (Ari + Aa) eae + 
+ (AEM AL) Abd al al at 
+ Ab all adapte + © (AE oft + 


27 i,j i ID LÉ AS ig (à La yoke 
As 3.) Vr+1i 2, (Ai CHE A, ai) PER 


Beim Übergang von is auf Be zur Verbes- 
serung der Lösungsfunktion mittels der Folge 


À 1+1 x 
(21) gehen zur Bestimmung von 6 + bereits 18 


dem Aufpunkt benachbarte le ein. 
Allgemein führt die Entwicklung 

T= (1+D+D?+...), Ah, die alle h-Werte 
des Raumgitters in dem betrachteten Integra- 
tionsgebiet erfaßt, auf die Bestimmung der 
lösenden Kerngitterfunktion des Differenzen- 
gleichungssystems im Punkte ri, dargestellt 
durch die Koeffizienten der h-Gitterwerte. s. III. 

Die Berücksichtigung der Randbedingungen 
erfolgt dabei automatisch, wenn der Operator 
D an den Randnachbarpunkten entsprechend 
den Koeffizienten von (16 b, c) definiert wird. 

Praktische Rechenversuche mit Hilfe der Re- 
laxations- oder Operatorenlösungsmethode er- 
gaben, daß bei willkürlicher Streuung der in- 
homogenen Terme, wie das bei tatsächlichen 
Wettersituationen auftritt, etwa $—10 Relaxa- 
tionen bzw. Operatorenschritte für jeden Git- 
terpunkt durchzuführen sind, bis eine aus- 
reichende Genauigkeit der Lösungsfunktion + 
erreicht wird. 

Eine wesentlich schnellere Konvergenz wäre 
sicherlich herbeizuführen, wenn mit Hilfe von 
Einfluß- oder Greenschen Funktionen eine wenn 
auch grobe Vorschätzung des Tendenzfeldes + 
durchgeführt werden könnte. 
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Unglücklicherweise hängen jedoch die Koef- 
fizienten des Differentialausdrucks in 16 wegen 
der darin vorkommenden o, 7-Terme eben- 
falls wie die inhomogenen Glieder der rechten 
Seite von (16) von der zu behandelnden Wet- 
tersituation ab. 

Um eine Vorstellung zu vermitteln, mit 
welch unterschiedlichem Gewicht inhomoge- 
nes Feld h und Koeffizientenfeld « des Differen- 
tialausdrucks die Lösung beeinflussen, sollen 
die Gleichungen (16) hinsichtlich der von z 
und inneren Ableitungen von z abhängigen 
Größen linearisiert werden: 

Wir machen für das stationäre Grundfeld 
(überstrichene Größen) folgenden Ansatz: 


z= L(p) : y+ I(p); dann ist 
n=; weiterhin sei fx & 0 
Für das Stérungsfeld 2’; z= z+ 2’, sollen keine 


einschränkenden Annahmen gemacht werden. 
Linearisierung der Koeflizienten der Vorticity- 


Gleichung: 
a (i) 
LS Ges er zn 
N dei EE 
ergibt wegen © = T= 0 
2 
pl tere 


len) (22) 


Linearisierung der Koeffizienten des I. Haupt- 
satzes (6) mit 


dq 
dt 


[A à] 


ergibt wegen Q=J (%,zZ)=0 


ps etl > 
0; © à (t+ Q) 


( 
g 


[A 
Dir 


Seis Pe SS 
Aa Re — en) = f 
0 ae oe 
CR + 82pp» 
worin wegen der speziellen Wahl des Grund- 


zustandes die Akzente auch fortgelassen werden 
können. 
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Elimination von @ aus (22) und (23) und Er- 
setzen der 7 und w enthaltenden Differential- 
ausdrücke durch entsprechende Differenzen- 
ausdrücke ergibt anstelle (16a) 


(4+ Xp+ + ak) Te Oey Te+1—OXp—TR-1 — 
- De = hip Je (2 n°) + we (ef) + 
+9 (ox Rr+3 — &- Qy-z) mit 

2 u? I 


N 


= = für 1 <k<n—1 (24) 
AY OLLE 


und (16b) und (16c) analoge Gleichungen fiir 
k= 1}, und k= n— 3/5. 

Nach Linearisierung treten in dem Differenzen- 
operator von t d. h. auf der linken Seite von 
Gleichung (24), keine StérgréBen mehr auf. 
Das legt den Schluß nahe, daß man auch mit 
einem «-Feld eines Grundzustandes, der von 
der speziellen Wetterlage unabhängig ist, noch 
zu einer vernünftigen Lösung gelangen kann. 
Die Tatsache, daß jetzt wegen „= f 


ARE = X41, + 


läßt weiterhin vermuten, daß auch der verti- 
kale Vorticitytransport wy einen geringeren 
Einfluß auf die Lösung ausübt als der Divergenz- 
ausdruck 70, 

Wir werden also berechtigt sein, uns zumin- 
dest einen Überblick über den vermutlichen 
Lösungsverlauf zu verschaffen unter Zugrunde- 
legung eines von den speziellen Topografien 
unabhängigen Koeffizientenfeldes «. 

Dies hat den entscheidenden Vorteil, daß 
man mit Greenschen Standardfunktionen G 
ohne einen Iterationsprozeß zu eine rNäherungs- 
lösung +, w gelangen kann, wobei bei einer nur 
von p abhängigen Standardatmosphäre diese 
Greenschen Funktionen sogar weitgehend in- 
variant gegen Koordinatenverschiebung in x, 
y sein werden, zumindest in hinreichendem Ab- 
stand von der seitlichen Begrenzungsfläche. 

Mit Hilfe solcher ebenfalls mit Hilfe der Re- 
laxations- oder Operatorenlösungsmethode ge- 
winnbarer Funktionen G kann man sich durch 
eine einmalige Integration bzw. Summation 
die zu dem Standard-x-Feld gehörige Lösung 
eines beliebigen Feldes inhomogener Terme 
verschaffen, die als Vorschätzung der Lösung 
des aktuellen «-Feldes betrachtet wird. 

Nach Auswechseln der «-Koeflizienten er- 
hält man sicherlich nach wenigen Relaxations- 
bzw. Operatorenschritten eine dem variablen 
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a-Feld zugeordnete Lösung, wenn mit der mit 
Hilfe Greenscher Funktionen gewonnenen 
Näherungslösung begonnen wird. 


III. Anwendungen 


Die Fig. 4—13 geben einige Einflußfunk- 
tionen für verschiedene Lagen des Quellpunktes 
wieder: Wir wollen unter einer Einflußgitter- 
funktion I (r’, r) eine Lösungsfunktion 7 ver- 
stehen, die in allen inneren Gitterpunkten die 
Differenzengleichungen (16a—c) mit hij =o 
erfüllt, mit Ausnahme des Quellpunktes r’(j7), 
wo hii’ + o sei. Entsprechend lassen sich Ein- 
Alußfunktionen I als Lösungsfunktion der Gin. _ 
(16) definieren, wo im Quellpunkt r’ nur ent- 
weder der Vorticitytransport u?] (z,n) oder die 
thermische Quantität Q von Null verschieden 
ist. 


Im Falle 22 J (z,7) +0, Q=0 
(2, 9) 


fiir r=r’ ist 
Ae = a pr 


undh=o fürr +r’ 


Im Falle. Q + 0,y?J.(z,4)=0, für re 
WY y= — oll Qs WIG = i IC 
und h=o fiir r +4774, 


Ein »punktförmiger« Term Q# kann demnach 
in seiner Wirkung auf das Lésungsfeld durch 
einen Dipol im Vorticitytransport: 


LY (2, nas Tor POP 
ersetzt werden. 

Setz man fiir die inhomogenen Bestandteile 
im Quellpunkt einen Einheitswert an, so hängen 
die EinfluBfunktionen lediglich vom Koeffizien- 
tengitterfeld des Differenzenoperators nach (16) 
ab. Gewählt wurde im Quellpunkt füru?](z, 7) 
bzw. aus Dimensionsgründen für 


1 dq 


R 
degre fe 20) 


der Wert von 100 Einheiten [m. std-1| 
Da die Differentialgleichungen (5, 6) in +, w 
bzw. die Differenzengleichungen (16) in t linear 
sind, gilt das Superpositionsprinzip. D. h. bei 
beliebiger Verteilung der inhomogenen Terme 
gibt die auf h angewandte Einflußfunktion 
I(r, x’) im Punkte r(#) den Anteil des Lö- 
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sungswertes t in % des inhomogenen Terms 
h im Punkte r(7) an. Zur Berechnung der 
in der Fig. 4—13 dargestellten Einflußfunk- 
tionen wurde ein «-Koctlizientengitterfeld einer 
Standardatmosphäre zugrundegelegt mit einer 
von x, y unabhängigen z- und damit T-Ver- 
teilung mit folgenden Daten: 


oT 
Eh = 283,2°; Er = — 6°kmtfürp> 225 mb; 


= o für p< 225 mb; 
als Maschenweite wurde w= 300 km, »= 50 
mb gewählt.! 

Wegen der automatisch durch (16a—c) er- 
füllten kinematischen Randbedingungen ei- 
nerseits und der Abhängigkeit der Koeffizienten 
& von p, x, y hängt selbstverständlich der Ver- 
lauf der Einflußfunktion I (r',r) wesentlich 
von der Lage des Quellpunktes r’ ab. 

Wegen der relativ geringen Variabilität von 
& mit x, y gegenüber der Abhängigkeit von 
p wird die Einflußfunktion jedoch weitgehend 
invariant gegenüber einer horizontalen Ver- 
schiebung des Quellpunktes r’ sein, zumindest 
soweit der Quellpunkt noch hinreichend weit 
vom seitlichen Rand entfernt ist. Wir wollen 
uns deshalb bei der Berechnung von Einfluß- 
funktionen auf eine Quellpunktslage im mitt- 
lerer Breite (p in etwa 50° geogr. Breite) be- 
schränken und lediglich die p’-Koordinate des 
Quellpunktes variieren. 

In den Fig. 4—13 sind einige derartige Ein- 
fluBfunktionen in einem den Quellpunkt in y’ 
= 51° enthaltenden Meridionalschnitt x = 
const.? wiedergegeben, wobei einmal w?J(z, n) 
= 100 Einheiten in den Quellpunkten 


Pp’ = 975, 725, 475, 275, 75 mb resp. 


und einmal pQ= 100 Einheiten in den Quell- 
punkten 


p’ = 950, 700, 500, 300, 100 mb resp. 


angesetzt wurde. 


$03 


Obwohl Lösungen von (16) naturgemäß nur 
Gitterfunktionen mit definierten Werten an den 
Gitterpunkten allein liefern können, wurde in 
den Abbildungen aus Gründen der Anschau- 
lichkeit eine durch Interpolation gewonnene 
kontinuierliche Darstellungsweise gewählt. 

Da jeder Gitterpunkt ein Volumelement y?» 
repräsentiert, so approximiert das t-Feld der 
Fig. 4—13 demnach übergreifende räumliche 
Mittelwerte der Lösungsfunktion 7 der Diffe- 
rentialgIn. (s—6) 


Bi Le 
X+-y+-p+- 
ae 


T na 
TDifferenzengl. ri TDifferentialgl. dxdydp 
wy 
bh u » 


ii; 


Der Differentialausdruck L [rt], den man durch 
Elimination von @ aus (5, 6) gewinnt, ergibt 


sich zu 
e 
L{t]=At+ (2) 
[7] alee 
ui 
o 
druck ist dann selbstadjungiert®, wenn 7, = 0, 


Dieser wegen => 0 elliptische Differentialaus- 


3 Liegt eine (elliptische) Differentialgleichung der Form 


3 3 
L[u] 2a; Haze, DE. ut a, h; dj jp Ay 
i,f=1 i=1 


= Funktionen von x, 9, x3 


b= 


vor, so wird der Differentialausdruck L [u] dann als 
selbstadjungiert bezeichnet, wenn gilt: 


3 3 
Lu] = DO lan)... — Do (a;u),. + ay 4 
j=1 i i=l z 
Für einen selbstadjungierten (elliptischen) Differential- 
ausdruck läßt sich dann die Lösung in der Form dar- 
stellen: 

U (XXX) = —J u(s) G, ds + 


se ff) G (org, 84. 85; X a, x) (201, 88, 83) dx dx dies, 


wobei das erste Integral über den Rand des Integra- 
tionsgebietes, das 2. Integral über den Integrationsraum 
zu erstrecken ist. Der Kern G wird als Greensche Funk- 
tion bezeichnet und genügt der Symmetrieeigenschaft 


1 Bei einer stereographischen Projektion mit 9, = 60° erhält man für die Vertikalsäule in einer Breite von  # 51 


folgende Werte für die dimensionlose Größe «: 


LI — 
p = 1 000| 950] 900] 850| 800| 750] 700| 650| 600| sso| 500) 450| 400) 350) 300| 250| 200| 150| 100] 50 


a= 0,25| 41,6] 37,7| 33,9| 30,4 27,0| 23,8] 20,8] 18,0| 15,3| 12,0) 10,6) 8,6 6,7 


wobei wegen n=f pr = bis ist. 


s,1| 3,6\0,92|0,520,23!0,06 


2 Der Verlauf der Einflußfunktionen ist nahezu rotationssymmetrisch. 
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a a a 
u 


1 
300 ö 
Al 


600! 
700! 
sal 
900} 
1000" 1060 “900 0 300 d 30 600 900 1200 km 


a = ” - Te 2 
G (X45 X X Xp Xa, X) = G (3, Xe, X35 X], 09,85) 


Im Falle, daß die Randwerte u(s) verschwinden, ver- 
schwindet das Randintegral. 

Die Lösung des Differenzengleichungssystems läßt 
sich dann analog zu 


ES GG ht 
= D GT Ph 
darstellen und die Einflußgitterfunktion I(r, r’) wegen 


der Symmetrie unmittelbar als Greensche Funktion 
interpretieren. 
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tool = ae : zen RTE 002 ee. 
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ENS, ie 


was wegen 7 =f im Falle der gewählten 
Standardatmosphäre erfüllt ist. 

Im vorliegenden Fall verschwinden die Rand- 
werte z an den Grenzflächen p = 0 und p = 
p* = 1000 mb jedoch nicht notwendigerweise, 
so daß auch die Symmetrie der EinfluBfunktion 
geringfügig gestört ist. Analytisch kommt die- 
ser Sachverhalt in der unterschiedlichen Form 
der Koeffizienten des Differenzenoperators der 


Gl. (16) für Randnachbarpunkte und die übri- 
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é à F d ae 
Fig. 4—8. Einflußfunktion 7 [sal mit einem nichtver- 
schwindenen dynamischen Jacobiausdruck pu? JE} (z, 4) = 
d 
= 100 [sal im Störpunkt r’ in @ # 51° Breite und 
p’ = 975 mb (Fig. 4), 725 mb (Fig. 5), 475 mb (Fig. 6), 
275 mb (Fig. 7), 75 mb (Fig. 8). T,= —6°/km in der 
Troposphäre; 7 = f. 


gen inneren Gitterpunkte zum Ausdruck. Auch 
im Falle der Asymmetrie von I, was insbeson- 
dere für 7,+ o gilt, kann man jedoch stets 
wegen der Gültigkeit des Superpositionsprin- 
zips bei durch die « gegebenem Differenzen- 


-ausdruck eine lösende Funktion G aus den Ein- 


Außfunktionen sich so zusammensetzen, daß 
Den Gi) hir, 

Summation über alle inneren Gitterpunkte. 

Man findet GYif = Ii?’. Das bedeutet, daß vor 

Aufstellung der Greenschen Funktion oder bes- 


ser des lösenden Kerns G die Einflußfunktionen 
für sämtliche Innenpunkte zu berechnen sind. 
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1000 r 7- T T T T D 
1200 900 600 300 0 300 600 900 1200 Ikm] 
Fig. Io. 
mb 
0 s—— —— lu 


10001  — 
2 1200 900 


Fig. 9. 


! : 4 d 
Fig. 9—13. Einflußfunktionen 7 [sal mit einem nicht- 
: 3 gs dm 
verschwindenden Jacobiausdruck Pa len 
dj o 
entspricht einem Werr( eu —v-V r) = 4 al 
2 dt kr Std | 
im Störpunkt r’ in @ & 51° Breite und p’ = 950 mb 
(Fig 9), 700 mb (Fig. 10), soo mb (Fig. 11), 300 mb 
(Fig. 12), 100 mb (Fig. 13). T, = —6°/km in der 
Troposphäre, n = f. 


Jedoch wird ein derartiger lösender Gitter- 
kern G bei Behandlung aktueller Wettersitua- 
tionen, also bei beliebiger Anordnung der « 
und h-Werte, nur dann mit Vorteil zu benutzen 
sein, wenn sich erweist, daß auch eine Lösung 
mit einem von der speziellen Wetterlage un- 
abhängigen Differenzenausdruck (seitlicher 
Mittelwert) die tatsächliche Lösung approxi- 
miert. 

Man wird in diesem Falle diese Approxima- 
tion als erste Näherung im Relaxations- bzw. 
Operatorenverfahren zur Abkürzung des Lö- 
sungsprozesses mit dem aktuellen Differenzen- 
operator verwenden können. 
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Liegt eine Lésung t vor, so gewinnt man 
durch Differenzenprozesse das Feld der Verti- 
kalgeschwindigkeit aus (17) 


Op = 2: (Te+4— + vQR) 
VOB 


und das Tendenzfeld der Temperatur sowie 
die horizontale Divergenz aus den beiden fol- 
genden Betimmungsgleichungen: 


(T,)r = — Hr (Tex — Te-4) 


— 4 (Div Vk = Wr+4— Wr-4 


Für Quellpunkte in Tropopausennähe sind 
diese letzten 3 Felder für ein u?J(z, 7) = 100 
(m/Std.) in p’ = 275 mb bzw. p Q = 100 (m/ 
Std.) im Quellpunkt p’= 300 mb in den Fig. 
14—19 dargestellt worden. 

Die EinfluBfunktionen lassen unmittelbar die 
Flächen gleichen Einflußgewichts und die prak- 
tische Einflußreichweite einer Störung in Form 
einer Vorticity-Advektion bzw. eines pQ-Wer- 
tes im Quellpunkt ablesen. 

Sie vermitteln andererseits auch einen Ein- 
blick in den atmosphärischen Mechanismus, 
wie er durch die definierten Störquellen aus- 
gelöst wird. Insbesondere liefern die zugehöri- 
gen Lösungsfelder eine quantitative Vorstellung 
von dem atmosphärischen Kompensations- 
mechanismus. 

Der Annahme quasistatisch verlaufender Pro- 
zesse liegt zwar die physikalische Conception 
unendlicher Signalgeschwindigkeit zugrunde, 
d. h. es existiert ein starrer Übertragungsmecha- 
nismus sowohl zwischen horizontal als auch 
vertikal benachbarten Gitterpunkten, der eine 
zeitlose Beeinflußung durch eine Störquelle 
auf den meteorologischen Zustand eines beliebig 
weit davon entfernten Gitterpunktes zuläßt. 

So enthält z. B. die Tendenzgleichung in der 
aus der Kontinuitätsgleichung ®„= — Div v 
und oberen und unteren Randbedingung ab- 
leitbaren Form 


p 
I 
OR ee 0 fo eat 
+ 
80 : 


die physikalische Aussage, daß eine Massenan- 
häufung in einer Schicht dp eine nicht von 
deren Entfernung vom Erdboden abhängige 
Tendenz am Boden hervorruft. Wegen der 


integralen Abhängigkeit muß jedoch stets die 
Möglichkeit unterstellt werden, daß eine Diver- 
genz der einen Schicht in ihrer Wirkung auf 
den Bodendruck durch entsprechende Vor- 
gänge in Schichten darüber oder darunter kom- 
pensiert wird, wobei a priori nichts darüber 
ausgesagt werden kann, ob derartige Massen- 
divergenzen willkürlich über die Vertikalsäule 
verteilt sein können, oder ob funktionelle Be- 
ziehungen zwischen diesen Größen verschiede- 
ner Schichten bestehen. 

Die synoptische Erfahrung erweist jedenfalls, 
daß allgemein das Kompensationsprinzip so- 
wohl in der Form 


Lf Div v dp| < f |Divv|dp 
als auch in der Form 
ifT.dlnp|<f |T;|dInp gilt; 


was auf eine Korrelation zwischen den Inte- 
granden innerhalb der betrachteten Vertikal- 
säule schließen läßt. 

Beide Quantitaten — Div v, T,— sind den 
Topographien eines Termins nicht zu ent- 
nehmen, zumindest nicht vollständig. Ableit- 
bar aus simultanen aerologischen Meßdaten ist 
lediglich ein gewisser Beitrag, nämlich die 
ersten Terme der rechten Seiten folgender 
Bestimmungsgleichungen für Div v und T;: 


nDivv=—v-Vn—4; 


9 9 
A=n+ Op = (5+ os) 4az 


T,=—v-VT+B,; B, = (j-°5)T 


T, = 5 pQ+ Ba; By= 500 1 


also die horizontalen Transporte der Vorticity 

und der Temperatur, bzw. mit der einschrän- 
= ; ; LA 

kenden Annahme über die Kenntnis von 1 die 


dt 


t Man erhält diese Formulierung für T, aus dem 
1. Hauptsatz in der Form: 


OTe Ada ET dp 
dt Cy dt p dt 


=D -V-VT+ oT, Dann ist 


kIT gpQ P 
Ti 4 WIR 7: ers = 
t Sr ne W/ (2) L © : + ow 
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Fig. 14. 


Fig. 14—16. Zur EinfluBfunktion t der Fig. 7 zu- 
“IS 
gehöriges Temperaturänderungsfeld T, | (Fig. 14); 


b 
Vertikalgeschwindigkeitsfeld — « | (Fig. 15) und 


Feld der horizontalen Massendivergenz — y Div v 


ee 
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Fig. 15. 
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Fig. 17—19. Zur Einflußfunktion 7 der Fig. 12 zu- 
One x 
gehöriges Temperaturänderungsfeld T; [sal (Fig. 17) 


b 
Vertikalgeschwindigkeitsfeld — « | (Fig. 18) und 


Feld der horizontalen Massendivergenz — v Div v 
mb 


=] (Fig. 19). 
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Größe Q; während die weiteren Beiträge A 
bzw. B,, B, zum Divergenz- und Temperatur- 
tendenzausdruck nur durch die Rechnung be- 
stimmt werden können. 

Man könnte annehmen, daß die Kompensa- 
tion dadurch zustande kommt, daß 


IYv: Vndo|< flv: vnldp und dass 
p* DA 
lJ Qdp| < / |Qp| bzw. 

| fv:v Tdinp|< fiv-VTldbp (25) 


Jedoch dürfte eine solche Annahme physika- 
lisch nicht begründbar sein. Eine physikalische 
Erklärungsmöglichkeit für die Kompensation 
würde darin bestehen, daß gezeigt werden kann 
daß auch bei einer völlig beliebigen vertikalen 
Verteilung der Transportgrößen, also im un- 
günstigen Falle bei Gleichheit der rechten und 
linken Seiten obiger Ungleichungen, noch eine 
Kompensation stattfindet, also (25) oder zu- 
mindest 


| f Div vdpl< if v- vndp| und 
px px g px | 
|f Tdinp|< |/v:V Tdlnp|bzw. pif Qdp|gilt. 


Eine derartige Kompensation käme dadurch 
zustande, daß die Advektionsgrößen v-Vn, 
v-V Tbzw. J(z,n), Q, die wir als primäre 
Divergenzen bzw. Temperaturtendenzen be- 
zeichnen wollen, automatisch sekundäre Diver- 
genzen und lokale Temperaturveränderungen 
A, B,, B, im Quellpunkt selbst und in dessen 
Umgebung hervorrufen in einer Weise, daß 
sekundäre Felder A, B und primäre Quellen 
(J(z, n), Q) sich weitgehend kompensieren, ein 
Mechanismus, den man als automatische Kom- 
pensation bezeichnen könnte. 

Gerade diesen Sachverhalt bringen aber die 
Einflußfunktionen zum Ausdruck, wie man 
leicht den Meridionalschnitten der Fig. 14— 
19 sowie den Bilanzen innerhalb der den Quell- 
punkt enthaltenden Vertikalsäule entnimmt, 
die in ‚Fig. 20, 21 dargestellt sind. 

Man entnimmt diesen Bilanzdiagrammen 
unmittelbar, daß primäre Divergenzen und 
Temperaturtendenzen an einer bestimmten 
Stelle r’ des Raumes sekundäre Divergenz- 


bzw. Temperaturänderungseffekte am Quell- 
punkt r’ selbst sowie in dessen Nachbarschaft 
auslösen, die einmal den primären Vorgang am 
Quellpunkt selbst abzuschwächen suchen und 
sich mit den primären Effekten soweit kom- 
pensieren, daß die Lösungsfelder aller met. 
Größen mit zunehmendem horizontalen und 
vertikalen Abstand vom Quellpunkt rasch ab- 
klingen, wobei noch das Ausmaß der Kompen- 
sation wesentlich von den Stabilitätsverhältnis- 
sen in der Atmosphäre abhängt". 

Bei Abänderung der Stabilitätsverhältnisse o 
und n werden die Einfluß- und damit Green- 
schen Funktionen insofern modifiziert, als nied- 
rige 7-Werte einen stärkeren »Durchgriff« des 
Einflusses in horizontaler Richtung, niedrige 
o bzw. s-Werte (hohe statische Labilität) einen 
stärkeren »Durchgriff« in vertikaler Richtung 
auf Kosten der anderen Koordinaten bewirken, 
wobei die Gewichtsverteilung vorwiegend 


durch die Kombination 1 des «-Koeffizienten 


bestimmt wird. 

Als Beispiel zeigt Fig. 22 die zu Fig. 6 analoge 
Einflußfunktion, unter Zugrundelegung eines 
Temperaturgradienten von 9°/km anstelle von 
6°/km in der Troposphäre, die eine merkliche 
Bevorzugung der vertikalen Koordinate an- 
zeigt. Wie man sich diesen automatischen Kom- 
pensationsmechanismus physikalisch vorzustel- 
len hat, sei kurz qualitativ erläutert, wobei das 
zeitliche Nebeneinander des quasistatischen Vor- 
ganges notwendigerweise in ein zeitliches 
Nacheinander eines nicht quasistatischen Pro- 
zesses aufgelöst werden muß. Im Falle einer 


1 Wegen empirischer Untersuchungen über den Kopp- 
lungsmechanismus advektiver und konvektiver Tempera- 
turänderungen vgl. CraDDOCK (1951) und W. REGULA 
(1952). 


Fig. 20. Einflußfunktion der Höhentendenz 7 | sal: 
t 
: BI : mb 
der Vertikalgeschwindigkeit w =]. der Temperatur- 
FIT 
änderung ie [sal und der horizontalen Massendiver- 


; mb |. 
genz — y Div v Ea innerhalb der den Störpunkt 


enthaltenden Vertikalsäule, ausgelöst durch einen nicht- 
verschwindeden Jacobi-Ausdruck ul (z, n) = 100 


m ‘ : 
El im Störpunkt r’ mit @ & 51° Breite und p = 


= 75 mb (obere Reihe), 475 mb (mittlere Reihe), 
975 mb (untere Reihe). 
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; ; 3 nl 
Fig. 22. Einflußfunktion t Ea mit einem nichtver- 


900 1200 km 


schwindenden Jacobi-Ausdruck wu? Ji (2,7) = 100 
d 
Ea im Störpunkt r’ mit @ & 51° Breite und p’ = 


= 475 mb bei hoher Labilität in der Troposphäre: 
re = —9°/km, 7 =f. 


positiven Störquelle J(z, 7) > o (primäre Kon- 
vergenz) entsteht vorerst im Störpunkt selbst 
eine horizontale Massenkonvergenz, die wegen 
der geringen Kompressibilität der Luft aus 
Kontinuitätsgründen eine vertikale Massendi- 
vergenz erzeugt mit aufsteigender Luft ober- 
halb und mit absinkender Luft unterhalb des 
Quellpunktes. Unterhalb des Quellpunktni- 
veaus wird wegen des Massenzuflusses im Quell- 
punkt positive Drucktendenz erzeugt, oberhalb 
des Quellpunkts infolge der positiven Vertikal- 


> : _ 
komponente — > o ebenfalls eine positive 


dt 
Drucktendenz, d. h. die Niveauflächen heben 
sich in der gesamten Störpunktsvertikalsäule 
relativ zur Umgebung an. Diese Anhebung 


Fig. 21. EinfluBfunktion der Höhentendenz 7 | , der 
: VINS mb 
Vertikalgeschwindigkeit w I der Temperatur- 
ax 
änderung T, [=| und der horizontalen Massen- 
t 


divergenz — v Div v [Fal innerhalb der den Stör- 


punkt enthaltenden Vertikalsäule, ausgelöst durch einen 
nichtverschwindenden thermischen Ausdruck (p Oe = 
I dq 


m ; : 
= 100 | entspricht einem Wert =4 —v- 
G 
p 


Std 


Breite und p’ = 100 mb (obere Reihe), soo mb (mitt- 
lere Reihe), 950 mb (untere Reihe). 
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kr 
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bewirkt wegen der Trägheit des Windes, 
der nur mit einer gewissen Verzögerung 
dem ansteigenden Gradientfeld sich anpas- 
sen kann, in allen Niveaus der Quell- 
punktssäule wegen Überwiegens der Gradient- 
kraft über die Corioliskraft ein Ausströmen, 
d. h. Divergenzen, die sich mit der primären 
Konvergenz im Quellpunkt weitgehend kom- 
pensieren. Temperaturmäßig gleicht sich das 
infolge des Aufwindes entstehende Abkühlungs- 
gebiet oberhalb mit dem infolge des Absinkens 
erzeugten Erwärmungsgebiet unterhalb des 
Störniveaus (bei stabiler Schicht) in seiner Wir- 
kung auf den Bodendruck ebenfalls weitgehend 
aus, siche Fig. 7, 14—16. 

Im Falle einer positiven Q-Störquelle findet 
vorerst infolge der Expansion ein Aufsteigen 
der Luft mit konvektiver Abkühlung und po- 
sitiver Tendenz T>o oberhalb des Quellpunktes 
statt; das mit diesem Druckanstieg verbundene 
Massendefizit führt zu einem Druckfall unter- 
halb des Quellpunktniveaus und damit zu Kon- 
vergenzen mit aus Kontinuitätsgründen auf- 
steigenden Komponenten und konvektiver Ab- 
kühlung. In der gesamten Vertikalsäule findet 
demnach konvektive Abkühlung statt, die die 
primäre Erwärmung im Quellpunkt selbst weit- 
gehend kompensiert. Die ausgelösten Massen- 
divergenzen mit entgegengesetzten Vorzeichen 
unterhalb und oberhalb des Quellpunkts kom- 
pensieren sich gleichfalls weitgehend in ihrer 
Wirkung auf die Bodendrucktendenz (vel. Fig. 
12, 17—19).Die exakten quantitativen simulta- 
nen Vorgänge in der zugrundegelegten Stan- 
dardatmosphäre zeigen die Fig. 20 und 21 für die 
betrachtete Störpunktssäule mitp’ in Bodennähe, 
im mittleren Niveau und in der Stratosphäre. 
Soweit die der Berechnung zugrundeliegenden 
Gesetzmäßigkeiten (16) die atmosphärischen 
Vorgänge adäquat beschreiben, muß aus den 
Lösungsergebnissen geschlossen werden, daß in 
der Atmosphäre trotz der Annahme quasista- 
tischer Vorgänge und unendlicher Signalge- 
schwindigkeiten doch ein Nahewirkungsprin- 
zip gültig ist!. 

Das bedeutet, daß in praxi eine Quellpunkts- 
störung im oben definierten Sinne Deforma- 
tionsfelder der meteorologischen Feldgrößen 
auslöst, die räumlich rasch mit wachsender Di- 
stanz vom Quellpunkt abklingen und daß ande- 


! Vergleiche dazu die Betrachtung von CHARNEY (1949) 
über die Gruppen-Geschwindigkeiten mit analogen Er- 
gebnissen. 
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rerseits die Ursachen für eine beobachtete De- 
formation in einem bestimmten Punkte eben- 
falls in unmittelbarer Umgebung derselben 
ihren Sitz haben müssen. Man muß natürlich 
berücksichtigen, daß mit zunehmendem Ab- 
stand r vom betrachteten Aufpunkt auch die 
Anzahl der beitragenden Gitterpunkte mit r? 
anwächst, so daß trotz abnehmenden Einfluß- 
gewichts des Einzelpunktes der Einfluß weit 
entfernter Flächen den Einfluß der Nachbar- 
punkte überdecken kann. 

Jedoch ist zu bedenken, daß mit zunehmen- 
der Entfernung auch die Verteilung der inho- 


mogenen Terme in zunehmendem Maße auf 
den Flächen gleichen Einflußgewichts die Form 
willkürlicher oder zufälliger Belegungen an- 
nehmen werden, so daß sich die einzelnen Git- 
terpunkte der Einflußfläche weitgehend in ihrer 
Wirkung auf den Lösungswert des Aufpunktes 
gegenseitig auslöschen. 

Für die während des Herbstsemesters 1951 
empfangenen wertvollen Anregungen und 
Ratschläge für die Abfassung der Arbeit sei 
den Angehörigen des Instituts für Meteorologie, 
Stockholm, besonders Prof. Rossby, Dr. Elias- 
sen und Dr. Bolin herzlich gedankt. 
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A Power Law Representation of the Heat Transfer in the 


Lowest 100 metres of the Atmosphere 
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(Manuscript received March 16, 1953) 


Abstract 


On the assumptions that the flux of heat in the air is governed by an interchange coefficient 
which can be represented by a simple power law of the height above the surface, that the flux 
of heat in the ground is by molecular processes and that during the daylight hours the variation 
of temperature can be represented by a sine function, the phase angle of the diurnal temperature 
wave is calculated and is shewn to be in good agreement with observations. The phase angle 
of the annual temperature wave is also discussed. 

Calculations of Ky the coefficient of eddy conductivity of heat and of its mode of variation 
with height in unstable conditions suggest that it is approximately equal to Ky the coefficient 


of eddy viscosity. 


I. Introduction 


Despite the fact that the mean temperature 
distribution within the lowest 100 metres of 
the atmosphere has been described in con- 
siderable detail in a series of geophysical 
memoirs (JOHNSON, 1929; BEST, 1935; FLOWER, 
1936; JOHNSON and Heywoop, 1938; Best, 
KNIGHTING, PEDLOW and STORMONTH, 1950) 
the theoretical conclusions from these obser- 
vations have been conflicting, partly owing 
to the complexity of the problem and lack 
of any clear cut picture of the physical pro- 
cesses involved in the mechanism of heat 
transfer, and partly to the effects of the differ- 
ent topography of the stations. 

In most discussions of heat transfer in the 
atmosphere (cf. SUTTON, 1949), it has been 
assumed that molecular conduction is neg- 
ligible, except within the lowest cm or so 
above the surface, that radiation may be 
neglected and that, at least in the mid period 
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of the day, convection is mainly responsible 
for the transfer of heat in the vertical, although 
there is disagreement as to the roles played 
by ‘free convection’ due to buoyancy effects 
and ‘forced convection’ due to frictional tur- 
bulence. On the other hand, Riper and Rosin- 
son (1951) have stated that in the lowest 
two metres above the ground the absorption 
of long wave radiation by the air is at least 
comparable in effect with convection. 

In view, therefore, of the uncertainty as 
to the mechanism of heat transfer, Richardson’s 
direct approach, which avoids the necessity 
of making any assumptions as to the nature 
of that mechanism, has been adopted in this 
paper. 

RICHARDSON (1919) defined K the coefficient 
of eddy diffusivity for any transferable pro- 
perty x of the atmosphere by 


Pr _ 24 2 


1.00 
Dt 22 dz 
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where the only limitation imposed upon x is 
that y must represent some property which is 
unchanged by transfer to a different level. 

Whilst the coefficient K in this equation was 
originally called the coefficient of eddy diffu- 
sivity, this interpretation is not essential and 
if in heat transfer part of the temperature 
change at any level is due to factors other than 
eddy motion, equation 1.00 can still be used to 
represent the whole of the observed variation, 
although in this case K is a composite diffusion 
coefficient which takes into account all the 
causative factors. 

It is evident that K as defined in this man- 
ner may be a function of y and z and hence 
may vary from property to property if the 
transferable property under consideration in- 
fluences the motion. Comparison of values 
of K for various transferable properties com- 
puted by means of this equation thus affords 
some evidence as to whether, as has been 
suggested by various writers on turbulent 
motion from Reynolds onwards, the transfer 
of heat momentum and water vapour is 
affected by the same mechanism. 

In the case of controlled wind tunnel ex- 
periments Pasquitt (1942) found that the 
coefficients of eddy diffusivity and momen- 
tum were practically identical, but in a second 
paper (1949) dealing with the diffusion of 
water vapour in the lowest two metres above 
the surface he concluded that whilst under 
stable conditions there was equality between 
the coefficients of eddy diffusivity for water 
vapour and heat, in unstable conditions the 
latter was substantially greater. 

Riper and ROBINSON (1951) however 
questioned Pasquill’s conclusions, and by a 
further series of measurements within the 
lowest two metres above the surface shewed 
that the wind and temperature profiles in 
nearly all cases had the same form and establish- 
ed a strong probability that the coefficients of 
eddy viscosity, eddy diffusivity and eddy 
conductivity were approximately equal under 
all circumstances. 


2. The present position of Ky 


It has been noted by many writers that near 
the ground the diurnal variation of tempera- 
ture can be represented approximately by a 
sine wave of period 24 hours, whilst at great 
heights it is negligible. If therefore, horizontal 
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advection can be neglected the problem of 
heat transfer in the lowest atmosphere reduces 
to the problem of obtaining a solution of 
20 9 = 
a 
where © is the potential temperature of the 


air at a height z above the surface at a time t, 
subject to the boundary conditions 


2.00 


O=0, + a singt whenz— o 2.01 
0=0, when z— oo 2.02 
27 
whereye- = = 7.3 10e, 
1 24 X 60 X 60 23 


and a is a constant. 

The solution of 2.00 which fits these conditions 
with Ky=a constant was given by G. I. 
TayLor (1915), and is 


9=0,+asin (gt -bz) exp-bz 2.03 
q 
Are 
where b ae 


Thus on the assumption that Ky =a constant, 
the lag in the time of maximum temperature 
with height should be a linear function of the 
height above the surface. 

Best (1935) however, from a very careful 
set of measurements between 2.5 cm and 17 m 
above the surface, found that the lag varied 
very closely as 222 whilst calculations, which 
he made on the assumption that Ky could be 
treated as a constant over certain shallow layers 
of the atmosphere, indicated that Ky increased 
with height according to the power law 


Krr= const. 21° 2.04 


so that the assumption that Ky is a constant 
is clearly unacceptable. Subsequently however 
CowLiNG and WHITE (1941) showed that 
methods which assumed that Ky could be 
regarded as independent of height through 
even a shallow layer led to erroneous results 
and could give values of Ky which were not 
even correct as regards order of magnitude. 
Values of Ky, which they computed by a 
variety of methods, indicated that Ky increased 
very rapidly with height above the surface, 
the rate of increase being greater by day than 
by night, and it was shewn by the present 
writer (1948) that their values of Ky under 
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night inversion conditions could be fitted 
very closely by a simple power law Ky= 
const. 22/3, 

In the case of momentum transfer it has 
been shewn by many writers, including the 
present writer (1948), that Ky can be repre- 
sented by a power law of the form 


Ky=czirm 2.05 


where the index m always lies between o and ı 
À I 

and varies from a value of about — under 

adiabatic conditions in the daytime to a value 


I I à : à 
of between = and 3 under night inversion 


I 
conditions, so that under inversion conditions 


the modes of variation with height of Kur 
and Ky are approximately the same. Under 
unstable conditions however the evidence 
from Leafield (JOHNSON and Heywoop, 1938) 
and Rye (BEST and OTHERS, 1950) is conflicting. 

SUTTON (1948) from an examination of the 
Leafield data found that in the hours around 


a became 
dt 

approximately constant with height and time 
above 1.2 m, and by making certain assump- 
tions as to the amount of heat convected from 
the surface and assuming that the horizontal 
flux of heat in the mid hours of the day was 
negligible (cf. however section 4) concluded 
that the variation with height of the vertical 
flux of heat could be ignored below about 


noon on a clear summer day 


90 
100 m. By plotting log = against log z, 
Sutton found that over the limited range 
90 
from 7 m to 45 m oe could be closely repre- 


sented by const. 27175, so that on the assump- 
tion that the vertical flux of heat was constant 
over this height interval it followed that over 
the same range 


Ky =const. 217% 2.06 


Sutton therefore argued that as the direct 
application of the solution of 2.00 with 
Ky= const. where m was a small 
positive fraction < 1, to the case of heat transfer 
in the hottest hours of the day led to results 
which were not in agreement with observa- 
tions, then the mechanism responsible for 
heat transfer must differ from that for momen- 
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tum transfer and suggested that this was due 
to buoyancy forces. 

On the other hand KNIGHTING and Best 
(1950), in a discussion of the Rye data, whilst 


making the assumption that 29 was approxi- 


ot 
mately constant and =A in the mid hours 
of a hot summer day, made no assumptions 
as to the constancy with height of the vertical 
flux of heat, and by integrating 2.00 when 
Kyp= cz!" obtained 


Azitm Bem 


c(1 +m) fi cm 


207 


where À, B and C are constants. By fitting 
this equation to a series of observations of 
potential temperature at four heights between 
I.1 m and 106.7 m and eliminating A, B and 
C, they derivedthe following equation form, i.e. 


0,-0, 
9-9, 
9-09, 
and by solving this equation graphically, 
found a value of m= 0.45. It will be noted 
that the value of m obtained in this manner is 


Lam 1+m m m 
24 21 24 zy 
Zu = zus zur it 210 =O 2.08 


1 1 
Ze me = +m mo er ae 


independent of A=, providing that A is 


a constant 0, but that the values of c and 
B depend upon the values of m and A. 


Unlike the negative value of m derived 
by SUTTON (1948) from the Leafield data on 
the assumption that the first term on the right 
hand side of 2.07 could be neglected in com- 
parison with the second, this value of m 
lies within the possible range of values derived 
from considerations of momentum transfer, 
although differing considerably from that 
appropriate to an unstable lapse rate. An ex- 
amination of the Rye temperature data on 


0 
clear summer days however shews that A 


falls from about 4.5 x1074° F. sec-1 at 1030 
to about 1.0 x 1074 °F. sec”! at 1230, so that 
calculations of m from small differences of 
potential temperature at Rye based upon the 
assumption of constancy of m: must accord- 
ingly be suspect. 

When the value of m was recomputed from 
the Leafield summer data using the observa- 
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tions at the four heights 12.4 m, 30.5 m, 
57.4 m and 87.7 m, using the method sugges- 
ted by KNIGHTING and BEST (1950), it was 
found that the potential temperatures over the 
period 1000 G.M.T. to 1300 G.M.T. could be 


fitted very closely with a value of m= > 2 


ÿ 
With this value of m together with the value 
of A—=33*x10,1PFisec 21 found by SuTTON 
(1948), the values of c and B calculated from 
2.07 are 22 and -3.1 respectively, so that the 
Leafield observations vield 


Ky = 22 2/7 2.09 


For comparison the observed values of the 
differences of potential temperature over 
various height intervals on a clear summer day 
and values calculated from 


OO Er ere 12.4087} — 
22087 N 
3-1 1/7 1/7 
= — — 12 
22 ars eg 


are shewn in Table 1. 

The potential temperatures used are means 
from 1000 G.M.T. to 1300 G.M.T. so that 
each entry in the second row of Table I re- 
presents the mean of about 30 observations. 


Table I. Observed and computed differences of 
potential temperature. 


Height interval inm 0e M Ape ore7 
A | 7 
Observed differences in °F. 0.31| — 0.46| — 0.52 
Calculated differences in °F| — 0.30] — 0.46| — 0.52 


The expression for Ky derived by the present 
writer (1948) was 


Ky =e Vhs m 2 Tel 


where V, was the velocity measured at a 
standard height h above the surface, z, was a 
parameter characteristic of the degree of 
roughness of the surface and m was a thermal 
parameter. 

The mean value of V;, during the period 
1000 G.M.T. to 1300 G.M.T. on clear sum- 
mer days at Leafield at a height of 12.7 m 
was 4.2 m sec-! and whilst the value of 2, 
over the grass covered surface is not known, 
it is unlikely to be widely different from the 
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value of z)=2.6 cm. found by a variety 
of methods for the rolling grassland of Salis- 
bury Plain. With this value of z, and a value 
of m= ig which was derived from velocity 


profiles measured near the surface under 


unstable conditions, 
Kr 2842 © 


Thus whilst with small differences of potential 
temperature, such as those in Table I, the 
possibility of observational error cannot be 
excluded, it would seem that, contrary to 
Sutton’s conclusions, the Leafield observations 
in the hottest hours of a clear summer day 
and at heights greater than 1.2 m are not 
inconsistent with predictions made from the 
heat transfer equation when Ky = cz!~™ where 


2.12 


m has a value of approximately >; and 


where c is of the same order of magnitude 
as mzy2"V,h-™. 

KNIGHTING (1950), in a mathematical discus- 
sion of the heat transfer equation, shewed 
that if any experimental results were analysed 
by means of solutions obtained from 2.00 
and this led to values of m which were negative, 
then the solutions were not applicable, or 
alternatively, if the experiments resulted in 
values of m which were negative, no use 
having been made of any solutions, then 
equation 2.00 was invalid and could not be 
used to represent the equation. 

In this connection, it is perhaps not without 
interest to note that if Sutton’s negative 
value of m is used to calculate c from 2.07 
using the Leafield observations at 12.4 m, 
30.5 m, and 87.7 m, then c= —5 x10-2 so 
that Ky is negative which is physically im- 
possible. 

KOHLER (1932), on the assumptions that near 
the ground the diurnal variation of tempera- 
ture could be represented by a sine wave of 
period 24 hours and that Ky could be repre- 
sented by cz!-" where c and m were constants, 
obtained a solution of 2.00 in terms of Bessel 
Functions of an imaginary argument and from 
this solution derived formulae which could 
be used for comparison with observations 
near the ground. With the aid of these for- 
mula Köhler computed the variation of the 
amplitude with height for various values of c 
and m over the height range from the surface 
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to about 30 metres and by comparing Jonn- 
SON’S (1929) measurements at Porton of the 
mean amplitudes at heights of 1.2 m, 7.1 m 
and 17.1 m with curves derived from these 


I 
computations for m=-, found that a close 
fit was obtained with 

Ky=13.828? 


Best (1935) who employed formula 12 of 
Köhler’s paper in order to discuss the variations 
of the lag with height of the time of maximum 
temperature found that for large values of z 
the lag tended to a law 


I+m 
lag = const. z 2 210 


and by comparing this with his empirical 
result obtained a negative value of m— -0.6 
which is incompatible with values derived 
from the momentum theory. 

Later however, Best (Meteorological Maga- 
zine, November 1935) concluded that the 
assumption that formula ı2 could be used for 
heights as low as those considered in his 
paper (1935) was not justified and that the 
deduction that m= -0.6 was invalid. 

If therefore we exclude Sutton’s conclusions, 
which appear to be open to some criticism, 
it appears that the Leafield observations both 
at night-time and in the daytime are not 
inconsistent with predictions made from solu- 
tions of the heat transfer equation on the 
assumption that Ky=const. Ky. 

In the following sections, where the con- 
sequences of this assumption are examined in 
more detail, it is shewn that the predicted 
values are also in very close agreement with 
observations at Ismailia and at Kew. 


3. The diurnal temperature wave at the earth’s 
surface 


It was shewn by the present writer (1948) 
that Ky=a const. z1-™, where m is a function 


of varying from a value of about yoke 
OZ 2) 


a 


in the early hours of the morning to a value 
of about — or less at midday and the Leafield 


observations suggest that a similar expression 
holds for Ky. Such a consideration in the 
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case of heat transfer however renders equation 
2.00 quite unmanageable, but since the varia- 
tion of m during the daylight hours is slight 
the solution of 2.00 with m=a const. can be 
used to discuss the temperature changes during 
this period. We thus require a solution of 


20 9 | eA | 


ioe |e oz) pe 
where c is independent of t and z 
subject to 
O = 0, + a sin gt wh 
9 + 4 sin gt when z>o | or‘ 


0-0, 


when z— oo j 


where a is the amplitude of the diurnal tem- 
perature wave at the surface and 


a TEN ER 2028 
IT 
Writing 
ung 
ae 3.02 

and 

O=0, + (exp igt) - f(é) 3.03 
in 3.00 gives 

N Lake 
(bn at ae f° 3.04 


The solution of this equation was given by 
the present writer (1948), and substituting 
this in 3.03 and selecting the real part, the 
solution of 3.00 which fits the required bound- 
ary conditions is readily found to be 


O=0,+ap(£) sin gt-ay(E) cos qt 3.05 
where 
p(£) = 
ee! fae 
Ve +m) la(2 + m)(3 +2m)(4 + 3m) 4 
mf Ee | 
CII ba ee aii + of + 
Tea ¢ _ 
ee een 0 


3.06 
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and 
& | 
w(é) = en tt 
= 2 
es a Jet Se) ae 
+ Céltm- 
on & 
‘V+ 2m Bene 
3.07 


and C and D are functions of m which were 
tabulated loc. cit. for various values of m 


I 
and whose values for m=- are 1.1161 and 


0.2221 respectively. Theoretical expressions 
for C and D derived by A. F. Crossley are 
given in an appendix to this paper, and it 
can be seen from equations 6.24 and 6.25 
of the appendix that 


I(2-n) NIT 
= ZN) sane 
@ Ta rer n) cos —, | | 
3.0 
er) —(1—n) 
TG aa n) sin 
where 
ee ) 
ı+m 2 
so that 
——— I[(2-n) 
2 IDe = — y\—(1—n) 
VC2 + ne n) 3.10 
and 
= tan Br tan a 2 
C 2 D es 


Equation 3.05 may be written 


@=0,+ay(é) sin (gt —A) 3.12 


where 


tan ine. & (8) =Vp(6) +pE) 3.13 


and these functions are tabulated and utilized 
in the following section. It can be seen from 
3.12, 3.13, 3.06 and 3.07 that very near the 
surface the lag with height of the time of 
maximum temperature is given by 
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ir à 
A= ran = tan! 25 = const. &m+1 
% I- CEm+1 


approximately = const. 2” 3.14 whilst at great 
heights, from the asymptotic expansion cf. 
2.13 the lag is given by 


1+m 


lag = const. z 2 3:28 


1 
4 I a == 
Thus with ae the lag varies as 27 near 


4 

the surface and as z7 at great heights so 
that it is not altogether surprising that BesT’s 
(1935) comparison of Köhler’s asymptotic 
solution for great heights with observations 
near the surface gave misleading results. 

The time of maximum temperature at the 
surface can be obtained from a consideration 
of the fluxes of heat in the air and ground as 
follows. The convective flux of heat C(o) 
from the surface 


à 
Es OCZ ME 2106 


and hence from 3.02, 3.05, 3.06, 3.07 and 3.08 


q m+1 
C(o) = mezgcal | 3.17 
{C sin gt + D cos qt} 
m 
g is /C2 + D2 si 
= mezgca| = = = VC?+ D? sin (gt + à) 
Bolin 
where 
a= ance I 
C 3-19 
_ but from 1 
=tan 2 
G m+I 2 
and hence 
m 7 
a= — 3.20 
m+I 2 


Thus the convective flux of heat from the 
surface can be represented by a sine wave and 
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is a maximum at a time f before the time of 
maximum temperature given by 


INTEL 2A 6m 
ı- — = hours 
M+I 2 


2 21 
27 M+I 


and it will be noted that this phase difference 
is a function of m only and is independent 
of the wind speed, amplitude of the diurnal 
wave of temperature or nature of the surface. 
The flux of heat in the ground may be 
derived by means of a similar analysis, with 
k,, the coefficient of thermal diffusivity of the 
soil, replacing Ky=cz!-”, the coefficient of 
eddy diffusivity of the air. The corresponding 
solutions may be derived very simply by 
writing m=1 in 3.16 and replacing c, (the 
specific heat of the air) by c, (the specific 
heat of the soil), @ (the density of the air) by 
0, (the density of the soil) and c by k, (the 
coefficient of thermal diffusivity) whenever 
they occur in equations 3.17 et sequitur. 
Hence the ground flux of heat G(o) is given 


by 


G(o) = aerka( À ) WG HD, sin (gt + x) 


2k, 
5 22 
where from 3.10 and 3.20 
VC,?+ D,2= V2 and =, 222 


and thus 


G(0) = ¢0,k, "ag": sin (a + =) 3.24 


Thus the ground flux of heat can be represent- 
ed by a sine wave and is a maximum at a 
time t before the time of surface maximum 
temperature given by 


pat ge hows 25 
A 29T 


and again it can be seen that this phase differ- 


_ence is independent of the amplitude of the 


diurnal temperature wave, nature of soil etc. 
This result was first obtained by SCHMIDT 
(1925) who also, on the assumption of a 
constant coefficient of eddy diffusivity in 
the air, found that the ground flux of heat 
was in phase with the convective flux of heat. 
Now the heat balance equation on the 
assumption that there are no horizontal 
Tellus V (1953), 4 
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gradients of temperature and water vapour is 
C(o) + E(o) + R(o) + G(o) =o 


where C(o) is the convective flux of sensible 
heat at the surface, E(o) is the convective 
flux of latent heat, R(o) is the net radiative 
flux of heat and G(o) is the flux of heat into 
the ground. Fluxes directed away from the 
surface are considered positive. 

Over a desert surface it may be assumed that 
E(o) is zero, and hence 


220 


R(o)= - C(o) - G(o) 3.27 
and 
q Mile 
R(o) = - miypea| | VC?+ De. | 
= TT 1, ly - IT 
lerne hr €,0,k, "ag" sin | gt + 5 
3.28 
= À sin (gt +y) 3.29 
where 
2m | 
/ [ q m-+1 
A=a \/ ge] | (C2+ D?) + 
\ c(t +m) 
les | 
jm+1 ( 
+ C12012k19 + msg] | : | 
CR Dinner le PAL 
VC? + Dock, lg cos fe | | 
3.30 
and 
q Re 5 | 
mise | MED | 
tan y= — — ze 
Pa Be ae roam 
a 2 de € 5 
moe] + 5j \ Laser 
| 
| TL cok tha’ sin = 
= N) I 2 12 == 
sin ar 1e if 1 4 
a, Ok, ':q’' cos = 
MES 27 | 


Thus the net radiative flux of heat at the 
surface R(o) is also represented by a sine wave 
which is out of phase with the temperature 
wave. 
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The net radiative flux of heat on a clear 
day is however a maximum at local noon, 
and hence the time of maximum temperature 
at the surface occurs at a time 12+f hours 
where 


where y is given by 3.31. 
It can be seen from 3.31 that y must lie between 


m ZT TT \ 
.-and-and hence over the desert it 
M+I 2 


follows that the time of maximum must occur 


5 om 
not earlier than ae hours and not later than 
m+I 


3 hours after local noon, the actual time de- 
pending upon the relative magnitudes of the 
convective and ground fluxes of heat, but 
being independent of the amplitude of the 
diurnal wave of temperature. With Ky const. 
with height, m=1 so that the time of maxi- 
mum would always occur 3 hours after local 
noon. If Kyi«Ky, so that cxmz,2”V,h-™, 
then an increase in wind speed, or to a lesser 
extent an increase of the surface roughness, 
should result in an earlier time of maximum. 

Over a grass covered surface from which 
evaporation is taking place the effect of the 
flux of water into the air is difficult to assess 
directly. According to Bowen (1926) how- 
ever, the ratio between the flux of latent heat 
and sensible heat is constant for any given 
surface and may be obtained from simul- 
tancous measurements of humidity and tem- 
perature at two heights, and although it is 
not possible to prove this theoretically, RIDER 
and ROBINSON (1951), who calculated values 
of Bowen’s ratio 4 near the surface on clear 
summer days at Kew, found that # was 
constant with height and shewed little varia- 
tion from day to day. If therefore the fluxes 
of sensible and latent heat are in phase, equa- 
tions 3.27 to 3.31 can still be used with (1 +4) 
C(o) replacing C(o) over a grass surface. 
Hence if the other variables remain unaltered, 
the time of maximum will be earlier over a 
grass surface from which evaporation is taking 
place than over a dry surface, but it will still 


ue om 
occur between the limits of SUR and 3 hours 
m+1 


after local noon. 
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4. Comparison with observations 


As this theory depends upon the funda- 
mental assumption that there is no horizontal 
flux of heat, it does not appear possible to 
use the results from either Rye or Leafield 
other than qualitatively. 


The station at Rye for example, is only 5 km 
from the sea and is affected by sea breezes 
especially on clear summer days, and the 
irregularities in the mean temperature curve 
for such days make it impossible to determine 
the time of maximum with any degree of 
confidence. Furthermore an examination of 
the individual daily records, from which the 
mean curve is constructed, shews that on 
many occasions in the hours around midday 
the temperature instead of rising is falling 
at all heights, so that any theory which 
depends upon the simple propagation of a 
sine wave upwards from the surface is clearly 
incapable of explaining the observed vertical 
thermal structure at Rye on such occasions. 


Leafield on the other hand, although free 
of complications such as a sea breeze, is 
situated on a spur of the Cotswold Hills 
with a rapid fall of 100 m in 3 km to the 
north and a fall of the same amount in about 
4 km to the south, and the meteorological 
observations at Leafield have been found 
to be unrepresentative of conditions over 
more level sites. JOHNSON and HEYwooD 
(1938) for example, list several anomalies 
in the meteorological elements at Leafield on 
clear summer days and ascribe them to the 
fact that during the hottest hours of the day, 
the surface air at Leafield mainly consists of 
air which has come from the relatively great 
height of about too m above the lower 
ground to the north and south of the spur. 


The desert site at Ismailia however is singu- 
larly free from topographical defects except 
in the summer months when it is sometimes 
affected by a slight ‘sea breeze’ from Lake 
Timsah 3 km to the south east, and hence the 
observations from here, and mainly those on 
clear December days (FLOWER, 1936) have 
been utilised to test out the theory. 


Now from 3.18 the maximum value of 


: TA 
C(o) with m== is 
yf 
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RS TO 0) LEE 
(25) a 


= ac? x 1.520 x 1075 g cal. cm?sec—1 4.00 


where a is the amplitude of the diurnal tem- 
perature wave at the surface, and from 3.21 
C(o) is a maximum 45 minutes before the 
time of maximum at the surface. 

Similarly, using the values of «1, o, and k, for 
Ismailia given by FLOWER (1936) i.e., c1=0.19, 
01=2.63 and k,=10-? the maximum value of 
G(o) is from 3.22 


G(0) = 0.19 x 2.63 x 0.1 x ax (7.3 x 10-812 


ea 4.05 X10-* cal. cm? sec.~1, 4.01 


and hence from 3.30 the maximum value of 


R(o) is 


R(o) =a{ (4.0)? + cl#(0.1 520)? + 2 x 4.05 x 


x 0.1520 x c7/8 cos =] 4.02 
16 
whilst from 3.31 
ET ler 
cl8 x 0.1520 sin 7 + 4.05 sin — 
tany = sers 4.03 


IT TT 
cl8 x 0.1520 COS — + 4.05 cos — 
16 4 


Whilst at Ismailia there appear to be no direct 
observations of R(o), C(o) or G(o) which 
can be used to calculate c and a, R(o) and 
G(o) can be estimated with a reasonable 
degree of confidence. 

Thus on a clear December day the maxi- 
mum downward flux of short wave radiation 
at Ismailia may be taken as approximately 
140 x 10-4 gram cal. cm ? sec- ! and the albedo 
of the desert surface may be taken as approxi- 
mately 20 %. The former value is based on 
American observations in California and is 
supported by observations at Kew Observa- 
tory made when the solar elevation is the 
same as that at Ismailia, whilst the latter is 
based upon mesasurements by Ängström and 
observations from American aircraft at high 
levels over Death Valley, California (1951). 
The incoming long wave radiation of at- 
Tellus V (1953), 4 
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mospheric origin, using Flower’s data of 
temperature and vapour pressure, is about 
75 gram cal.cm-? sec-! whilst the outward 
long wave radiation from a desert surface 
at a temperature of 300° A is about 112 x 10-4 
gramcal.cm ?sec-!, Hence the net radiative 
flux of heat, which Dr G. D. Robinson of 
Kew Observatory, to whom I am indebted 
for this estimate, considers may possibly be 
in error by as much as 10%, is 75 x 10-* gram 
cal.cm-?sec-1. 

Dr Robinson has also drawn my attention to 
unpublished measurements of the ground flux 
of heat in arid soil in Palestine by Ashbel 
which give G(o) = so x 10"? gramcal.cm-?sec-1 
and by W. Swinbank in semi-desert soil 
near Melbourne which also give values of 
G(o) of about 5ox 10* gramcal.cm-?sec.”1. 

With a value of G(o)=s50x10- gram 
cal.cm?sec-1, the amplitude a of the diurnal 
wave at the surface is, from 4.01, 12.4° C. 
With this value of a=12.4° C and a value of 
R(o) =75 x 1074 gram cal.cm~-*sec-! the value 
of c from 4.02=21.9, whilst from 4.00 with 
C=2 ro and 12,4 © 2 C(0)= 284107: 
gramcal.cm- sec". 

Whilst there is some uncertainty as to the 
value of C(o) at midday, the above value is 
in accordance with the usual concept of its 
value, and it is of interest to note that the 
value of c is very nearly the same as that derived 
from the Leafield observations in the hottest 
hours around noon on a clear summer day 


with similar wind speeds. 
With this value of c=21.9 from 4.03, 


rg at 
2.27 sin 7 + 4.05 sin — 
I 


tan y= = tail 33.17. 


TT 
a und OSICOSES 
2.27 COST 4.05 C 7 


4.04 


and hence the time of maximum temperature 
at the surface at Ismailia on clear December 
1724800 
days should occur 33 — x — 
y 3360 360 
after the time of local noon, i.e. at 1413 local 
time. 


The values of x(&) and À from 3.06, 3.07, 


= 133 minutes 


and 3.13 with m=- for selected values of 


E are shewn in Table II. 
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I 
Table II. Values of y(£) and À when er 


en 710252117022 Iron 2702215271022 5.1022 
HAE) 73710-650) 5337| 391 .272| .224 
Ay 1325371 5 53 | 9°324 17.4221. 29/07/1003 52551 


The lag in the time of maximum temperature 
LR: À . ÀxX24 x 60 

in minutes from equation 3.12 is Bares 4A, 
and hence from 3.02 using the value of c= 
21.9, Table II shewing the amplitude and 


lag at various heights can be constructed. 


Table III. Percentage amplitude and lag in minutes 
at various heights. 


Height 3.07 | 22.04] 1.65 | 12.40] 50.78] 92.94 
8 emal.cm m m m m 
Percentage 
Amplitude] 73.7| 65.0] 53.7] 39.1| 27.2] 22.4 
Lag in min- 
tes... LO 24 38 Fas Ro E77: 


Table IV shews the time of maximum tem- 
perature at four heights above the surface as 
determined by Flower from the December 
clear day observations and those derived from 
Table III together with the time of maximum 
temperature at the surface i.e. 1413 hours, 
and it can be seen that the agreement between 


the two sets of times is remarkably close. 


Table IV. Times of maximum temperature at 
various heights. 


Height in m | TAT 16.2 | 46.4 | 61.0 
Observed times...| 1447| 1524| 1605) 1623 
Calculated times..| 1447| 1530| 1605] 1616 


It is of interest to note that in view of the 
doubtful accuracy of determining times of 
maximum temperature at various heights 
from smooth curves drawn through mean 
hourly values, FLOWER (1936) carried out a 
harmonic analysis of the variation of tem- 
perature at cach height and determined the 
time of maximum from the phase angle of the 
diurnal harmonic. 

Flower also calculated the times of occur- 
rence of the maximum temperature at each 
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of the four heights for each month of the 
year and, noting that the time of maximum 
especially at a feat of 1.1 m varied very 
closely with the equation of time, suggested 
that the variation of temperature at the surface 
was equal to the equation of time. It will 
be seen that provided the variations in the 
mean monthly fluxes of convective heat are 
small, this follows from the fact that the radia- 
tive flux of heat is a maximum at local noon. 

If, as might be expected by analogy with 
momentum transfer, c is proportional to the 
mean wind speed measured at some standard 
height, then the effect of the wind speed upon 
the time of maximum temperature can readily 
be determined. The value of c in 4.03 was 
determined at a time when the mean wind 
speed at a height of 15.2 m was 3.8 m sec-! 
and hence if the other factors remain un- 
changed the values of y over a desert surface 
are, from 4.04, given by 


( =) Se UT 
2.27| —— sin — + 4.05 sin — 
3.8 4 


Fi 16 
pt VAE A 
2.27| — — + 4.0 = 
(7) cos = 4 ss ie 


so that an increase in the wind speed results 
in an earlier time’ of maximum. 

The values of y and the corresponding times 
of maximum temperature for various values 


of V, are given in Table V. 


Table V. Variation with wind speed of time of 
maximum temperature over the desert. 


Wind speed 

at 15.2 min o 2%, 5 7% 10 
al Sees 

DET EST: 452 |136%797|737°212| 23°32226270% 


Time of maxi- 
mum at sur- 


1425| 1406) 1354| 1345 


mum at I.I 


1504| 1439| 1425| 1414 


Time of maxi- 
mum at 16.2 


1556| 1519| 1459| 1445 


Time of maxi- 

mum at 46.4 
1635| 1548| 1523] 1507 
Time of maxi- 


mum at 61.0 
1558| 1532| 1515 
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There are unfortunately no observations in 
the Ismailia report (1936) which can be used 
to check this dependence of the time of maxi- 
mum upon the wind speed, but that at least 
the order of magnitude of the effect is correct 
can be seen from tables given by JOHNSON and 
Heywoop (1938) shewing the times of maxi- 
mum at various heights at Leafield in summer 
for various wind speeds. 


Table VI. 
Wind speed at 12.7 m in 
mn, GES = 47 75 
Time of maximumati.2 m| 1505} 1455| 1430 
ee ” ” 12.4m| 1546| 1535| 1500 
Be = ” 30.5m| 1600| 1550] 1512 
BE: ” 57.4m| 1616| 1606| 1522] 


A curious feature of the Leafield observations, 
which has no counterpart in either the Rye 
or Ismailia observations, is the fact that on 
clear summer days at Leafield the times of 
maximum temperature at all heights occur 
very nearly an hour later than the correspond- 
ing times of maximum on a mean June day, 
and in fact occur outside the theoretical times 
of maximum based on the assumption of 
very light winds and a very high value of the 
thermal conductivity. This is consistent with 
other anomalous features of the Leafield obser- 
vations on clear summer days noted by 
JoHNsoN and Heywoop (1938) and lends 
support to their suggestion that these were 
due to advection and lateral mixing. 

As already stated, the value of V, at Is- 
mailia at a height of 15.2 m during the midday 
hours was 3.8 m, and assuming that the 
parameter z, for Ismailia is of the same order 
of magnitude as that for the rolling down- 


land of Salisbury Plain, then Ky with m= 


is given by 
I 380 
je 2/7 __ 2°? cn 
a =e ) (1520)? 4.06 
= 25.6287 
I 
whereas with m or 
Kp= 21.927 4.07 
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which also lends some support to the suggestion 
that Ky=const. Ky, where the constant is 
approximately unity. 

From 3.05, 3.06, 3.07 and 3.08, it can be seen 
that around midday and within a few metres 
of the ground © is given to a good degree of 
approximation by 


m 


O=0,+a sin gt —a&™+1(C sin gt+ D cos qt) 


4.08 
so that from 3.02 and 3.11 
m | 
m+1 
O=0, + a sin qt az" 1 ] . | 
(m ip 1)c 4.09 
- VC? + D? sin (a+ 2 2) 
m+i 2 | 
and hence near the surface 
0=0, - Bz™ 4.10 


where ©, is the value of © at the surface and 


m 


D= | vor sin (gt = 2) 


m+1 +I 2 
ARLE 


It follows therefore that if Ky and Ky, are 
represented by similar power laws then the 
temperature and velocity profiles near the 
ground should be the same, and conversely 
that if the temperature and velocity profiles 
near the ground are the same then Ky and 
Ky are represented by similar laws. 

Now Riper and ROBINSON (1951), from 
observations at Kew Observatory, shewed 
that when the winds and temperatures meas- 
sured at several heights between 10 cm and 
200 cm were superposed with adjustment for 
scale, so that the profiles as drawn by hand 
coincided at 25 cm and 100 cm, the wind 
and velocity profiles in nearly all cases had 
the same form. 

Table VII shews the velocities from 25 cm 
to 200 cm for all unstable cases which they 
listed for clear summer days at Kew 1949 
at times between 1100 and 1400 G.M.T. 
reduced to a standard velocity of 200 cmsec-1 
at a height of 2.02 m. The wind velocities 
are mean values over a 30 minute period 
which the authors considered was sufficiently 
long to give smoothymean profiles and suffi- 
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Table VII. Velocities near the surface at Kew on clear sunny days. 
Case No. I 2 3 4 5 6 7 8 9 10 1/oth law 
a Io 
Date 18/6 | 20/6 | 20/6 | 21/6 | 21/6 | 23/6 | 23/6 | 24/6 | 1/7 27 200[ >) 
Height \ Time | 17217. 1721, 114102177702. 7400, 111032 |.04322|074795 10171222 217728 = 
IRC 
25 159 161 158 152 150 155 153 163 158 156 159 
37 162 165 166 163 166 165 163 173 148 161 166 
58 173 177 171 174 174 173 175 179 168 173 170 
76 176 180 179 178 177 179 180 176 171 175 179 
102 187 187 189 188 185 188 186 189 184 185 185 
202 200 200 200 200 200 200 200 200 200 200 200 


ciently short to exclude the effects of the 
diurnal variation. For comparison the veloc- 
ities computed from the simple power law 


A EN Pa ote 
V=200 ( =) which gives a closer fit than a 


„ch law are also tabulated. 


It can be scen that over this height range the 
observed velocities can be represented rea- 
sonably closely by a simple power law where 


; I : 
the index m=-, and hence since the ve ocity 


and temperature profiles when adjusted for 
scale coincide, it follows that over this height 
range the temperature data can be represented 
by 

0 =Og - Bz1/9 


where © is the value of © at the surface, and 
where 


Hes 


ß = V202 à Oo, = O202 
(202)"® Voge — Vos 
The values of 8 in °F. cm-¥/ for the above 


occasions have been computed and are listed 
in Table VIII. 


4.13 


Table VIII. Values of f on clear summer days. 


20 
C(o) = Lt. - Ku0% ay 


PS E x age 
800 5-5 00 


I 


I 


1.06 x 1072 xc gram cal..cm see 
Ars 


The mean value of C(o) for the occasions listed 
in Table VII and computed by Rider and 
Robinson from the heat balance equation is 
C(Q)=37.3 x 10 ~* gram'cal. cm 2 secure 


and hence from 4.15 and 4.16 


CES 
This value of c is subject to some uncertainty, 
as the computations of C(o) from the heat 
balance equation depend upon the funda- 
mental assumption that there is no horizontal 
flux of heat, and in view of the topography 
at Kew it is possible that this does not hold. 
Thus Riper and ROBINSON (1951) noted that 
the temperature fluctuations which they 
observed near the ground appeared to be 
inconsistent with this assumption, whilst 
values of C(o) computed by LANDER and 
ROBINSON (1952) from simultaneous observa- 


| tions of wind and temperature fluctuations 
| Case No. | 1 2|3|#|5|s|7|s|s |» + I 

| | | | | | | were consistently only about = to — the 

B....... 6.315-815.31 5.61 5.01 6.11 4-2] 5-815-914-91 values computed from the heat balance 


The mean value of f from Table VIII is 5.5, 
and hence near the surface where © is given 
by 4.12 


Ky=cz8!9 4.14 


The convective flux of heat at the surface is 
thus 


equation, so that the value of c lies somewhere 
between about 12 and 35. 

Now the mean value of V,, at 2.02 m for 
these observations was about 3 m. sec-! and 


: Ry 
with a value of m=- and a value of z= 


2.6 cm, the value of Ky=22.929/9. Whether 
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however this value of z, is approximately true 
for Kew is uncertain, for whilst the actual 
observations were made over a. surface of 
level grass which was kept mown to an 
average length of 1 to 2 cm, this grass surface 
was surrounded by a solid fence about 1.5 m 
high! which was about so m away at the 
nearest point, whilst in some directions 
there were very large trees and buildings, 
including the observatory itself which was 
on rising ground to the north west of the 
site. In view of the uncertainty of both c and 
Zo, all that can be said with any degree of 
confidence is that under unstable conditions 
the modes of variation of Ky and Ky at Kew 


à aE R 
are the same, with m=— approximately, 


whilst the observed values of the convective 
heat flux are not inconsistent with those pre- 
dicted on the assumption that Ky=const. Ky, 
where the constant of proportionality is 
approximately unity. 

The general agreement between the observ- 
ed and calculated values at Leafield, Ismailia 
and Kew suggest that, except perhaps in the 
lowest few cm, the mechanism for heat 
transfer is the same as that for momentum 
transfer and that the coefficients of eddy 
viscosity and eddy conductivity are approxi- 
mately the same. 


5. The annual temperature wave 


The results of the previous sections can be 
used with slight modification to discuss the 
annual variation of temperature. For example 
during the summer months at Ismailia the 
annual variation of temperature can be closely 
represented by a single sine wave of period 
one year. Similarly during the period May to 
September the average monthly lapse rate 
of potential temperature between 1.1 m and 
61 m is positive, 1.€. the potential temperature 
decreases with height so that during this 
period the index m can be assumed to be 
approximately constant and equal to about 1/,. 

The solution of 3.00 with m=1/, will there- 
fore be utilized to discuss the variation in the 
date of occurrence of the maximum air 
temperature at Ismailia. 

The phase lag of the annual temperature 
wave against the annual convective flux wave 
is from 3.20 given by 
Tellus V (1953), 4 
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[0 A 


m TT 
2 


5 .00 
INT 3 


and since for annual values =, the con- 
305 


vective flux of heat from the surface is a maxi- 
mum at a time f, before the time of maximum 
temperature at the surface given by 


ie 12 | 


= : months 5 


= 111/, days approximately | 


The phase lag of the annual temperature wave 
against the annual ground flux wave is from 
3.23 given by 


and hence 


5.02 


= 11}, months 


and it will be noted that, as in the case of the 
diurnal wave, these times are independent of 
the amplitude of the annual temperature wave 
and of the physical properties of the air and 
soil, etc. 

The phase lag of the annual temperature 
Wave against the net radiative flux wave is 
given by 3.31 with q, replacing q. Thus 
using the same values of the soil constants at 
Ismailia as in 4.01 and the same value of c, 
the time interval from 4.03 is 


5.03 
= 33. days 


Now from 3.12, the lag in days at various 
Ade 2056. : ; 
heights is ay ae approximately, whilst 


from 3.02, assuming that c does not increase 
with the increased time over which the means 
are now taken, the heights corresponding 
to the specified values of & in Table III have to 
be multiplied by (365)7*. 


The dates of maximum air temperature 
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Table IX. Dates of maximum air temperature at Ismailia at various heights. 
Height in m | TAT 16.2 | 46.4 | 61.0 
Observed date of maximum. a. eu see July roth July 29th July 30th July 31st 
Calculated date of maximum............... July 26th July 28/29 | July 30th July 30/31 


from Flower’s (1936) harmonic analysis of 
the monthly mean temperatures at 4 heights, 
during the period October, 1931 to October, 
1932, together with those derived by the 
present theory for the same heights, are given 
in Table IX. 

The agreement between the observed and 
predicted dates is very good except at 1.1 m 
where Flower’s date of maximum, in the 
light of his values at greater heights, appears 
to be somewhat early, and suggests that c is 
independent of the period over which means 
are taken. 

It is possible that this discrepancy is due to 
the limited period over which means are 
taken, or it may be that the present theory 
has oversimplified the problem of the annual 
temperature wave. 

Acknowledgment: I am indebted to the 
Director of the Meteorological Office for 
permission to publish this paper. 


6. Appendix by A. F. Crossley, M. A. 


Equation 13 (Frost, 1948) c.f. also 3.04 of 
the present paper, is 


2W IW 
er + FE 


iW=o 6.00 


Substituting 
m in 


W = yé? and x = 2(1 - n)!2e 4 E12 +. 


where 
m 
= 6.02 
Tat 
gives 
ayia dy ie 
= = Q 
EE REE (x =)! 2 Ge 
which is Bessel’s equation of order n 
Since o<m<sI 6.04 
therefore o<n <}/, 6.05 


so that (except for n=o) y=J,(x) and y=J_„(x) 


are independent solutions of 6.03 and hence 


= CAES (x), and W,=é 27, (x) 6.06 


are independent solutions of 6.00. 


Now 
wie Ke 
Ja) = 2" (1 + Alt | 2%{r1+n) i 
a 6.07 
di 24/2(1 + n)(2 +n) a 
and hence the two solutions of 6.00 are 
(1 —n)—"!2¢ 
— En] NEE 
ee. 
1262 
à +iE+ EG Em + 
and - 
(1 ne 4 | 
= En] (x) = 
tl ST aa 
ig i2€2 | 
sie "Irzm |2(1 + 2m)(2 + 3m) FA -| 


Except for the constant factors these solutions 
coincide with equations 17 and 18 (Frost 1948). 
We require asymptotic expressions for the 
two series which (following Frost) we sepa- 
rate into their real and imaginary parts. Thus 


ins 


AE +if© =P (x -ne * (a — meer] (x) 
6.10 
and | 
KO + (9) = Ts + met (x = ment, (x) 
6.11 
Now from 6.01 amp x= -”, and the corre- 


sponding asymptotic expression for J,(x) is 
(Gray, MATTHEwS and MACROBERT, Bessel 
Functions, p. 61) 
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2, 1/2 
In(x)> (2) cos (x-7 =) 6.12 
in 1 
e8 E « 
D RE sue = ] . 
ER 6.13 


+ COS {a — n)Y2(ı - i)EU2 — vn a 
4 


2 
or putting 
in = and u={2(1-n)Ë}12 6.14 
LOS 
ese 4 
JG) > Tr — n)14 : 6.15 
+ {cos (u — a)cos hu +i sin (u — &) sin hu} 
ce 
PC Pane ARE PE 
Te ar Ze À 2) 6.16 
2° (1 —n)* 
Similarly 
1 
~ 4 ou fe ie 
ea B 2) 6.27 
2907 (1 —n)* 
Therefore 
f (6) + ifp()> - 
Az) 271 ty a nn 
OS i nn) 0,75 if st) 
a 
27° 
6.18 
and 
fs(&) + ifs(€)> - 
(1+2n) (2n—1) i(u a) 
SOE a ae 
rt 
TR 
6.19 


Equations 26 and 27 (Frost, 1948) therefore 


become 
__(1—2n) 


Dee BG, (Te n)(1= nu # 


(1+2n) 


+ COS (« = . + Pr +n)(1 - n) 4 


8 
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and 
(1—2n) 


0=-9GFlı-niı-n * 


__(1+2n) 
sin (u 94) + rt + m) ; 


6-27 
whence 
C cos („-3-%) — D sin (5-7 = 
& 4 6 A 
ee) D. A 
0 Gr nt eo nn 
Co (x nu n)" cos (« 5 =) 
6.22 
C sin („-5-®) + D cos (en) = 
fh Sy. fll 
(1 -n) En T nm 
= Cor, :: a —n)" sin (« à +) 
6.23 
whence 
I(i-n) 
C= Cor, - ae —n)"cos— 6.24 
and 
I(1-n) 
D=0 Tat mie —n)t sin 6.25 
These forms are suitable for computation. 
The figures in Table II (Frost 1948) are con- 


firmed (apart from some minor differences 
in the fourth place of decimals). 


Also 
D NIT 
tan Wo ir Mm 6.26 
nt m m 
u er, ce 
which confirms equation 30 (Frost 1948). 
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A Remark on the Energy Transport from a Warm River 
Surface into Cold Air 


By A. NYBERG and L. RAAB, Swedish Meteorological and Hydrological Institute 


(Manuscript received September 10, 1953) 


Abstract 


Temperature measurements in very cold air above an open river surface indicate that in 
spite of very high temperature lapse rate (temperature of river, in the 2 cm, Io cm, and 
12 m levels are + 0.5, — 14.0, — 15.3, and — 18.0° C respectively), the eddy conduc- 
tivity is of the same order of magnitude as in extremely stable conditions (e. g. above a 


snow surface during nocturnal cooling). 


The vertical transfer of energy in the at- 
mosphere is due to vertical variations in the 
mean kinetic energy (vertical velocity shear 
caused by friction or other effects) or due to 
vertical statical instability (thermal convection). 


The transport through a horizontal surface of 
99 


dz 
Cp is the heat capacity of air at constant pressure, 


unit area can be expressed as — c,A,—— where 


A, the eddy conductivity and = the lapse 


rate of potential temperature. The eddy con- 
ductivity is a function of the lapse rate. In 
stable conditions the dynamical turbulence is 
damped whereas in unstable conditions there 
is a tendency for convection cells to form. In 
unstable conditions the thermal convection is 
generally supposed to be the dominating 
factor. When in an air layer the density 
increases with height auto-convection is ex- 
pected to set in, i.e. the colder upper air should 
break through to the bottom of the layer and 
the warm air rise to the top of the layer. 
This condition is reached when the lapse rate 
amounts to 3.4°/Ioo m. 

Some measurements, recently made at the 
river Indalsälven near Östersund in central 
Sweden, seem to indicate, that considerable 


Tellus V (1953), 4 


thermal convection does not start as soon as 
the lapse rate 3.4° in 100 m is reached. 

Fig. 1 shows the location of measurements 
and the estimated prevailing wind. The 
surface of the river is not frozen due to the 


Direction ot 
water flow 


100 200 300m 
1 iL 1 


un 


Fig. 1. Observational site. 
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regulation of the water flow system for the 
water power industry. The water.from the 
lake Storsjön which is covered by a sheet of 
ice almost half a meter thick has a temperature 
of about +0.5° C when leaving the lake. The 
flow in the river is so rapid that no ice cover 
is formed until more than 2 km from the lake 
even at air temperatures of — 20°C. 

The measurements on 12 February 1953 
were made at the point A on land 100 m on 
the windward side of the river and at the point 
B on the bridge over the river. 


al Over open river Over snow tield 


Over open river 


20 


-20 


Temperature height curves in the two cases 
refered to in the text. 


Bien 2. 


Fig. 2 shows the observed temperature 
profile. (In fig. 2 is also shown the temperature 
curve obtained at an extremely stable situa- 
tion, see NYBERG (1939). The values are given 
also in Table 1 which in addition shows some 
values obtained on 11 February at a point C 
over the river 2 m from the shore. See fig. 1. 
The instrument used here was a direct in- 
dicating resistance thermometer with three 
thermistors (Western Electric 14 A) coupled 
to a standard ohmmeter circuit. The ther- 
mistors were mounted equally spaced on a 6 
m bamboostick. The measurements on the 6 
m high bridge (point B in Fig. 1) were made 
by directing the bamboostick down to the 
water surface and touching the water surface 
with the top thermistor. The bamboostick 
was then successively raised and the tempera- 
ture variation followed. The height above the 
water surface was measured with an estimated 
absolute accuracy of s mm. The accuracy of 
the height difference was about 2 mm. 
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From these measurements approximate val- 
ues of the vertical eddy conductivity may be 
derived under certain assumptions. 

Consider a horisontal surface of unit area 
in the level z. The vertical energy transport 
through the surface due to eddy motion is in 
unit time 

VE 

Let this surface move with the horizontal 
wind u; in the x-direction during the time #. 

The total vertical eddy transport of energy 
through the area is then 

ee 
- [64% dt 
0 


Assuming at first that the horizontal wind 
u is independent of z this energy transport is 
equal to the increase in energy in a vertica 
air column above the surface from the level 
z to the top of the atmosphere. We did 
estimate the wind in the level 6 m being about 
2 m/sec but no measurements took place. 
However, it is clear from similar cases that 
the wind must have increased with height 
roughly at a logaritmic rate. (The wind 
profile must be fairly complicated in the 
lowest centimeter layer as the speed of the 
water flow is as high as about half a meter per 
second.) The column in which the energy has 
been stored above the level z,, is then not 
vertical but in a twodimensional picture forms 


Uz 

an area between the curves u.f, and ut, + — 
u 

The temperature increase from A (t=o) to 
B (t=t,) in the level z is A T.. Then we have 


fy oo 


90 E u 
- ford Baim foenr.-Ede 
0 Z1 
The data available are not sufficient for a 
direct computation of A, from this equation. 
We therefore make the following simplifica- 
tions. At first we put @ equal to the constant 
value 1.2: 10 g/cm’. We further assume 


that A, is independent of time, that = is 
2 


equal to zero at f=o and that = has a 
FR 


measured value atıx=x;. 
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dz 
time f, in each level has the value measured 
at «=x, in that level. We then get the follow- 
ing estimates of A, at different levels. t, = so 
sec. See Table I. From the measurements at 


C we can see that = in the lowest levels 


rapidly approaches the values at x=x, but 
in higher levels this is not the case. The values 
of A, in the lowest layers must be considered 
more reliable than the values in the upper- 
most levels. 


C 


According to Exner A,= 30° 

4-10: — 
3-4 1 FE 
This is an empirical law valid in certain cases. 
However, if the absolute value of the lapse 


rate is greater than 3.4 : 1074 it is obviously 
not valid. 


Sverdrup, see LETTAU (1938), gives the 
expression 


Zs, 
Ae (1 + ßRi) 


where A‘ is the eddy conductivity at in- 
different stratification, ß is a constant and Ri is 
the Richardson number 
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ENERGY TRANSPORT FROM A WARM RIVER SURFACE INTO COLD AIR $31 
Table I 
RS a EB EIN eee eee 
Height above Profile “NY oa ice 20 
Eur == 
water surface A B C par Je DO dz As 
1,200 cm — 17.9 — 18.0 (— 0.1) 
900 — 18.1 — 18.0 (0.1) 
600 — 17.8 —197.9 — 16.5 0.05 20 0.02°/m 12 
400 0.25 80 0.10°/m 1.0 
300 — 17.6 — 16.6 0.4 150 0.18°/m 1.0 
200 245 0.24°/m 12 
50 415 1.04 °/m 0.34 
30 — 17.0 1.0 445 4.5°/m 0.18 
20 — 16.3 1.7 470 8.5°/m 0.12 
Io — 15.3 27 505 16 °/m 0.072 
4 15.0 3.0 525 0.4°/cm 0.048 
3 — 14.2 — 16.0 3.8 530 0.65°/cm 0.030 
2.5 — 14.7 33 532 1.4 /cm 0.018 
2 — 14.2 3.8 533 2.0°/cm 0.013 
1.5 — 13.5 4.5 534 2.6°/cm 0.011 
I — II.0 7.0 536 5.5°/cm 0.006 
0.5 — 6.5 — 14.0 175 539 11.0°/cm 0.0035 
o + 0.0 18.0 541 
x=x]=100 m at z=6 m corresponds to g 90 
2 ae 
t=t,. We now assume that ae and T at the ome 


ar 5 
dz 


Ri is negative over the warm water surface 
and the absolute value is so large that A, gets 
unlimited. Thus even Sverdrup’s formula does 
not hold. 

Jeffreys, see Brunt (1944), p. 244, has given 
another criterion on the greatest possible 
gradient. If the actual lapse rate is large com- 
pared with the dry adiabatic lapse rate, the 
condition is: 

Of re CL 
dz ~ 4gh*o? 


Here g is the acceleration of gravity and h 
is the thickness of the layer. 
The formula is valid in a homogeneous free 


j DIE 
layer with ie and A constant throughout 


OZ 
the layer. The formula is thus not directly 
applicable when there is a great variation with 


SRE) 
height of the two quantities = and À. How- 


aie 3 
ever mean values of = and A in the layer 


considered have been used here. 
The formula is valid for a layer with free 
upper and lower boundaries. In our case the 
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lower boundary is the water surface and the 
4 

constant ~~ should be replaced by the 


constant I,100. 
Applying the formula for three different 


thicknesses of the bottom layer we obtain: 


Table II 
oT 1,1004?T 
h 4 = “Jr | °Icm 
2 mean °8° | 9 mean ghto? 
2 0.01 7 1,000 
20 0.07 0.8 5 
50 0.08 0.34 182 


From Table II we see that the actual lapse 
rate never reaches the critical value at which 
convection starts, according to Jeffreys formula. 
The discrepancy between the theoretical and 
actual gradients may be due the above described 
limitations of Jeffreys’ formula. 

Studies of cold air outbreaks over warm 
water by G. I. Taylor and others, see BRUNT, 
p- 228 (1944), have given values of A, amount- 
ing to $0—100 cgs in the lowest 300 m deep 
layers whereas values of A, amounting to 
1-3 cgs have been obtained in a similar layer 
when warm air flows over cold water. That 
is to say A is about so times larger in slightly 
unstable conditions than in slightly stable 
conditions. 


A. NYBERG AND L. RAAB 


We have obtained about the same low 
values of A in extremely cold air over a 
very warm water as NyBERG (1939) found at 
extremely stable conditions. Nyberg found in 
his group 6 and 7, which had the same winds 
as in this case, Ayy=29 : 1074 and 37: 1074 cas. 
Here Aj, is equal to 60 : 1074 cgs. The value 
of the molecular heat transport coefficient 


corresponding to A is —=2-10~* cgs and 
c 


p 
this value of A; is reached at about the level 
1 mm over the snow field and at a level less 
than half a mm over the water. 

Our suggestion is that the air over the 
warm water is not in a very unstable state and 
that A, is a function of time, or of x, that is 
to say that convection cells do not start on 
too small a scale. They have to have a con- 
siderable space dimension. An observation of 
a Swedish glider pilot Mr Söderholm in- 
dicates that a length of a warm bay of about 
half a km is needed to get warm bubbles that 
can be used by gliders. 
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Letters to the Editor 


Remarks on »Electron Microscope Study of Cloud and Fog Nuclei» 


Dear Sir, 


We have read with considerable interest the pa- 
per by S. Ogiwara and T. Okita on an »Electron 
Microscope Study of Cloud and Fog Nuclei in the 
August 1952 issue of this Journal, and feel obliged 
to write this letter because of questions concerning 
certain observations and conclusions made by the 
authors. 

A. Because the studies were made on only 75 
fog droplets, we feel that any results are inconclu- 
sive because of the difficulties in obtaining a re- 
presentative sample from so few droplets. 

B. The conclusions on changes in form drawn 
from the electron microphotographs shown in the 
paper are questionable. The authors state that par- 
ticles in Figures 1, 2, 3, 4, 6, and 7 did not change 
their from after electron bombardment and wett- 
ing, while 10, 11, 12, and 13 changed. Inspection 
of the pictures shows no appreciable change in any 
of the particles. The slight differences seen in com- 
paring the microphotographs appear due to a change 
in the focus of the already poorly focussed electron 
microscope and perhaps slight aging of the collodian 
film. These apparent changes appear just as notice- 
able in the »unchanged» pictures. 

C. Because of the growth of a droplet which 
had no visible nucleus in its center, and the fact 
that hygroscopicity changed when external particle 
shape did not; it is felt that the particles investiga- 
ted were not the nuclei concerned with droplet for- 
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mation and the true nucleı were too small to be 
observed with these optical-electron microscope 
techniques. 

D. The investigation of a so-called sulfuric acid 
droplet in the electron microscope and the finding 
that it did not change after electron bombardment 
is improbable. Sulfuric acid, in spite of its low va- 
por pressure, will evaporate quite rapidly at the low 
pressures in the electron microscope. Even if it were 
possible to focus an electron beam on the droplet, 
the increased temperature due to bombardment with 
the electron beam would hasten its evaporation. 

E. The statement that sea salt does not change 
after bombardment with the electron beam in the 
microscope is contrary to our finding, using the 
RCA electron microscope. When we observed sea 
salt crystals in our electron microscope the change 
from a contrasting opaque material to a transpa- 
rent crystal was quite spectacular. 

We feel that these discrepancies between our 
findings and those of the authors should be brought 
before the public so that other interested readers 
can help clarify some of these disputed points. 


Cambridge, March 26th 1953 


S. J. BIRSTEIN 


Geophysics Research Directorate 
Air Force Cambridge Research Centre 
Air Research and Development Command 
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Reply 
Dear Sir, 


We thank Dr Birstein for his discussion of 
our study which was done when we were at 
Sendai. We studied the nature and origin of con- 
densation nuclei and concluded that the greater 
part of cloud and fog nuclei are composed of 
both hygroscopic and non-hygroscopic substances 
which originate in combustion and that sea-spray 
was not a significant source of condensation nuclei. 
It was also shown that the greater part of the 
hygroscopic substances contained in the nuclei 
were not pure sulphuric acid. 


Dr Birstein asked us some questions mainly 
on the techniques of electron microphotographing. 
We should like to make the following statements: 

A. It is clear that in order to obtain any quan- 
titative results, for instance regarding the size 
distribution of nuclei, the number of droplets 
which were observed by us was too small. We 
only studied the nature and origin of the nuclei 
in which case our samples may be of interest. 

B. In fig. (13) filmy substances are observed as 
if they flew out from the upper part of the large, 
well-defined dark particle, though it is not so 
distinct in the reproduction owing to the con- 
traction of the photograph. 

However the change in fig. (11) appears due 
to the poorly focussed electron microscope as Dr 
Birstein states. In the other figures, namely (2), 
(4), and (6), such changes are not perceptible. 
Therefore the original statement should be revised 
as follows. “Most nuclei retained their original 
form after wetting, in some nuclei which retained 
their hygroscopic nature, however, the filmy 
substance changed its form.” 
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C. Only one cloud sample was obtained in 
which no nuclei were found neither by the optical 
nor by the electron microscope. We do not know 
what kind of nucleus the droplet contained. It 
may be such a small volatile particle that it has 
evaporated perfectly by the electron bombardment. 

From the theory and the current opinion on 
condensation, however, it is unreasonable to con- 
sider that most nuclei are too small to be detected 
in the electron microscope. 

D. Before we made our observations we had the 
same opinion as Dr Birstein. However, observa- 
tions yielded a different result. Through observa- 
tion of a fuming sulphuric acid droplet by electron 
microscope we found that the droplet did not 
change when the intensity of bombardment was 
about the same as when microphotographing the 
nuclei. When the intensity became strong or when 
the time of bombardment was long enough, the 
droplet evaporated and left a honey-comb-formed 
image on the collodion film. 

E. In the previous paper we only inferred the 
change of hygroscopicity but not the change of 
transparency of sea-salt. 

We think that Dr Birstein’s finding with regard 
to E is a very interesting and important contribu- 
tion to our knowledge of nuclei in the atmos- 
phere. 

S. OGIWARA 
Institute of Low Temperature Science, Hokkaidö 
University, Sapporo, Japan 
and 
T. Oxira 


Geophysical Institute, Asahigawa Branch, Hokkaidö 
Educational University, Asahigawa, Japan 
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